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ABSTRACT 


The  last  week  of  April  and  first  week  of  May  1988  along  the  central  California  coast 
were  characterized  as  a  period  of  strong  coastal  upwelling  produced  by  moderate  to 
strong  northwesterly  winds  present  throughout  the  period.  A  product  of  this  upwelling 
event  was  the  manifestation  of  a  southward  geostrophic  jet  which  extended  to  a  distance 
of  approximately  50  km  from  the  coast. 

During  the  period  from  08  to  11  May  1988,  a  hydrographic  survey  consisting  of  17 
stations  each  occupied  four  to  five  times  was  conducted  within  the  Monterey  Bay. 
Internal  waves,  with  amplitudes  of  up  to  30  m  were  present  throughout  the  period  and 
effectively  masked  the  mean  signal,  implying  that  averaging  is  essential  to  avoid  aliasing. 
The  CTD  data  were  averaged  to  estimate  the  mean  field  during  this  time  frame.  ADCP 
data,  also  acquired  during  this  period,  were  also  averaged.  The  mean  (low  field  and  dy¬ 
namic  topography  implied  anticyclonic  surface  flow  with  cyclonic  flow  at  200  m  depth, 
both  roughly  centered  over  the  axis  of  the  Monterey  Submarine  Canyon  near  the  mouth 
of  the  Bay.  Flow  at  depth  appeared  to  be  Canyon  "trapped."  Surface  current  speeds 


were  on  the  order  of  10  to  15  cm  with  somewhat  slower  speeds  at  depth  of  approxi¬ 
mately  10  cm  r1.  ADCP  derived  mean  flows  compared  favorably  with  geostrophic  mean 
flows  in  all  areas  except  one,  the  deep  outflow  region  along  the  northern  wall  of  the 
Canyon.  Calculated  standard  errors  were  smaller  than  the  mean  signal  for  both  data 
sets.  The  smaller  standard  errors  calculated  for  the  ADCP  data  implied  that  this  method 
resolved  the  mean  signal  better  than  did  the  geostrophic  calculations. 

Application  of  the  geostrophic  mean  field  to  the  open  boundary  of  a  two  layer, 
primitive  equation  numerical  ocean  model  yielded  flows  similar  to  those  described  above. 

The  inclusion  of  linear  bottom  friction  was  shown  to  be  particularly  important  in  limit¬ 
ing  deep  flow  to  the  confines  of  the  Canyon.  Wind  stress  forcing  experiments  indicated 
that  a  strong  wind  field  may  influence  surface  circulation  in  the  Bay. 

Mean  flow  fields  compared  favorably  with  results  of  Klinck's  (1989)  three-level 
model  of  geostrophic  adjustment  of  flow  over  a  "narrow"  submarine  canyon.  Inter-  n 
actions  between  the  coastal  upwelling  geostrophic  jet  and  the  Monterey  Submarine  ,'41 
Canyon  is  believed  to  have  been  a  major  mechanism  responsible  for  producing  the  ob-  ,;ed 
served  mean  flow.  ltion- 
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I.  INTRODUCTION 


A.  THE  MONTEREY  BAY  AND  VICINITY 

1.  Dimensions 

The  Monterey  Bay,  located  on  the  central  California  coast,  is  basically  orien¬ 
tated  north-south  with  the  city  of  Santa  Cruz  located  on  the  northern  shore,  Moss 
Landing  to  the  east,  and  the  Monterey  Peninsula  situated  on  the  southern  shore.  To  the 
west,  the  bay  is  open  to  the  sea.  A  line  drawn  between  Santa  Cruz  and  the  Monterey 
Peninsula  measures  37  km,  and  Moss  Landing  lies  approximately  20  km  east  of  this  line. 
The  bay  contains  two  bights  (Breaker  and  Broenkow,  1989),  one  situated  along  the 
northern  coast  between  Santa  Cruz  and  Aptos,  the  other  between  the  Monterey  Penin¬ 
sula  and  Fort  Ord  (Figure  1  on  page  2).  The  area  of  the  bay  is  approximately  550  kmJ 
constituting  a  volume  of  approximately  420  km3  (Martin,  1964). 

Monterey  Bay  is  essentially  a  shallow  water  bay  (less  than  100  m  deep)  split  into 
northern  and  southern  shelf  regions  by  the  Monterey  Submarine  Canyon.  This  canyon 
is  the  largest  bathymetric  feature  on  the  west  coast,  comparable  in  shape  and  relief  to 
the  Grand  Canyon  of  the  Colorado  (Shepard,  1966).  The  canyon's  axis  runs  roughly 
east-west  within  the  bay,  originating  near  the  surface  at  the  canyon  head  near  Moss 
Landing  where  it  is  5  km  in  width,  then  widens  to  15  km,  rapidly  deepens  to  1500  to  2000 
m  and  begins  to  meander  past  the  mouth  of  the  bay,  northwest  of  the  Monterey  Penin¬ 
sula.  The  two  largest  tributaries  entering  the  Monterey  Canyon  are  the  Soquel  and 
Carmel  Canyons,  the  former  intersecting  the  axis  near  the  center  of  the  bay  from  the 
northeast  and  the  latter  entering  it  from  the  south  just  seaward  of  the  Monterey  Penin¬ 
sula.  Past  the  entrance  to  the  bay,  the  axis  runs  southwest  for  20  km,  then  heads 
seaward  around  a  large  fan  feature  to  the  north  and  continues  deepening  to  over  3000 
m.  The  relief  of  this  feature  is  truly  dramatic  (Figure  2  on  page  3). 

2.  Currents  of  the  Central  California  Coast 

Because  circulation  within  the  bay  »s  undoubtably  affected  by  the  flow  offshore 
(Breaker  and  Broenkow,  1989),  a  brief  description  of  major  current  systems  in  the  vi¬ 
cinity  of  the  California  coast,  and  their  generation  mechanisms,  is  provided.  The  three 
major  currents  influencing  the  California  coast  arc  the  California  Current  (CC),  the 
California  Undercurrent  (CUC),  and  the  Davidson  Current.  These  names  are  usually 
applied  to  features  existing  in  the  long  term  or  climatological  sense  (Hickey,  1979). 


Figure  1.  Hie  Monterey  Bay  (Source:  Breaker  and  Broenkow,  1989). 
a.  The  California  Current 

'I  he  California  Cui rent  refers  to  the  broad,  diffuse,  cquatorward  wind  driven 
flow  of  the  eastern  boundary  of  the  North  Pacific  Subtropical  anticylonic  gyre.  This 
current  dominates  (in  the  mean)  the  eastern  Pacific  flow  and  extends  south  from  the 
North  Pacific  to  approximately  25’  N,  where  it  turns  to  the  southwest.  This  relatively 
fresh  (Pacific  Subarctic  Watci)  cuncnt  is  usually  strongest  within  a  few  hundred  kilo- 
mctcis  of  the  coast.  The  western  boundary,  of  inteicst  in  this  study,  has  arbitrarily  been 
set  at  1000  km  from  shore.  Mean  velocities  range  fiom  1  to  10  cm  j1,  typical  for  most 
eastern  boundary  currents.  Of  more  interest  to  this  study  aic  the  CUC  and  the 
Davidson  current,  both  located  in  the  near  shore  legions,  typically  cast  of  the  CC 
(Hickey,  1979). 


Figure  2.  Bottom  topography  of  the  Monterey  Bay  Submarine  Canyon  and 
vicinity:  lop  figure  extends  seaward  to  112  km,  bottom  to  60  km  (fig¬ 


ures  derived  from  NOAA/NOS  250  m  resolution  gridded  data). 
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b.  The  California  Undercurrent  and  Davidson  Current 

The  CUC  is  the  mean  poleward  flowing  current  with  mean  core  velocities 
of  10  cm  s-1  that  is  typically  found  along  the  continental  slope  at  depths  ranging  from 
100  to  300  m,  The  Davidson  Current  refers  to  the  northward  flowing  surface  current 
found  along  the  coast  with  velocities  on  the  order  of  10  cm  s~l  and  a  width  of  approxi¬ 
mately  80  km  (McCreary  et  al.,  1987). 

Prevailing  winds  out  of  the  northwest,  generated  by  a  quasi-stationary  at¬ 
mospheric  high  pressure  system  (Figure  3  on  page  5)  over  the  north  central  Pacific,  are 
present  year  round  over  the  California  coast.  The  North  Pacific  High  typically 
strengthens  during  the  summer  months  while  the  thermal  low  pressure  over  the 
California  coast,  particularly  in  the  central  valley,  deepens  as  the  land  warms.  The 
combined  effects  of  these  phenomena  results  in  a  tighter  atmospheric  pressure  gradient 
producing  strong  northwesterly  winds  along  the  California  coast  during  the  summer 
months  (Huyer,  19S3).  Halliwell  and  Allen  (1987)  observed  that  high  pressure  frequently 
builds  northeast  across  California  after  the  passage  of  an  atmospheric  cyclone  and  its 
associated  cold  front  generating  periods  of  strong  northwest  winds  and  upwelling  events 
along  the  coast.  They  also  noted  that  wind  maxima  tended  to  propagate  northward 
along  the  coast.  Chelton  (1984)  noted  that  the  spatial  structure  of  the  wind  stress  field 
is  highly  coherent  in  the  alongshore  direction  and  that  the  magnitude  varies  seasonally 
with  the  maximum  intensity  in  May  (Figure  4  on  page  7).  Since  the  strongest  compo¬ 
nents  are  found  approximately  200  km  offshore,  a  positive  wind  stress  curl  (curl  r)  is 
present  year  round  in  the  nearshore  regions  (Nelson,  1977).  These  winds,  through  the 
Ekman  transport,  produce  near  shore  upwelling  resulting  in  dynamic  height  fields  which 
slope  downward  towards  the  coast.  The  result  is  a  geostrophically  balanced  surface  jet 
which  flows  in  the  direction  of  the  wind  (i.e.  equatorward)  and  a  poleward  undercurrent 
(the  CUC)  at  approximately  100  to  300  m  in  depth.  Hickey  (1979)  observed  that  this 
surface  flow  is  strongest  along  the  Pacific  northwest  coast  during  the  spring  with  maxi¬ 
mum  flow  typically  located  within  the  internal  Rossby  radius  of  deformation,  between 
5  and  25  km  from  the  coast,  and  noted  that  variations  in  width,  location,  and  strength 
can  occur  quite  rapidly  (on  the  order  of  days)  in  response  to  local  changes  in  wind  stress 
(t).  She  also  observed  that  during  periods  of  weak  northerly  or  southerly  winds,  this 
surface  jet  either  weakens  or  disappears  altogether.  Hickey  (1979)  finally  noted  that  the 
dynamics  of  the  nearshore  southward  flow  observed  off  the  California  coast  are  likely 
similar  to  the  flows  noted  along  the  Pacific  northwest  coast. 
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Figure  3.  Long-term  mean  sea  level  atmospheric  pressure  over  the  no.th 
Pacific:  Top  figure  is  for  January,  bottom  figure  is  for  July  (lourcc: 
Iluyci,  1983). 


During  the  summer,  observed  current1:  along  the  coast  arc  in  general 
agreement  with  i'kman  transpoit  theory  but  during  the  winter  poleward  flow  (the 
Davidson  Current)  occurs  along  the  coast,  against  the  surface  wind  sticss  (McCreary  ct 
al.,  1987).  McCreaiy  et  al.  (1987)  cite  four  possible  mechanisms,  along  with  pros  and 
cons  of  each,  which  might  explain  the  appearance  of  the  Davidson  cuncnt,  namely; 

•  Poleward  propagation  of  coastally  trapped  waves. 

•  Relaxation  of  local  winds.  Sustained  northwesterly  winds  produce  a  poleward 
pressure  gradient.  When  these  winds  relax,  a  poleward  cuircnt  will  be  driven  by 
this  pressure  gradient.  'Ibis  theory  is  consistent  t’  c  appearance  of  the  Davidson 
Current  during  the  winter,  when  mean  equatorwaid  suifacc  wind  stress  decreases. 
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•  Positive  wind  stress  curl.  Sverdrup  (1947)  surmised  that  the  vertically  integrated 
transport  in  the  meridional  direction  (My)  is  proportional  to  the  vertical  component 
of  the  curl  of  the  wind  stress,  or: 
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where  fi  =  djjdy  is  the  meridional  change  in  the  Coriolis  parameter  (/),  ^  is  the 
stream  function  ( diff/dx  >  0  implies  northward  flow),  and  tx  and  ry  are  the  x  and  y 
components  of  the  wind  stress  respectively.  Thus,  positive  curl  r  should  drive  a 
poleward  current,  according  to  Sverdrup  theory.  The  problem  with  this  mechanism 
however,  is  the  fact  the  Davidson  Current  is  usually  found  during  the  winter 
months,  a  time  period  when  curl  t  is  typically  the  weakest. 

•  Thermohaline  forcing.  The  north-south  density  gradient  between  subtropical  and 
polar  regions  creates  lower  dynamic  heights  over  the  poles  and  relatively  higher 
dynamic  heights  in  the  equatorial  region.  This  situation  drives  a  mid-ocean  east¬ 
ward  geostrophic  current  which  should  turn  cyclonically  (poleward)  when  it  en¬ 
counters  an  eastern  boundary  such  as  the  California  coast.  They  argue  that  this 
mechanism  is  likely  to  be  responsible  for  driving  the  Leeuwin  Current  off  western 
Australia;  however,  since  wintertime  sea-surface  topography  off  tilt  California 
coast  differs  significantly  from  that  observed  off  the  west  Australian  coast,  (i.e. 
dynamic  heights  bend  toward  the  equator  vice  the  pole  off  California),  this  mech¬ 
anism  appears  to  play  a  role  which  is  subordinate  to  wind  forcing  in  this  area. 


Using  both  wind  stress  curl  and  curl  free  wind  stress  as  inputs  into  a  coastal  model, 
McCreary,  et  al.  (1987),  concluded  that  positive  wind  stress  curl  was  the  most  likely 
mechanism  responsible  for  driving  the  Davidson  Current  during  the  winter.  During  the 
summer  months  the  along  shore  wind  stress  (t,)  is  strong  enough  to  negate  the  effects 
of  Sverdrup  transport  and  therefore  dominates  the  surface  flow'  field.  Hickey  (1979) 
states  that  oceanic  responses  to  sign  changes  in  curl  t  w'ould  be  slow,  i.e.  on  the  order 
of  two  to  three  months. 

Chelton  (1984),  through  the  averaging  of  23  years  of  hydrographic  data  ac¬ 
quired  through  the  CalCOFI  program,  found  the  surface  flow  field  to  be  closely  coupled 
to  the  w'ind  field.  Analysis  of  the  derived  dynamic  height  fields  at  the  surface  and  at  150 
dbar  relative  to  500  dbar  at  Monterey  (station  67.57)  showed  surface  flow  to  be  south¬ 
ward  from  February  to  June  and  from  August  to  November  with  polew'ard  flow  from 
November  to  February  and  again  in  June  to  July.  Subsurface  flow  (CUC)  towards  the 
north  w’as  present  year  round  except  for  weak  southward  flow  in  April  and  May.  Further 
south,  at  Point  Conception,  no  reversal  in  the  CUC  w'as  noted,  and  thus  Chelton  con¬ 
cluded  that  subsurface  convergence  must  be  taking  place  somewhere  between  Pt  Con- 
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Figure  4.  Long  term  average  wind  stress  (top)  and  wind  stress  curl  (bottom)  along 
the  California  coast:  Averages  were  obtained  from  100  years  of  ship 
observations.  The  shaded  region  corresponds  to  positive  wind  stress  curl 
(Source:  Nelson,  1977). 

ception  and  Pt.  Sur,  although  the  CalCOFI  grid  spacing  was  too  coarse  to  resolve  this 
convergence.  In  general,  Chelton  (1984)  concludes  that  his  findings  support  the  general 
description  of  an  equatorward  CC,  a  poleward  flowing  CUC,  and  a  coastal  surface  cur¬ 
rent  whose  direction  of  flow  is  seasonally  dependent. 

In  summary,  since  coastal  circulation  undoubtably  influences,  and  possibly 
drives,  flow  within  the  Monterey  Bay,  full  knowledge  of  the  near  shore  processes  which 
take  place  along  the  central  California  coast  is  crucial  prior  to  attempting  to  understand 
(low  patterns  within  the  bay  itself.  For  instance,  surface  flow  within  the  Monterey  Bay 
during  winter,  when  the  Davidson  Current  is  present,  more  than  likely  differs  from  flow 
patterns  within  the  bay  during  spring  upwclling  periods  when  the  southward  geostrophic 
jet  is  present.  Also,  deep  flow  within  the  Monterey  Bay  Submarine  Canyon  is  probably 
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influenced  to  some  extent  by  the  CUC,  although  the  effects  of  this  current  on  circulation 
within  the  bay  does  not  appear  to  be  very  well  understood  at  present. 

B.  SYNOPSIS  OF  PREVIOUS  MONTEREY  BAY  STUDIES 
1.  Direct  Observations 

Numerous  studies  have  been  conducted  to  ascertain  the  characteristics  of  the 
circulation  within  Monterey  Bay  utilizing  a  wide  variety  of  methods  including  current 
meters,  drift  bottles,  drogues,  hydrographic  data,  thermistors,  satellite  imagery,  aircraft 
observations,  and  sediment  and  nutrient  distributions.  Breaker  and  Broenkow  (1989) 
provide  an  excellent  in-depth  summary  of  past  research  conducted  in  the  bay,  including 
comments  relating  these  works  to  their  own  findings.  Highlights  of  some  of  these  major 
studies  are  provided. 

a.  Early  Studies 

One  of  the  earliest  studies  of  flow  patterns  within  Monterey  Bay  was  con¬ 
ducted  by  Bigelow  and  Leslie  (1930)  using  hydrographic  data  obtained  throughout  the 
bay  over  the  period  30  June  to  24  July  1928.  They  observed  higher  sea  surface  temper¬ 
atures  over  the  shelf  areas  of  the  bay  and  somewhat  depressed  temperatures  over  the 
canyon  implying  that  local  heating  played  a  major  role  during  the  summer  months. 
Analysis  of  the  cross  sectional  temperature  field  also  led  them  to  believe  that  upwelling 
from  within  the  canyon  occurred  and  subsequently  spread  out  into  the  shelf  regions  of 
the  bay.  Dynamic  height  analysis  (0  dbar  referenced  to  500  dbar)  revealed  the  presence 
of  high  dynamic  heights  near  the  center  of  the  canyon  implying  anticyclonic  flow.  They 
concluded  that  circulation  within  Monterey  Bay  varies  greatly  and  is  dramatically  influ¬ 
enced  by  week  to  week  changes  in  small  scale  local  events  (heating,  land-sea  breeze, 
winds,  canyon  upwelling,  etc.). 

Skogsberg  (1936),  through  analysis  of  hydrographic  data  acquired  between 
1929  and  1933  in  the  southern  portion  of  Monterey  Bay,  concluded  that  three  separate 
cycles  occurred  annually.  The  first  he  described  as  a  "cold  water"  or  upwelling  phase 
which  lasted  from  February  through  November,  a  "warm  water  phase"  attributed  to  light 
winds  and  onshore  movement  of  warm  oceanic  waters  from  August  through  November, 
and  finally,  a  "low  thermal  gradient  phase"  which  lasted  from  December  to  February  and 
coincided  with  the  manifestation  of  the  northward  flowing  Davidson  Current.  He  em¬ 
phasized  the  highly  variable  nature  of  flow  vithin  the  bay  stating  instances  where  coastal 
flow  was  northward  through  the  bay  and  other  instances  where  it  was  southward. 
Through  these  thermal  studies,  Skogsberg  finally  noted  that  temperature  trends  at  lower 
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levels  were  quite  often  opposite  to  those  encountered  at  shallower  depths  suggesting  that 
flows  may,  at  times,  be  independent  and  opposite  at  separate  levels  within  the  bay. 
b.  Surface  Current  Studies 

Utilizing  nutrient  displacement  and  mass  fields,  Lasley  (1977)  concluded 
that  surface  flow  within  the  bay  is  predominantly  to  the  north,  frequently  entering  the 
bay  around  Pt.  Pinos  and  exiting  near  Pt.  Santa  Cruz.  Utilizing  property  distributions 
assimilated  over  a  27-month  time  span,  Smethie  (1973),  and  Broenkow  and  Smethie 
(1978),  also  inferred  surface  flow  within  the  bay  to  be  northward  the  majority  of  the 
time,  in  agreement  with  Lasley.  They  expected  northward  flow  during  the  fall  and  winter; 
however,  in  light  of  the  northwesterly  winds  present  during  the  spring  and  summer, 
prevailing  northward  surface  flow  within  the  Monterey  Bay  was  contrary  to  their  ex¬ 
pectations  during  this  time  frame.  Current  meter  data  acquired  by  Engineering  Science, 
Inc.  (ESI,  1978)  in  the  area  southwest  of  the  Salinas  River  at  depths  of  9  and  30  m  over 
the  period  from  January  1976  to  January'  1977,  indicated  a  northward  flow  65%  of  the 
time  and  southward  flow  35%  of  the  time.  Reise  (1973),  using  drift-bottles  over  a  14 
month  period  in  1963  and  1964  in  the  southern  bight,  and  M corny  (1973),  employing 
drogues  over  a  two  day  period  in  August  1972,  both  determined  surface  flow  within  the 
bay  to  be  generally  cyclonic  in  nature  (Figure  5  on  page  10). 

Pirie  and  Stellar  (1974)  inferred  surface  flow  within  the  Monterey  Bay  and 
vicinity  through  sediment  distribution,  obtained  from  visual  satellite  imagery-  and  aircraft 
observations.  Their  flow  patterns  suggested  large  current  variability  but  seasonal  trends 
indicated  northward  flow  along  the  coast  from  October  through  February-,  with  south¬ 
ward  flow  predominating  from  March  through  July.  Northward  flow  appears  to  corre¬ 
late  well  with  the  classic  Davidson  Current  regime  from  October  through  February,  with 
southward  flow  occurring  during  the  typical  coastal  upwclling  period.  Breaker  and 
Broenkow  (1989),  using  AVHRR  imagery  taken  during  the  upwelling  spring  transition 
in  March  1980  observed  cold,  upweiled  water  extending  along  the  entire  central 
California  coast  from  a  strong  thermal  front  at  Pt.  Pinos.  During  this  event,  sea  surface 
temperatures  at  Granite  Canyon  (located  between  the  M  ontercy  Peninsula  and  Pt.  Sur, 
see  Figure  1  on  page  2)  dropped  4°C  over  a  six  day  period  (Breaker  and  Mooers,  1986). 
Coincident  with  this  event,  current  direction  was  observed  to  reverse  from  north  to  south 
in  the  vicinity  of  Pt.  Sur  (Wickham  et  ah,  1987).  Other  satellite  imagery  reviewed  by 
Breaker  and  Broenkow  (1989)  suggest  instances  of  warm  water  intrusions  entering  the 
bay  from  the  open  sea.  They  concluded  that  during  strong  upwelling  events,  cold  water 
was  often  adveacd  across  the  entrance  of  the  bay,  then  south  towards  Pt.  Sur,  forming 
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Tigure  5.  Surface  densib- analysis  and  drogue  trajectories  (21-22  Aug  1972):  Note 
inferred  cyclon. .  circulation  (Source:  Moomy,  1973). 

a  continuous  cold  frontal  feature.  Another  common  characteristic  found  in  satellite 
images  was  that  generally  warmer  waters  were  found  inside  the  bay  providing  support 
to  Bigelow  and  Leslie's  (1930)  findings  that  local  heating  plays  a  large  role  in  modifying 
water  masses  within  the  bay  and  coastal  upwclling  occurs  on  an  infrequent  basis  within 
the  bay  itself  (Breaker  and  Broenkow.  1989). 
c.  Flow  at  Mid-Depths 

Flow  at  intermediate  depths  (30  to  120  m)  was  investigated  by  Lammers 
(1971),  using  data  acquired  from  Skogsberg,  1929-1933;  Bolin  ct  al.,  1951-1955;  and  the 
CalCOFI  program,  1954-1967  to  construct  a  gridded  temperature  field  over  Monterey 
Bay.  Gcostrophic  flow  patterns  wcic  then  implied  tluough  analysis  of  monthly-mean 
isotherm  topographies  (Figure  6  on  page  1 1).  '1  he  study  concluded  that,  in  the  mean,  a 
warm  core  antic} clonic  feature  within  the  thcimoclinc  persisted  over  the  center  of  the 
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rigure  6.  Monthly-mean  isotherm  topographies  (in):  Figures  were  obtained  from 
30  ycats  of  data  and  show  a  highly  developed  thermal  high  situated  over 
central  Montcicy  Bay  (Source:  Adapted  from  hammers,  1971). 


bay  and  tended  to  intensify  duting  the  winter.  This  waim  core  feature  is  in  contrast  to 
Bigelow  and  Leslie's  (1930)  findings  of  slightly  depressed  temperatures  over  the  canyon 
during  a  one  week  period  in  summer  but  is  in  agreement  with  hammers"  inferred 
anticyclonic  circulation.  One  possible  mechanism  thought  to  be  responsible  for  this 
feature  is  momentum  and  heat  transfer  from  the  northward  flowing  CUC  (Breaker  and 
Broenkow,  19S9). 

d.  Deep  Flow  Studies 

A  series  of  deep  flow  studies  were  undertaken  from  1965  to  1975  by  students 
at  the  Naval  Postgiaduatc  School  (NFS).  These  studies  weic  primarily  conducted  near 
the  canyon  head  at  depths  ranging  from  91  to  485  m  (T  able  1  on  page  12)  and  were 
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Table  1.  SUMMARY  OF  NPS  NEAR-BOTTOM  CANYON  CURRENT 


STUDIES. 


Reference 

Bottom 

Depth(s) 

Duration 

Comments 

Gatje  and 
Pizinger,  1965 

130  m 

77  observa¬ 
tions  on  nine 
different  days 

! 

Currents  tended  to  follow  the 
canyon  axis.  Downcanyon  flow 
was  associated  with  incoming 
tide  and  upcanyon  flow  with 
outgoing  tide. 

Dooley,  1968 

104- 165m 

5-162  h 

Flows  were  predominantly  in 
direction  of  canyon  axis.  Re¬ 
versals  in  flow  direction  showed 
a  strong  semidiurnal  compo¬ 
nent.  Sudden  increases  or  bursts 
in  current  speed  also  were 
noted.  Spectral  analysis  re¬ 
vealed  6  and  4  h  peaks  in  addi¬ 
tion  to  semidiurnal  component. 

Xjus,  1968 

150-202  m 

168  h 

Current  direction  was  generally 
along  the  canyon  axis.  Re¬ 
versals  in  current  direction  oc¬ 
curred  primarily  at  high  and  low 
tide.  Spectral  peaks  occurred 
at  12.5,  6.  and  4  h. 

Caster,  1969 

91-366  m 

3-171  h 

1 

Cross-canvon  How  was  fre¬ 
quently  observed.  Flow  at  an 
adjacent  shelf  location  showed 
"eaker  tidal  oscillations  and 
’enerallv  slower  speeds. 

Hollister,  1975 

485  m 

262-371  h 

Flow  at  30  m  above  bottom  was 
predominantly  along  canyon 
axis  and  apparently  tidally 
driven.  Strong  cross-canyon 
flow  was  observed  at  60  m 
above  bottom  and  was  poorly 
correlated  with  flow  at  30  m. 

Source:  Adapted  from  Breaker  and  Broenkow,  1989. 


primarily  of  short  duration  ranging  from  a  day  (Gatje  and  Pizinger,  1965)  to  approxi¬ 
mately  two  weeks  (Hollister,  1975)..  Most  current  meter  depths  ranged  from  5  to  16  m 
above  the  bottom  (Hollister's  study  used  meters  placed  30  and  60  m  above  the  bottom). 
Data  acquired  near  the  canyon  head  revealed  predominantly  along-axis  flow  with 
semidiurnal  periods.  Castor  (1969)  and  Hollister  (1975),  on  the  other  hand,  observed 
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frequent  cross-canyon  flow  in  the  deeper,  wider  areas  of  the  canyon.  The  spectral  ana¬ 
lyses  done  by  Dooley  (1968)  and  Njus  (1968)  indicated  super-tidal  periods  at  6  and  4  h 
in  addition  to  the  expected  semidiurnal  period.  Gatje  and  Pizinger  (1965)  correlated 
downcanyon  flow  with  the  surface  flood  tide  and  upcanyon  flow  with  the  ebbing  tide. 

Shepard  (1979)  conducted  extensive  research  on  flows  in  numerous  subma¬ 
rine  canyons,  including  the  Monterey  Canyon.  His  work  in  the  Monterey  Submarine 
Canyon  was  at  much  greater  depths  (up  to  1445  m)  than  those  undertaken  by  NPS 
(Figure  7  on  page  14).  Of  the  nine  records  obtained,  five  showed  virtually  no  directional 
correlation  with  the  canyon  axis,  unique  among  all  of  the  canyons  investigated.  One 
suggestion  for  the  cause  of  this  lack  of  directional  control  was  related  to  the  meandering 
of  the  axis  along  with  flows  into  and  out  of  canyon  tributaries.  Another  interesting  as¬ 
pect  noted  by  Shepard  (1979)  was  that,  unlike  other  California  stations,  the  net  flow  was 
upcanyon  in  6  out  of  the  9  records.  Shepard  stated  that  current  direction  and  surface 
tides  were  not  always  coherent  (Figure  8  on  page  15),  at  times  in  phase  with  the  tidal 
cycle,  at  other  times  out  of  phase.  By  observing  discrepancies  noted  when  time-velocity 
plots  of  different  stations  were  overlayed  (in  particular  stations  34  (at  357  m)  and  35  (at 
384  m)  and  stations  58  (at  1061  m)  and  59  (at  1445  m)  (see  Figure  7  on  page  14)),  then 
shifting  the  profiles  in  time  to  match  the  observed  up  and  down-canyon  flows,  Shepard 
concluded  that  internal  wave  patterns  were  propagating  upcanyon  between  stations  at 
shallower  depths  but  downcanyon  between  stations  at  greater  depths. 
e.  Internal  Waves 

Almost  every  study  done  in  Monterey  Bay  cites  evidence  indicative  of 
internal  waves.  Shea  and  Broenkow  (1982)  investigated  the  occurrence  of  semidiurnal 
internal  waves  using  hydrographic  and  thermistor  data  and  resolved  internal  waves  with 
heights  of  50  to  120  m  (Figure  9  on  page  16).  Spectral  and  cross-spectral  analysis  of  7 
days  of  surface  tidal  heights  and  thermistor  data  acquired  near  the  canyon  head  sug¬ 
gested  that  these  waves  have  a  semidiurnal  periodicity  but  lag  the  tides  by  about  7  h  at 
the  semidiurnal  period.  They  cited  the  Monterey  Submarine  Canyon  topography  as  the 
probable  generation  mechanism.  As  the  semidiurnal  long  wave  field  approaches  the 
shoaler  and  narrower  canyon  head,  the  period  remains  constant  and  the  amplitudes  in¬ 
crease  as  the  wave  energy  becomes  focused,  as  small  amplitude  wave  theory  predicts,  and 
assumes  the  characteristics  of  an  internal  tidal  bore.  At  high  internal  tides  they  found 
evidence  that  denser  water  from  below  the  main  thermocline  within  the  canyon  moved 
up  and  onto  the  shelves  and  then  became  pinched  off  as  the  internal  tides  fell 
(Figure  10  on  page  17).  Baroclinic  tidal  currents  were  also  investigated  by  Broenkow 
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stations  34,  35,  58,  and  59  were  used  to  imply  internal  wave  propagation 
(Source:  Shepard,  1979). 


and  McKain  (1972)  who  found  that  tidal  amplitudes  generally  decreased  seaward  from 
the  canyon  head  and  noted  that  the  baroclinic  and  the  barotropic  tides  were  generally 
not  in  phase. 

2.  Numerical  Studies  of  Monterey  Bay 
a.  Tidal  A  [odds 

Schomakcr  (1933)  used  an  implicit,  two-dimensional  (1  km  grid)  program 
with  realistic  topography  to  model  homogenous  barotropic  tidal  forcing  in  the  bay. 
Resulting  height  fields  showed  the  clear  progression  of  the  tides  into  the  bay;  however, 
the  inferred  currents  into  and  out  of  the  bay  associated  with  these  heights  were  relatively 
weak.  Consequently,  other  effects,  such  as  wind  forcing  and  ofl'shoie  currents,  were 
determined  to  dominate  the  flow  field,  rather  than  barotropic  tidal  influences. 
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Figure  8.  Along-axis  current  \el«city  vs.  time  in  the  Monterey  Canyon:  Figures 
show  the  incohctcnce  between  tidal  height  (smooth  line)  and  current 
velocities  (shaded).  Current  meter  locations  arc  on  the  previous  figure 
(Source:  Shepard,  1979). 

b.  Mean  Flow  Models 

Using  a  one-layer  baiotiopic  numerical  cavity  flow  model,  with  boundaries 
representative  of  the  bay,  Gaicia  (1971)  attempted  the  first  numerical  study  of  the  bay. 
The  model  was  forced  by  momentum  transfer  (frictional  shear)  from  northward  flow 
adjacent  (west  of)  to  the  cavity  and  included  the  effects  of  bottom  topography 
(smoothed  version  of  actual  topography),  friction,  and  the  Coriolis  force.  Model  out¬ 
puts  were  in  the  form  of  mass  transport  stream  functions.  The  results  of  this  model  led 
to  the  conclusion  that  the  presence  of  a  submarine  canyon  can  induce  a  closed  circu¬ 
lation  within  the  bay  and  that  the  ^-effect  is  negligible  (Figure  11  on  page  18).  He  also 
noted  that  the  effects  of  bottom  friction  on  the  flow  were  negligible  within  the  canyon 
itself  but  significantly  influenced  flow  in  the  shallower  regions.  The  most  important  re¬ 
sult  obtained  fiom  this  model  was  the  fact  that  the  rotational  nature  of  flow  within  the 
bay  was  changed  by  suppressing  or  adding  the  bottom  topography  of  the  Monterey 
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Hgure  9.  Time  series  of  temperature  distributions  near  the  canyon  head:  Data 
were  collected  1  km  (top)  and  4  km  (bottom)  west  of  the  canyon  head 
and  imply  internal  wave  heights  or  50  to  120  m.  Isotherms  in  'C 
(Source:  Shea  and  Brocnkow,  1982). 


Submarine  Canyon.  For  cavity  flow  withou'  canyon  topography,  the  northward  mass 
tiansport  stream  simply  entered  the  bay  from  the  south  and  exited  to  the  north  resulting 
in  cyclonically  curved  (low  stemming  from  Coriolis  foicc  considerations.  When  canyon 
topography  was  added,  momentum  transfer,  through  shear  stress  from  the  north  flow 
adjacent  to  the  cavity,  drove  an  anticyclonic  gyre  which  effectively  filled  the  Monterey 
Bay  (Figure  II  on  page  18). 

Using  output  from  a  two  layer,  semi-implicit,  primitive  equation  model  and 
forcing  across  a  north-south  oriented  boundary,  Bruner  (1988).  showed  that  the  upper 
layer  was  insensitive  to  bottom  topography  and  to  lower  layer  boundary  forcing  for  ve¬ 
locities  less  than  5  cm  .s  '.  When  lower  layer  flow  was  increased  above  5  cm  s‘\  the 
upper  layer  basically  mirrored  the  lower.  The  lower  laver  was,  of  course,  highly  sensitive 
to  bottom  topogiaphy.  Also,  when  the  inflow  width  was  decreased  along  the  northern 
boundary,  and  southern  boundary  outflow  was  increased,  the  surface  flow  changed  dra- 
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Figure  10.  Conceptual  model  (top)  of  temperature  distributions  (bottom)  near  the 
edge  of  the  Monterey  Canyon:  Top  figures  show  the  simple  2-level 
model  used  by  Shea  and  Brocnkow  (1982)  to  describe  how  denser  water 
from  within  the  canyon  is  lifted  out  of  the  canyon  at  high  internal  tide 
(a),  then  pinched  off  onto  the  shelves  during  low  internal  tide  (b). 
Bottom  figures  show  actual  temperature  cross  sections  along  the  shelf 
during  high  internal  tides  (a)  and  low  internal  tides  (b)  (Source:  Shea 
and  Brocnkow,  1982). 
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Figure  11.  Garcia's  one  layer  cavity  model:  Top  figure  represents  model  output 
with  no  canyon  topography;  bottom  figure  shows  output  with 
topographic  eficcts  included.  Contours  are  volume  transport  stream 
functions  with  velocity  arrows  overlaid  (Source:  Garcia,  1971). 

matically  from  cyclonic  over  the  inner  bay  to  totally  anticyclonic  flow  while  the  lower 
layer  flow  was  inhibited  from  entering  the  bay  and  turned,  instead,  almost  immediately 
to  the  south  of  the  bay  where  ii  exited  (compare  Figure  12  on  page  19  and  Figure  13 
on  page  20).  Model  outputs  also  showed  that  flow  within  the  bay  was  highly  sensitive 
to  vertical  shear.  When  no  shear  was  present  between  the  layers  (barotropic  case)  the 
bottom  layer  had  a  great  effect  on  the  surface  layer.  However,  when  vertical  shear  was 
included  (baroclinic  case)  the  layers  were  essentially  decoupled.  Finally,  Bruner  noted 
that,  when  inflow  was  forced  over  the  southern  boundary  simulating  inflow  of  the 
Davidson  Current,  model  runs  depicted  cyclonic  circulation  within  the  bay  (Figure  14 


Figure  12.  Bruner's  two-layer,  baroclinic  model  output  #1:  Surface  flow  (top)  and 
lower  layer  flow  (bottom)  with  forced  inflow  on  the  northern  boundary 
and  outflow  over  the  southern  boundary.,  North  is  to  the  left.  Isotachs 
arc  overlaid  with  velocity  vectors  (Source:  Bruner,  1988). 

on  page  21).  These  results  contradict  Garcia's  barotropic  model  outputs,  which  pre¬ 
dicted  anticyclonic  circulation  within  the  confines  of  the  bay  from  an  adjacent  northward 
flow.  Boundary  forcing  in  both  Garcia's  and  Bruner's  models  was  purely  hypothetical, 
intended  solely  to  test  model  responses  (and  the  physics  involved)  to  various  flow  sce¬ 
narios. 

Two  final  model  studies,  although  not  conducted  exclusively  for  Monterey 
Bay,  arc  mentioned  here  because  of  their  obvious  relevance.  Both  models,  the  first 


Figure  13.  Bruner's  *«o  layer,  baroclinic  model  output  #2:  'lhe  same  as  in 
Figure  12  on  page  19  except  that  inflow  area  was  constricted  and  out¬ 
flow  area  was  widened.  North  is  to  the  left  (Source:  Bruner,  1988). 

solved  analytically  and  the  second  numerically,  attempt  to  represent  the  adjustment 
process  that  takes  place  when  a  gcostrophically  balanced  flow  over  a  shelf  encounters  a 
submarine  canyon,  and  arc  from  Klinck  (1988  and  1989). 

In  the  first  model.  Klinck  (1988)  initially  solved  the  shallow  water  equations 
by  the  Laplace  Transform  method  for  a  two  layer  fluid  on  an  /-plane,  and  then  applied 
the  boundary  conditions  of  a  flat  bottom  ocean  intciscctcd  by  an  infinitely  long  rectan¬ 
gular  trench.  He  then  applied  the  initial  condition  of  an  arbitrarily  wide  flow  in  the 
upper  layer  perpendicular  to  the  canyon  with  no  initial  flow  within  the  canyon  itself. 


Figure  14.  Bruner's  two-layer,  baroclinic  model  output  #3:  Surface  flow  (top)  and 
lower  layer  flow  (bottom)  with  forced  inflow  on  the  southern  boundary, 
outflow  over  the  northern  boundary.  North  is  to  the  left  (Source: 
Bruner,  1988). 

The  solution  to  this  problem  yielded  standing  wave  motion  over  the  canyon  with  waves 
radiating  from  cither  side.  In  the  absence  of  friction,  these  waves  would  continue  to 
oscillate  indefinitely  as  standing  gravity  modes.  Results  also  showed  that,  in  the  steady 
state,  there  was  no  net  flow  across  the  canyon,  implying  that  submarine  canyons  may 
inhibit  barotropic  gcostrophie  flow.  For  shelf  currents  on  the  order  of  the  external 
Rossby  radius  of  deformation  ( R rf„),  these  oscillations  have  periods  on  the  order  of  0.1 
to  1.0  times  the  local  inertial  ficqucncy.  For  "narrow"  canyons  (defined  in  this  case  as 


a  canyon  whose  width  is  less  than  one  half  of  the  smaller  of  either  the  width  of  the  cur¬ 
rent  or  Rd0),  Klinck  showed  that  these  oscillations  would  have  periods  on  the  order  of 
0.2  times  the  inertia  period.  He  notes  that  these  frequency  ranges  are  observed  in  most 
canyons.  Breaker  and  Broenkow  (1989),  citing  this  work,  calculated  a  period  of  4  h  for 
Monterey  Bay  and  noted  the  close  agreement  with  Dooley  (1968)  and  Njus'  (1968)  ob¬ 
served  4  h  periods.  They  concluded  that  both  inertial  and  non-linear  effects  associated 
with  the  interactions  of  the  internal  baroclinic  tides  and  steep  topography  are  both 
possibly  at  work  in  generating  supertidal  oscillations  in  Monterey  Bay. 

Klinck  (1989)  conducted  further  experiments  on  the  geostrophic  adjustment 
process,  this  time  utilizing  a  three  layer  numerical  model  (Figure  15  on  page  24)  where 
both  the  internal  and  external  modes  can  act  (the  two  top  layers  modeled  shelf  flow,  the 
third  layer  modeled  flow'  within  the  canyon).  As  he  states,  the  importance  of  using  a 
three  layer,  stratified  fluid  is  the  fact  that  pressure  gradients  at  the  bottom  of  the  upper 
layers  will  now  influence  the  fluid  within  the  canyon  and  that, ". .  .density  redistribution 
can  cancel  the  effect  of  a  surface  slope  such  that  there  is  a  steady  geostrophic  flow,"  re¬ 
sults  which  were  not  possible  in  his  original  model.  In  addition  to  the  coastal  current 
width,  canyon  width,  and  RM,  the  first  internal  Rossby  radius  of  deformation(Rdl)  now 
becomes  another  important  controlling  factor  in  the  geostrophic  adjustment  process. 
Klinck  initialized  this  model  in  much  the  same  w'ay  as  his  barotropic  model,  i.e. 
geostrophicallv  balanced,  steady,  barotropic  10  cm  s~}  flow  over  an  infinitely  wide  flat 
bottom  in  the  positive  x  direction  (this  would  be  the  same  as  looking  shoreward  at  a 
cross  section  of  Monterey  Bay  from  sea  with  an  initial  shelf  current  flowing  equatorward 
aiong  the  coast)  which  crosses  a  rectangular  canyon  at  right  angles.  Flow  in  the  canyon 
was  initially  at  rest.  A  50  km  wide  current  was  used  and  runs  were  made  using  varying 
canyon  widths  including  200  and  10  km.  In  the  200  km  wide  case,  the  induced  disturb¬ 
ance  at  the  canyon  walls  are  independent  of  each  other  (Figure  15  on  page  24).  Flow 
over  the  shelf  decays  away  from  the  wall  on  the  order  of  the  current  width,  while  flow 
within  the  canyon  decays  away  from  the  wall  on  the  order  of  Rdl  for  the  internal  mode 
and  on  the  order  of  the  current  width  for  the  external  mode.  Also,  the  largest  density 
gradients  are  located  within  the  canyon  which,  as  Klinck  points  out,  indicates  that  the 
strongest  currents  are  along  the  walls  of  the  canyon  and  are  bottom  trapped.  He  ex¬ 
plains  these  results  in  terms  of  vortex  stretching,  stating; 
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During  adjustment,  the  water  in  the  canyon  is  pushed  down  the  initial  pressure 
gradient  which  causes  upwelling  at  places  where  the  surface  elevation  in  the  initial 
state  was  lowest  and  downwclling  where  it  was  highest.  Upwelling  stretches  the 
vortex  tubes  in  the  canyon  creating  cyclonic  vorticitv  at  the  edges  of  the  canyon. 
The  mass  redistribution  reduces  the  slope  of  the  free  surface  creating  cyclonic 
(anticyclonic)  vorticity  in  level  1  on  the  shelf  side  of  the  canyon  wall  in  places  where 
the  free  surface  rises  (falls).  The  isopycnal  surface  over  the  canyon  near  the  wall 
rises  sharply  compressing  vortex  tubes  creating  anticyclonic  vorticity  right  at  the 
wall. 

In  the  narrow  canyon  case  (10  km  wide),  the  external  mode  is  inhibited  because  the 
width  is  now  smaller  than  any  of  the  decay  scales  (current  width,  Rd0,  or  Rdl).  Isopycnals 
within  the  canyon  create  vortex  stretching  (positive  vorticity)  resulting  in  a 
topographically  trapped  cyclone  within  the  canyon  and  vortex  compaction  (negative 
vorticity)  which  results  in  anticyclonic  rotation  in  the  fluid  above  the  canyon.  Finally, 
Klinck  noted  that  as  the  canyon  width  continued  to  be  shortened  (narrowed),  it  ceased 
to  have  much  of  a  perturbing  effect  on  shelf  flow  when  it  reached  approximately  one  half 
of  Rdl.  Limitations,  however,  should  be  kept  in  mind  when  applying  these  model  results 
to  Monterey  Bay.  First,  the  Monterey  Bay  Submarine  Canyon  is  roughly  triangular  in 
shape  (in  plan  view  and  in  its  vertical  cross-section);  second,  the  shelf  region  north  and 
south  of  the  bay  is  relatively  narrow  (approximately  6  km  to  the  north  and  only  2  to  3 
km  to  the  south),  and;  finally,  the  canyon  is  located  in  a  bay  and  does  not  simply  extend 
seaward  across  a  coastal  shelf. 

Typical  values  of  the  first  internal  Rossby  Radius  (Rdl)  were  calculated  for 
the  M  onterey  Bay  using  the  relationship: 
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Figure  15.  Geostrophic  velocity  adjustment  over  a  submarine  canyon:  Side  view 
of  Klinck's  3-levcl  canyon  model  (top).  Bottom  figures  show  cross 
canyon  (x-direction)  outputs  forced  by  a  50  km  wide  current  over  200 


km  (bottom-left)  and  10  km  (bottom-right)  wide  canyons  (Source: 
Klinck,  19S9). 


Substituting  equation  (4)  into  equation  (3),  and  solving  for  the  first  internal  mode  (n=  1) 
yields: 
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Brunt- Vaisala  frequencies  were  obtained  through  2  m  finite  differencing  techniques  and 
averaged  throughout  the  water  column  to  obtain  values  for  A',  which  ranged  from  ap¬ 
proximately  11  X  10"3  Hz  over  the  shelf  area  (very  stable)  to  4  X  10-3  Hz  over  the  can¬ 
yon  axis  (Appendix  D).  Values  for  RdU  obtained  from  equation  (6),  ranged  from  2  km 
over  the  shelf  regions  to  15  km  over  the  axis  of  Monterey  Bay  (Appendix  D).  Canyon 
widths  range  from  approximately  10  to  15  km  at  the  mouth  of  the  canyon,  7  to  9  km 
across  the  center,  and  2  to  3  km  near  the  head  of  the  canyon.  Calculated  values  of  R 
obtained  using  the  relationship: 
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ranged  from  200  km  over  the  shelf  regions  to  1000  km  near  the  center  of  the  bay  (Ap¬ 
pendix  D).  These  values  are  in  close  agreement  with  results  obtained  by  Breaker  and 
Broenkow  (1989)  through  the  averaging  of  3  months  of  hydrographic  data.  Thus  the 
Monterey  Canyon  could  be  defined  as  "narrow"  in  both  the  barotropic  and  baroclinic 
sense,  although  portions  of  the  bay  do  approach  the  width  scale  of  Rdx  . 

C.  PURPOSE  OF  STUDY 

In  light  of  recent  sewage  and  agricultural  run  off  problems  experienced  within  the 
Monterey  Bay,  a  solid  grasp  of  the  mean  circulation  is  essential.  This  study  attempts  to 
quantify  the  mean  flow  field  experienced  within  the  bay  over  a  four  day  period  in  May 
1988  (upwelling  period)  using  conductivity,  temperature,  and  depth  (CTD)  data,  acoustic 
Doppler  current  profiler  (ADCP)  measurements,  current  meter  data  acquired  by  the 
Monterey  Bay  Aquarium,  and  available  NOAA-11  Advanced  Very  High  Resolution 
Radiometer  (AVHRR)  satellite  sea  surface  temperature  imagery.  Phase  relationships 
between  the  tides  and  the  baroclinic  flow  field,  along  with  suggested  internal  wave  fields 
in  the  Monterey  Submarine  Canyon,  are  also  examined. 

Finally,  the  circulation  response  within  the  bay  to  forcing  by  the  observed  mean  flow 
determined  from  the  CTD  and  ADCP  data  will  be  examined  using  an  updated  version 
of  the  baroclinic  model  used  by  Bruner  (1988).  Wind  data,  acquired  both  before  and 
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during  the  period  of  this  study,  will  also  be  used  to  force  this  model,  a  factor  not  previ¬ 
ously  considered. 
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II.  DATA  ACQUISITION  AND  PROCESSING 


A.  COLLECTION 

1.  May  1988  NPS  Student  Cruise 

The  vast  majority  of  the  data  utilized  in  this  study  were  acquired  during  the 
second  half  of  an  NPS  student  cruise  conducted  within  the  Monterey  Bay  during  the 
period  08  to  11  May  1988  onboard  the  Research  Vessel  Point  Sur,  homeported  at  Moss 
Landing,  California.  Two  cross-sections  were  run  perpendicular  to  bottom  topography, 
the  first  running  roughly  north-south  across  the  mouth  of  the  bay  and,  the  second, 
roughly  parallel  but  inshore  of  the  first.  Each  line  contained  CTD  stations  which  were 
sampled  around  the  clock  (Figure  16  on  page  28).  The  seaward  section  intersected  the 
mouth  of  the  Monterey  Submarine  Canyon  with  station  7  situated  directly  over  the  axis 
(hereafter  referred  to  as  "Line  1"),  while  the  landward  section  crossed  the  Monterey 
Canyon  axis  at  station  13  and  the  axis  of  the  Soquel  Canyon  at  station  15  (hereafter 
referred  to  as  "Line  2").  The  first  cross-section  contained  10  CTD  stations,  each  sampled 
a  total  of  five  times,  while  the  shoreward  Line  acquired  data  from  7  hydrographic 
stations,  each  sampled  four  times,  resulting  in  the  overall  collection  of  78  CTD  casts  at 
17  stations  throughout  the  period.  The  cruise  was  designed,  in  part,  to  determine  the 
mean  currents.  Thus,  to  negate  the  effects  of  the  tidal  motions  in  the  water  column, 
each  station  was  occupied  at  approximately  15.5  hour  intervals,  resulting  in  exactly  a 
125%  phase  shift  of  the  semidiurnal  tidal  cycle  each  time  the  station  was  sampled  (Fig¬ 
ure  17  on  page  29).  Each  station  was  therefore  occupied  on  all  phases  of  the  surface 
tide. 

In  addition  to  hydrographic  data,  the  ADCP  was  operated  continuously  for 
each  line  throughout  the  duration  of  the  cruise.  The  RV  Point  Sur  is  outfitted  with  an 
RD  Instruments  DR0150  ADCP  which  utilizes  a  four  beam  JANUS  array  operating  at 
a  frequency  of  150  kHz.  Beams  are  aligned  fore  and  aft,  and  port  and  starboard,  of  the 
RV  Point  Sur.  Raw  data  were  acquired  in  128  four  meter  bins  averaged  over  a  three 
minute  ensemble  interval.  The  ship's  position  was  recorded  from  Loran-C  at  the  end 
of  each  ensemble. 

2.  Other  Supporting  Data 

Cloud  free  NOAA-11  AVHRR  satellite  imagery  was  obtained  from  4,  8,  and  11 
May  1988  and  processed  at  the  Xa\al  Oceanographic  and  Atmospheric  Research  Lab- 
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Figure  16.  CTD  stations  occupied  during  May  1988  NPS  student  cruise:  Stations 
1-10  (designated  "Line  1")  were  each  sampled  5  times  and  stations  10-17 
(designated  "Line  2")  were  each  sampled  4  times.  The  sampling  interval 
was  approximately  15.5  hours. 

oraton  (NOARL)  -  West,  located  in  Monterey,  CA.  GOES  -  West  visual  and  infra-red 
imager}'  were  also  obtained  for  the  period  of  the  cruise.  Finally,  synoptic  weather  charts 
were  collected  from  the  National  Meteorological  Center  (NMC)  for  the  months  of  April 
and  May  1988  via  facsimile  broadcast. 

A  continuous  (5  minute  averaged)  record  of  sea  surface  temperature,  wind  speed 
and  direction,  tidal  height,  and  current  speed  and  direction  was  obtained  from  the 
Monterey  Bay  Aquarium  (MBA)  data  base  for  the  period  21  April  through  20  June  1988 
(Figure  21  on  page  37).  The  MBA  weather  station  is  situated  on  the  roof  of  the 
Aquarium  which  is  located  on  the  shore  of  the  southern  bight  at  36°  37.  F  N,  12F  54.0" 
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Figure  17.  CTD  sampling  interval  overlayed  on  observed  surface  tidal 
heights:  Each  station  is  represented  by  a  separate  symbol.  Stations 
1-10  represent  Line  1,  stations  10-17,  Line  2.  The  tide  data  were  ob¬ 
tained  from  the  MBA  data  base. 

W  at  an  elevation  of  21  m  above  sea  level.  The  MBA  current  meter,  an  electromagnetic 
IntcrOcean  Model  S4,  is  situated  approximately  305  m  seaward  of  MBA  at  a  depth  of 
6.1  in  below  Mean  Lower  Low  Water  (MLLW)  at  36°  37.3'  N,  121°  53.6'  W.  Tide  data 
were  obtained  from  a  tide  gauge  located  on  Monterey  Wharf  #2  (36°  36.3'  N,  121°  53.5' 
W).  The  gauge  is  operated  and  maintained  by  the  National  Oceanic  and  Atmospheric 
Administration  (NOAA).  Raw  tidal  heights  were  corrected  to  MLLW  by  subtracting  a 
2.5  m  offset. 

B.  PROCESSING 

1.  Neil  Broun  Instrument  Systems  (NBIS)  Mark  IIIB  CTD  Probe 

Casts  were  made  at  each  station  utilising  a  NBIS  Mark  II IB  CTD  probe  with 
manufacturer  sensitivities  listed  in  Tabic  2  on  page  30.  Water  samples  were  taken  at  the 
bottom  of  each  cast  and  sea  surface  temperature  measurements  were  obtained  from 
bucket  measurements.  Pressure  offset  (the  pressure  recorded  by  the  CTD  while  sitting 


on  deck)  values  were  recorded  just  prior  to  deploying,  and  upon  recover)’  of,  the  instru¬ 
ment.  Finally,  in  order  to  insure  accurate  salinity  measurements,  the  conductivity  probe 
was  rinsed  with  fresh  water  and  covered  after  each  hydrocast. 


Table  2.  NBIS  MARK  IIIB  CTD  MANUFACTURER  SPECIFICATIONS 


Variable 

Range 

Accuracy 

Resolution 

Pressure 

0  to  3200  dbar 

+  3.2  dbar 

0.05  dbar 

Temperature 

-3  to  +32  0  C 

+  0.005  0  C 

.0005  0  C 

Conductivity 

1  to  65  mmho 

±  0.005  mmho 

0.001  mmho 

Calibrations  were  conducted  prior  to  the  cruise  by  checking  the  CTD 
conductivity,  temperature,  and  pressure  readings  against  standards  in  the  laboratory. 
Differences  obtained  in  this  manner  were  then  averaged  and  fit  to  a  linear  regression 
scheme  in  order  to  obtain  the  coefficients  necessary  to  adjust  the  measurements  made 
by  the  CTD  to  the  reference  standards.  Temperature  coefficients  had  a  slope  of 
0.999363  and  an  intercept  of  0.003435  while  the  pressure  slope  was  calculated  to  be 
1.000000. 

Post-analysis  of  data  in  the  laboratory  consisted  of  first  measuring  the  salinity 
from  the  collected  water  samples  to  obtain  the  "true"  salinity,  then  determining  the  dif¬ 
ference  between  salinities  derived  from  the  CTD  data  and  those  obtained  in  the  labora¬ 
tory  (Figure  18  on  page  31).  Points  more  than  two  standard  deviations  from  the  mean 
(outliers)  were  considered  to  be  the  result  of  random  errors  (human  errors,  ship  motion 
in  areas  of  large  salinity  gradients,  etc.)  and  were  removed.  The  mean  and  standard  de¬ 
viation  were  then  recalculated  and  outliers  greater  than  two  standard  deviations  from  the 
mean  were  once  more  removed.  A  linear  fit  to  the  data  obtained  by  this  method  had  a 
slope  of  1.036714  and  an  intercept  of  -0.347450  which  resulted  in  a  salinity  offset  of 
0.868  psu.  Temperature  and  pressure  values  were  once  again  compared  to  values  ob¬ 
tained  by  laboratory  instruments.  Finally,  a  program  which  used  both  pre-cruise  and 
post-cruise  calibration  coefficients,  including  recorded  pressure  offsets,  was  applied  to 
t'nc  raw'  data  acquired  by  the  CTD  to  obtain  the  values  used  in  this  study.  A  more  de¬ 
tailed  description  of  these  calibration  procedures  is  provided  by  Tisch  (1990). 

2.  Acoustic  Doppler  Current  Profiler  (ADCP) 

After  calibration  of  the  raw  data  and  correcting  for  the  ship's  velocity,  current  velocities 
to  a  depth  of  approximately  400  m  were  calculated.  During  collection,  the  ADCP  data 
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Figure  18.  Scattergram  of  initial  salinity  difference  for  STMAY88:  Plot  repres¬ 
ents  differences  between  C  I'D  and  water  sample  salinity  at  the  bottom 
of  each  cast.  Positive  values  indicate  the  instrument  is  reading  too  low. 
Mean  (/*)  and  standard  deviation  (sd)  are  indicated  (Source:  Tisch, 
1990). 

were  averaged  and  stored  every  3  minutes  along  with  the  ship’s  position.  In  this  way,  a 
unique  position  is  associated  with  each  ensemble  average.  To  recover  cross  sectional 
velocity  profiles,  data  is  again  averaged  at  specified  increments  along  the  entire  line. 
Inaccuracies  inherent  in  this  method  include  errors  induced  through  bad  navigation  data 
and  through  the  averaging  techniques.  For  instance,  if  the  ship  was  stationary  at  one 
of  the  specified  increments,  but  underway  at  the  next  increment,  the  previous  station 
contains  more  data  than  the  latter,  so  averaging  is  somewhat  biased  in  this  respect.  For 
this  reason  the  ADCP  data  acquired  during  this  study  was  averaged  at  4  km  intervals, 
which  approximately  corresponds  to  the  hydrographic  station  spacing,  to  minimize  both 
the  errors  induced  through  navigation,  and  unequal  ensemble  averaging.  A  more  de¬ 
tailed  description  of  converting  raw  ADCP  relative  velocities  into  meaningful  absolute 
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velocities  is  given  by  Kosro  (1985).  Accuracies  in  current  velocities  obtained  using  this 
method  have  been  cited  to  be  on  the  order  of  4  to  5  cm  s*1  (Kosro,  1985). 

Built  into  the  ADCP  system  software  are  constraints  which  determine  whether 
data  is  to  be  considered  "good"  or  not,  namely: 

•  The  shallowest  bin  is  no  good  due  to  side  lobe  interference  reflected  from  the  ship's 
hull  and  the  ocean  surface. 

•  There  must  be  at  least  a  60%  return  in  signal  strength. 

•  Error  velocity,  defined  as  the  difference  in  vertical  velocities  between  orthogonal 
pairs  of  beams,  must  be  within  specifications  in  order  for  data  within  a  certain  bin 
to  be  labeled  "good." 

•  If  the  Automatic  Gain  Control  (AGC),  or  reflected  energy,  increases  greater  than 
10%,  the  system  classifies  this  depth  as  the  bottom  and  negates  data  below  this 
level.  Typically,  if  no  bottom  is  present,  the  AGC  signal  will  simply  decrease  with 
depth  (beam  distance).  In  contrast,  if  a  bottom  is  present  within  range  of  the 
ADCP  beam,  a  jump,  or  increase  in  the  AGC  signal  will  result  and  data  below  this 
depth  will  be  flagged  as  invalid. 

Data  from  the  deepest  good  bin,  determined  by  the  above  constraints,  to  the  ocean 
surface  is  considered  "good."  During  this  cruise,  ADCP  data  were  suppressed  below' 
approximately  75  to  100  m  in  the  areas  near  the  northern  half  of  the  Monterey  Canyon 
due  to  the  AGC  constraint.  During  replay  of  the  acquired  ADCP  data  over  this  area, 
it  was  noted  in  three  out  of  the  five  passes,  that  a  strong  anomalous  "spike"  occurred  in 
the  AGC  plots  at  the  depths  mentioned.  As  expected,  the  AGC  plot  showed  gradual 
decreasing  values  from  the  surface  to  approximately  75  m  where  a  sharp  increase  was 
noted.  Below  this  depth,  values  again  gradually  decreased.  This  spike  wras  assumed  to 
be  the  product  of  biological  activity  of  unknown  origin  which  persisted  near  the  bottom 
of  the  thermocline  in  these  areas  throughout  most  of  the  cruise.  At  depths  below  the 
AGC  "spike,"  data  exhibited  low  propagation  loss  (greater  than  60%  return  signal 
strength)  lending  some  credibility  to  the  fact  that  data  acquired  below  this  layer  may,  in 
fact,  have  been  good,  although  to  what  extent  it  was  biased  or  influenced  by  the  "scat¬ 
tering"  layer  is  unknown.  At  any  rate,  ADCP  data  in  these  areas  were  re-processed  by 
Mr.  Paul  Jessen.  an  NPS  Oceanographer,  with  the  AGC  constraint  suppressed,  to  re¬ 
cover  the  velocities  in  the  area  along  the  northern  wall  of  the  Monterey  Canyon. 


III.  METHODS  AND  RESULTS 


This  chapter  describes  the  results  and  observations  obtained  through  the  analysis 
of  collected  data.  Methods  used  in  obtaining  stated  results  are  provided  in  a  separate 
subsection  unless  they  are  self-explanatory',  such  as  in  the  section  titled  "synoptic 
meteorological  situation."  Evaluation  of  each  data  type  is  discussed  separately  with 
some  limited  comparison  of  results  provided.  An  effort  to  tie  all  results  together  is  given 
in  Chapter  IV.  Results  are  presented  in  the  following  order:  internal  waves  inferred 
through  CTD  data,  mean  geostrophic  velocities,  mean  ADCP  velocities,  description  of 
NOAA  AVHRR  sea  surface  temperature  imagery',  analysis  of  the  synoptic 
meteorological  situation  both  prior  to  and  during  the  May  1988  student  cruise,  and 
finally,  the  methods  and  results  of  applying  the  mean  flow  field  to  a  primitive  equation 
numerical  model  of  the  Monterey  Bay  are  presented.  Some  comparisons  between  the 
MBA  data  and  results  obtained  from  this  study  are  made,  when  applicable. 

A.  INTERNAL  WAVE  FIELD 

1.  Methods 

To  determine  the  extent  of  aliasing  by  internal  wave  processes  and  relate  them 
to  surface  tidal  oscillations,  various  depictions  of  the  fluctuations  of  the  IO  C  isotherm 
were  generated.  The  10'C  isotherm  was  chosen  because  of  its  apparent  proximity  to  the 
center  of  the  thermocline.  Data  were  taken  from  CTD  measurements.  As  a  conse¬ 
quence  of  the  sampling  interval  (approximately  15  h),  only  waves  with  periods  greater 
than  30  hours  (frequencies  smaller  than  the  Nyquist  frequency)  could  be  resolved.  The 
total  structure  (period,  wave  number,  and  phase  speed)  of  the  internal  wave  field, 
therefore,  could  not  be  determined  from  these  data.  First,  a  graph  of  the  10‘C  isotherm 
pressure  across  the  mouth  of  the  bay  (Line  1),  along  with  the  mean  isotherm  pressure 
was  generated  (Figure  19  on  page  34).  To  obtain  a  relationship  between  surface  tidal 
height  and  this  internal  wave  field,  plots  of  the  surface  tidal  height  versus  isotherm 
pressure  were  then  produced  for  each  station  (Figure  20  on  page  35). 

2.  Results 

The  plots  of  the  10'C  isotherm  pressure  changes  during  each  pass  of  Line  1 
(Figure  19  on  page  34)  show  the  active  internal  wave  structure  encountered  during  data 
collection.  The  largest  amplitudes  (up  to  30  m)  were  observed  to  be  in  the  vicinity  of  the 
shelf  breaks  (stations  3  and  9  at  9  and  30  km  on  ihc  graph),  over  the  axis  of  the 
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Figure  19.  Isotherm  variance  along  Line  1:  Top  figure  shows  the  pressure  of  the 
10’C  isotherm  across  the  mouth  of  the  Monterey  Bay  during  each  pass 
(see  legend).  North  is  to  the  left,  south  is  to  the  right.  Bottom  figure 
shows  the  mean  10'C  isotherm  depth  and  the  calculated  isotherm  vari¬ 
ance  at  each  station  along  the  line.  Note  the  large  amplitude  fluctu¬ 
ations  near  the  shelf  break  {variances  of  94.36  m2  and  69.43  m2)  and 
over  the  center  of  the  Monterey  Submarine  Canyon  (variance  of  78.14 
m2). 
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change  in  pressure  of  the  IO'C  isotherm  as  a  function  of  the  observed 
surface  tides,  ’fide  data  were  obtained  from  MBA  data  base. 


Monterey  Submat ine  Canyon  (near  the  24  km  mark  on  the  graph),  and  near  shore  over 
the  shelf  at  station  1.  What  appear  to  be  nodal  points  can  be  seen  near  stations  2  and 
6.  The  plotted  mean  10'C  isotherm  pressute  (Figure  19  on  page  34)  shows  gradual 
deepening  of  the  isotherm  from  the  northern  part  of  the  bay  (station  1)  to  the  center  of 
the  canyon  (station  7)  then  steep  shoaling  from  station  7  to  the  southern  shelf  at  station 
10.  Calculated  \ariance  shows  high  values  in  the  aieas  of  large  amplitude  fluctuations 
picviously  mentioned.  In  general,  the  observed  internal  wave  amplitudes  were  smaller 
in  Line  2  (stations  11-17)  with  most  amplitudes  falling  within  the  10  to  15  m  range.  This 
was  somewhat  surptising  considering  past  studies  (Shea  and  Broenkow,  19S2,  and 
Biocnkow  and  McKain,  lc)72)  which  showed  higher  amplitude  internal  waves  nearer  the 
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head  of  the  canyon  than  towards  the  mouth.  This  again  may  simply  be  the  result  of  the 
relatively  large  station  sampling  interval,  i.e.  the  internal  waves  in  Line  1  may  have  had 
a  different  frequency  or  phase  relation  with  the  surface  tidal  heights  than  those  en¬ 
countered  in  Line  1,  resulting  in  poorer  resolution  of  these  waves  in  Line  2  data. 

The  10’C  isotherm  pressure  versus  tidal  height  at  each  station  of  Line  1 
(Figure  20  on  page  35)  show  no  direct  surface  tidal  correlation.  Large  differences  in 
pressure  can  be  seen  even  though  the  surface  tidal  height  is  at  nearly  the  same  level. 
For  instance,  at  station  8,  when  the  tide  is  approximately  0.5  m,  the  pressure  of  the 
10'C  isotherm  ranges  from  20  to  45  dbar  (note  1  dbar^I  m  in  these  ranges).  Similar 
cases  of  varying  isotherm  pressure  at  constant  tidal  heights  can  be  seen  at  other  stations 
(station  1  had  isotherm  pressures  which  varied  by  as  much  as  25  dbar  during  two  sepa¬ 
rate  periods  when  the  tide  was  at  the  same  height).  If  internal  tides  were  in  phase  with 
the  external  tides,  these  graphs  should  basically  show  a  positive  sloping  line  from  the 
bottom  left  to  the  upper  right  (i.e.  the  isotherm  would  rise  and  fall  with  the  external  tidal 
height).  Station  9,  near  the  southern  shelf  break,  is  the  only  station  which  comes  close 
to  this  description,  where  all  points  nearly  lie  on  a  straight  line.  A  cross-spectral  analysis 
between  one  week  of  continuous  thermistor  data  near  the  canyon  head  and  surface  tidal 
heights,  conducted  by  Shea  and  Broenkow  (1982),  revealed  a  7  h  phase  lag.  A  similar 
analysis  of  approximately  3  months  of  thermistor  data,  from  mid  April  to  mid  June  1988, 
collected  from  the  MBA  sensor  at  a  depth  of  6  m  in  the  southern  bight  (Figure  21  on 
page  37),  revealed  an  8  h  phase  lag  at  the  predominant  semi-diurnal  12.8  h  period. 
However,  phase  relationships  between  surface  tidal  heights  and  internal  waves  near  the 
canyon  head  and  over  the  shelf  regions  more  than  likely  differ  from  those  in  the  deeper 
regions  of  the  bay,  especially  over  the  canyon. 

In  summary,  the  internal  wave  field  encountered  during  the  08  to  11  May  1988 
data  collection  period  was  highly  energetic,  especially  along  Line  1.  Amplitudes  varied 
across  the  bay,  but  appeared  to  reach  a  maximum  over  the  shelf  breaks  which  lends 
support  to  Shea  and  Brocnkow's  (1982)  proposed  model  of  processes  occurring  near  the 
edge  of  the  Monterey  Canyon  (see  Figure  10  on  page  17).  Although  the  direction  of 
propagation  of  these  internal  waves  could  not  be  determined  from  the  available  mass 
field  data,  the  high  variability  across  the  mouth  of  the  bay  would  seem  to  suggest  at  least 
some  cross  canyon  component  exists.  Possibly  the  most  useful  information  gained  from 
the  above  analysis  is  the  fact  that  the  isotherm  slopes  created  by  the  internal  waves  ex¬ 
ceed  the  mean  slopes  due  to  geostrophic  adjustment  by  two  orders  of  magnitude.  This 
implies  that  geostrophic  velocities  derived  from  one  pass  alone  is  virtually  useless  be- 
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Figure  21.  Acquired  MBA  ocean  and  atmospheric  data  21  April  -  20  June 
1988:  MBA  data  is  averaged  and  stored  ever}’  5  minutes.  Above 
curves  were  generated  using  a  3  hour  running  mean  filter.  Dashed  lines 
represent  total  mean. 
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cause  the  internal  wave  "noise"  will  quite  effectively  mask  the  mean  signal.  Therefore, 
some  sort  of  averaging  scheme  is  a  necessity.  The  next  sections  describe  the  geostrophic 
flow  obtained  from  the  mean  mass  field  and  the  mean  ADCP  field. 

B.  MEAN  TEMPERATURE,  SALINITY,  DENSITY,  AND  GEOSTROPHIC  FLOW 
1.  Method 

Analysis  of  the  hydrographic  data  acquired  during  the  May  1988  NPS  student 
cruise  was  conducted  in  two  steps.  First,  data  from  each  pass  were  run  through  a  com¬ 
puter  program  which  utilizes  horizontal  density  gradients  between  stations  via  the  ther¬ 
mal  wind  equation: 


dv  8  cp 
dz  jPq  cx 

to  calculate  geostrophic  velocities  throughout  the  water  column.  The  second  part  con¬ 
sisted  of  averaging  the  mass  field.  This  was  accomplished  by  averaging  the  temperatures 
and  salinities  at  each  station  (from  5  samples  in  Line  1  and  4  samples  in  Line  2).  From 
this  mean  field,  geostrophic  velocities  and  cross-sections  were  then  produced  (using 
equation  (8))  for  both  Lines  (Figure  24  on  page  42,  Figure  26  on  page  47,  and  Appen¬ 
dix  A).  Because  of  the  sharp  topography,  the  assumed  level  of  no  motion  used  in  all 
geostrophic  velocity  calculations  was  the  deepest  common  depth  between  stations. 

As  previously  stated,  the  sampling  interval  was  designed  so  that  tidal  oscillations 
could  be  averaged  out  of  the  water  column.  To  measure  the  success  of  this  particular 
aspect  of  the  cruise,  the  mean  tidal  height  for  the  period  was  calculated  using  MBA  data, 
and  the  variance  of  tidal  heights  (from  observed  heights  at  each  station  during  sequential 
passes)  about  this  mean  were  then  derived  (Figure  22  on  page  39).  The  calculated  var¬ 
iance  was  relatively  small  at  3.45  x  10'3m2  for  Line  1,  and  1.12  x  10_3m2  for  Line  2 
(Figure  22  on  page  39).  Initially,  averaging  schemes  using  various  pass  combinations 
were  tested  in  an  attempt  to  minimize  the  variance  of  station  tidal  height  about  the 
MBA  tidal  mean.  Although  averaging  all  line  passes  did  not  produce  the  minimum 
variance  about  the  MBA  mean,  it  was  decided  to  use  all  five  because  it  was  apparent 
that  the  internal  tides  were  not  in  phase  with  the  surface  tidal  heights  (discussed  in  sec¬ 
tion  on  internal  waves). 

To  further  measure  the  "success"  of  this  averaging  process,  the  net  mean  mass 
transport  across  each  line  was  calculated.  In  the  mean,  the  net  geostrophic  mass  trans¬ 
port  into  and  out  of  the  bay  should  be  zero,  assuming  that  the  barotropic  transport  was 
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Figure  22.  Effective  tides:  'I  lie  plot  shows  surface  tidal  fluctuations  from  5  min¬ 
ute  averages  of  observed  tidal  height  obtained  from  the  MBA  data 
base.  Solid  line  shews  mean  MBA  tidal  height,  dotted  line  shows  cf- 
fccti.e  surface  tide  derived  thiough  averaging  of  tidal  heights  observed 
at  each  hydrographic  station  (see  also  Figure  17  on  page  29).  Variance 
wul.  respect  to  the  MBA  mean  shown  at  the  bottom  right. 


minimal,  Calculated  net  volume  tiansport  was  .152  Sv  for  Line  1  and  .012  Sv  for  Line 
2,  both  considered  balanced  within  the  resolution  of  the  methods  used  within  this  study 
for  reasons  discussed  shortly. 

Finally,  the  standard  error  of  the  mean  gcostrophic  velocity  (aj)  was  calculated 
(Figure  24  on  page  42  and  Figure  26  on  page  47)  using  the  relationship: 


°x  ~ 


(9) 


where  the  standard  deviation  (a)  is: 
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The  derived  geostrophic  velocities  (X)  between  each  station  on  each  individual  pass  was 
used  as  the  population  (n). 

In  addition  to  geostrophic  cross-sections,  maps  of  the  temperature  fields  at  0, 
50,  and  200  dbar,  dynamic  heights  at  0  dbar  relative  to  50  dbar  and  100  dbar  and  200 
dbar  relative  to  500  dbar,  as  well  as  cross  sections  of  temperature,  salinity,  and  density, 
were  generated  (Figure  27  on  page  48,  Figure  28  on  page  49,  Figure  23  on  page  41,  and 
Figure  25  on  page  46),  The  dynamic  height  topographies  were  then  used  to  infer  the 
geostrophic  flow  patterns. 

In  summary,  three  parameters  were  calculated  for  each  data  set  to  determine  the 
effectiveness  of  the  averaging  process  in  extracting  the  mean  flow:,  station  tidal  variance 
about  the  MBA  mean  tidal  height,  section  mean  volume  transport,  and  the  standard 
error.  Experimentation  with  the  averaging  of  various  pass  combinations,  to  minimize 
the  tidal  variance  of  the  data  set  about  the  MBA  mean  tidal  height,  did  not  reduce  the 
net  volume  transport  (i.e.  increase  the  implied  accuracy  of  the  mean)  through  the  sec¬ 
tion,  This  further  illustrates  the  fact  that  the  internal  wave  field  was  not  in  phase  with 
the  surface  tidal  heights.  If  they  were  in  phase,  the  most  accurate  mean  flow  (smallest 
net  volume  transport)  would  naturally  correlate  with  the  averaging  scheme  which  yielded 
the  smallest  variance  about  the  mean  MBA  tidal  height  during  the  period.  Since  this 
was  not  the  case,  the  decision  to  average  all  5  passes  in  Line  1  and  all  4  passes  in  Line 
2  was  arrived  at,  from  a  purely  statistical  viewpoint,  to  maximize  the  number  of  degrees 
of  freedom.  Statistically,  this  method  worked  well  since  calculated  standard  errors  were 
smaller  than  the  derived  mean.  With  this  in  mind,  the  overall  data  acquisition  method 
was  successful  at  minimizing  the  variance  of  the  tidal  heights  and  creating  a  nearly  bal¬ 
anced  mean  volume  transport  by  virtue  of  pass  averaging  alone.  The  method  was  not 
entirely  successful  however,  at  systematically  (deterministically)  reducing  the  high  fre¬ 
quency  internal  motions. 

2.  Results 

a.  Mean  Property  Cross  Sections,  Geostrophic  Flow,  and  Mass  Transport 

(1)  Line  1.  The  Line  1  mean  temperature  and  salinity  cross  sections 
(Figure  23  on  page  41)  indicate  an  area  of  warm,  fresh  water  at  the  surface  (greater  than 
12'C)  extending  across  the  mouth  of  the  bay  from  station  5  to  station  9.  This  warm, 
fresh  "pool"  of  water  appears  well  imbeded  in  the  vertical  structure  and  extends  down 
to  approximately  100  dbar,  as  the  8.5  and  9'C  isotherms  indicate.  At  deeper  depths,  over 
the  Monterey  Submarine  Canyon,  upwelling  of  cold  water,  most  notably  along  the 
northern  shelf  break,  appears  to  dominate  the  temperature  field  as  shown  by  the  8  and 
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Figure  23.  Line  1  mean  temperature,  salinity,  and  density  cross  sections:  Figures 
depict  mean  temperature  (°C,  top),  salinity  (psu,  middle),  and  density  ( 
kg  nr1,  bottom)  sections  for  Line  1.  Caution  should  be  noted  in  inter¬ 
preting  deep  isolines  adjacent  to  the  canyon  walls.  These  steep  slopes 
more  than  likely  do  not  exist  and  are  merely  the  product  of  the  con¬ 
touring  program  used  to  develop  these  graphs.  Station  positions  are 
indicated  across  the  top  of  the  salinity  section  (sec  also  Figure  16  on 
page  28  for  locations). 
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Figure  24.  Line  1  mean  grostrophic  flow  field  and  standard  error:  Top  figure  is  the 
mean  geostrophic  flow,  bottom  is  the  standard  error..  Hydrographic 
stations  are  notated  at  the  top  of  the  section  and  are  situated  across  the 
mouth  of  Monterey  Bay  (Line  1,  see  also  Figure  16  on  page  28  for 
station  locations).  Solid  lines  depict  flow  into  the  bay,  dashed  lines 
show  flow  out  of  the  bay.  Geostrophic  velocities  were  calculated  using 
the  deepest  common  depth  between  stations  as  the  level  of  no  motion. 
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8.5’C  isotherms.  The  salinity  cross  section  shows  the  fresher  water  (33.65  psu)  located 
between  stations  6  and  9.  Again,  some  indication  of  upwelling  of  denser,  saltier  water 
from  the  canyon  onto  the  shelf  region  is  indicated  (again  most  notably  over  the  northern 
shelf)  by  the  34.05  and  34.00  isohalines.  Density  anomaly  cross  sections  (Figure  23) 
reflect  the  temperature  and  salinity  fields  by  indicating  relatively  light  water  located  in 
the  top  100  dbar  of  the  water  column  and  centered  nearly  over  the  canyon  axis  with 
denser  water  extending  up  onto  the  northern  shelf  along  the  canyon  wall  from  a  depth 
of  200  m. 

The  calculated  mean  geostrophic  flow  (Figure  24  on  page  42,  and 
Appendix  A)  reveals  two  apparently  separate  flow  patterns.  At  the  surface  (0  to  ap¬ 
proximately  75  dbar)  flow  is  split  into  two  regimes:  broad  inflow  across  the  northern 
half  of  the  Monterey  Bay  from  stations  1  to  7;  and,  outflow  across  the  southern  half  of 
the  bay  from  stations  7  to  10.  Inflow  extends  to  a  somewhat  greater  pressure  (125  dbar) 
than  the  outflow  region  (50  dbar).  Maximum  surface  inflow  was  calculated  at  15.6  cm 
5_I  between  stations  3  and  4  and  a  maximum  outflow  of  -18.6  cm  s-’,  situated  between 
stations  9  and  10,  was  indicated.  Corresponding  standard  errors  of  the  mean  at  these 
locations  were  9.9  and  4.8  cm  .s-1  respectively  (Figure  24  on  page  42).  Analysis  of  the 
current  speed  and  direction  from  the  MBA  current  meter  in  the  southern  bight  area 
(closest  to  station  10)  over  the  period  showed  northwest/southeast  flow  reversals  with 
the  seaward  flow  predominating  (Figure  21  on  page  37).  Also  evident  was  the  corre¬ 
lation  between  the  strongest  flow  (highest  velocities)  and  the  seaward  direction.  Mean 
current  speed  and  direction,  derived  from  these  data  during  the  period  of  the  May  cruise, 
was  4.5  cm  s~‘  towards  245T,  or  seaward.  This  outflow  is  smaller  than  the  noted 
geostrophic  flow  between  stations  9  and  10,  not  unusual  considering  the  current  meter's 
close  proximity  to  shore.  The  second  flow  pattern  is  located  over  the  Monterey  Sub¬ 
marine  Canyon  between  100  and  600  dbar  and  consists  of  outflow  along  the  northern 
wall  of  the  canyon,  with  inflow  indicated  over  the  southern  half  of  the  canyon.  The 
cores  are  found  at  250  dbar  (downcanyon  flow  of  10.6  cm  .r1)  between  stations  5  and 
6,  and  at  150  dbar  (upcanyon  flow  of  10.1  cm  s*1)  under  station  8.  Standard  errors  cal¬ 
culated  at  these  positions  and  depths  yielded  values  of  8.3  and  9.1  cm  5'1  respectively, 
showing  that  the  deep  flows  were  barely  resolved  with  this  averaging  technique.  The 
mean  flow  pattern  shows  strong  vertical  shears  near  stations  6  and  8,  along  with  weaker 
horizontal  shear  zones  over  the  center  of  the  bay.  Standard  errors  at  each  level 
(Figure  24  on  page  42)  indicate  that  the  derived  surface  mean  flow  is  much  more  re- 


solved  (especially  the  outflow  in  the  southern  half  of  the  bay,  past  the  Monterey  Penin¬ 
sula)  than  the  mean  flow  at  depth,  which  barely  surpasses  the  standard  error  values. 

The  net  mass  transport  calculated  for  the  mean  flow  was  .152  Sv. 
This  result  was  not  surprising  and  appears  to  be  within  the  resolution  of  the  geostrophic 
methods  exploited  for  a  number  of  reasons,  namely: 

•  The  geostrophic  method  does  not  take  into  consideration  any  non-linear  effects, 
which  can  produce  errors  -up  to  10%  in  some  areas  of  the  ocean  such  as  near 
eddies,  high  shear  zones,  etc.  Because  of  the  Monterey  Bay's  geometry,  and  ap¬ 
parently  active  internal  wave  system,  non-linear  processes  more  than  likely  are  at 
work  and  are  neglected  in  geostrophv.  Rough  estimates  of  the  extent  of  non¬ 
linearity  were  obtained  through  calculations  of  the  Rossby  number  (R„): 

Ro=jr>  on 

where  /,  the  coriolis  parameter  is  0(  10~^);  velocities  (V)  are  0(.l  m  5"‘);  and  typical 
length  scales  (L),  estimated  by  the  computed  first  internal  Rossby  radius  of  defor¬ 
mation  (Appendix  D),  of  0(5  to  15  km).  Values  range  from  approximately  .1  to 
.2,  implying  that  up  to  20%  of  flow  may  be  affected  by  non-linear  processes  in 
Monterey  Bay. 

•  Recent  studies  (Tisch,  1990)  indicate  that  geostrophic  velocities  calculated  between 
stations  with  close  spacing  typically  have  higher  errors  than  those  calculated  at 
larger  grid  spacing.  When  station  spacing  is  less  than  approximately  4  to  5  km, 
which  was  the  case  in  this  study  for  nearly  all  the  stations,  small  errors  in  posi¬ 
tioning  and  measurement,  along  with  the  ship  drift,  tend  to  propagate  larger  errors 
in  geostrophic  velocity  calculations  than  values  derived  from  stations  situated  far¬ 
ther  apart. 

•  Geostrophic  velocities  were  calculated  using  the  common  deepest  depth  between 
stations.  Thus,  no  deep  flow  (greater  than  1000  dbar)  is  resolved  in  this  study. 

•  Approximately  4  km  of  the  bay,  on  b  ‘h  ends  of  Line  1,  were  not  covered.  This 
was  particularly  important  in  the  southern  part  of  the  bay  where  strong  outflow 
was  indicated  right  up  to  station  10.  If  strong  outflow  continued  right  up  to  the 
Monterey  Peninsula,  a  mass  transport  rate  of  roughly  -.036  Sv  ( 
4000  m  x  50  m  x  0.18  m.r1)  was  missed  in  Line  1. 

•  The  barotropic  flow  was  not  resolved  through  geostrophv.  Rough  estimates  of  this 
component  were  obtained  by  first  calculating  the  surface  area  within  the  bay  and 
bounded  by  Line  1  (882  kmJ),  then  computing  the  flow  that  would  be  required  to 
raise  the  water  level  1.5  m  (approximate  tidal  range  of  the  bay)  and  the  ensuing 
transport  rate.  From  these  calculations,  barotropic  flow  was  determined  to  be  on 
the  order  of  .1  to  .3  cm  y'  winch  would  account  for  a  mass  transport  of  ±  ,06  Sv. 

Therefore,  errors  on  the  order  of  .1  Sv  can  be  accounted  for  by  reasoning  of  missed 
coverage  across  the  bay,  deep  flow,  and  barotropic  transport.  Errors  of  O(.l)  Sv  are  also 
probably  a  realistic  estimate  of  "rrors  propagated  through  geostrophv  (i.e.  non-linearity, 
station  -pacing,  assumed  level  of  no  motion,  etc.;.  In  view  of  this,  and  calculations  of 
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the  standard  error,  the  derived  mean  flow  field  appears  balanced  within  the  limitations 
of  available  data  and  methods.  Within  the  confines  of  this  study,  the  term  "mean  flow" 
refers  to  the  flow  during  the  period  from  08-11  May  1988  and  may  not  represent  a 
"typical"  flow  pattern  in  the  Monterey  Bay  (if  such  a  pattern  even  exists).  An  attempt 
will  be  made  later  to  relate  this  "mean  flow"  to  processes  believed  to  have  been  occurring 
at  the  time  these  particular  data  were  collected. 

(2)  Line  2.  The  mean  Line  2  temperature  and  salinity  cross  sections 
(Figure  25  on  page  46)  show  smaller  perturbations  in  the  thermocline  than  the  Line  1 
data.  Somewhat  warmer  water  can  be  seen  at  the  surface  over  the  northern  and  south¬ 
ern  shelf  regions  but  overall  isotherms  are  relatively  level.  This  is  in  contrast  to  the  Line 
1  isotherms  which  show  a  temperature  maximum  over  the  center  of  the  canyon.  This 
may  be  a  consequence  of  increased  heating  of  the  shallower  shelf  waters  which  rapidly 
dominate  the  bathymetry  as  the  canyon  head  is  approached.  This  warmer  water  situated 
over  the  shelf  regions  of  the  inner  portions  of  the  bay  can  also  be  seen  on  the  OS  May 
1988  satellite  imagery  (Figure  32  on  page  58),  discussed  later.  The  salinity  field  also 
displays  relatively  level  isohaline  contours  within  the  thermocline  except  for  somewhat 
fresher  water  indicated  over  the  southern  shelf  region  (Figure  25  on  page  46).  The  cross 
section  of  mean  density  shows  somewhat  lighter  (less  dense)  water  situated  over  the 
northern  and  southern  shelf  regions  as  indicated  by  the  25.6  isopycnal. 

The  mean  currents  derived  for  Line  2  shows  broad  inflow  over  the 
center  of  the  section  between  stations  15  and  12  with  a  maximum  at  the  surface  between 
stations  12  and  13  of  7.4  cm  s~'  and  a  subsurface  maximum  at  50  dbar  of  6.2  cm  s~' 
(Figure  26  on  page  47).  Standard  errors  were  calculated  at  2.7  and  1.4  cm  s-1  respec¬ 
tively..  The  calculated  mean  outflow  was  again  narrower,  but  stronger,  than  the  inflow, 
and  confined  mainly  to  the  southern  portion  of  the  bay,  over  the  shelf  region,  between 
stations  11  and  10.  Maximum  outflow  was  calculated  at  15.9  cm  s'1  in  this  area  with  a 
standard  error  of  6.2  cm  .r1  As  noted  earlier,  both  mean  velocities  (hence,  mass  trans¬ 
port)  and  velocities  calculated  during  each  pass,  were  smaller  than  those  calculated  for 
Line  1,  implying  that  much  of  the  circulation  derived  in  Line  1  is  closed,  and  does  not 
extend  significantly  into  the  heart  of  the  bay  itself.  The  net  mass  transport,  derived  from 
mean  flow  values,  was  .012  Sv  (an  order  of  magnitude  smaller  than  that  of  Line  1). 
Again,  similar  arguments  for  not  attaining  mass  balance  closure  in  Line  1  can  also  be 
applied  to  Line  2. 

(1/  Mean  Temperature  and  Dynamic  Height  Maps.  The  mean  surface 
temperature  map  (Figure  27  on  page  48),  interpreted  by  hand,  depicts  two  pools  of 
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Figure  25.  Line  2  mean  temperature,  salinity,  and  density  cross  sections:  Figures 
depict  mean  temperature  (°C,  top),  salinity  (psu,  middle),  and  density 
(kg  nr3,  bottom)  sections  extending  roughly  north-south  across  the 
center  of  the  Monterey  Bay  (Line  2)  and  shows  both  the  Soquel  (left) 
and  Monterey  (right)  Submarine  Canyons.  Caution  should  be  noted 
in  interpreting  aecp  isolines  adjacent  to  the  canyon  walls.  These  steep 
slopes  more  than  likely  do  not  exist  and  are  merely  the  product  of  the 
contouring  program  used  to  develop  these  graphs.  Station  positions 
arc  indicated  acioss  the  top  of  the  temperature  section  (see  also 
Figure  16  on  page  28  for  locations). 
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Figure  26.  Line  2  mean  geostrophic  flow  field  and  standard  error:  Top  figure  is 
mean  geostrophic  flow,  bottom  is  standard  error.  Hydrographic 
stations  are  notated  at  the  top  of  the  section  and  run  roughly  north- 
south  across  the  center  of  Monterey  Bay  (Line  2,  see  Figure  16  on  page 
28  for  station  location).  Solid  lines  depict  flow  into  the  bay,  dashed 
lines  show  flow  out  of  the  bay.  Geostrophic  velocities  were  calculated 
using  the  common  deepest  depth  between  stations  as  the  level  of  no 
motion. 
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Figure  27.  Mean  temperature  fields:  Mean  temperature  fields  ('Q  on  the  0  dbar 
(top-left),  50  dbar  (top-right),  and  200  dbar  (bottom)  pressure  surfaces. 


Figure  28.  Mean  dynamic  height  fields:  Mean  dynamic  heights  at  0  dbar  relative 
to  50  dbar  (top-left)  and  to  100  dbar  (top-right)  and  at  200  dbar  relative 
to  500  dbar  ('bottom'). 
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warm  water  at  the  surface.  The  first  is  centered  over  the  mouth  of  the  canyon  with  a 
12.45'C  maximum  at  station  7,  and  the  second  appears  to  be  contained  in  the  northern 
bight  area  with  a  12.43'C  maximum  at  station  17.  The  pool  of  warm  water  centered  at 
station  7  maintains  its  identity  down  to  50  dbar  (the  approximate  base  of  the  mixed 
layer),  with  station  7  retaining  the  maximum  50  dbar  temperature  at  9.61’C.  Evidence 
of  the  warmer  water  present  at  the  surface  in  the  vicinity  of  the  northern  bight  disap¬ 
pears  at  depth.  This  may  be  a  consequence  of  the  previously  mentioned  surface  heating 
of  shelf  waters  in  the  northern  bight,  which  typically  does  not  extend  down  through  the 
thermocline.  Other  satellite  SST  data  (Tracy,  1990)  also  show  enhanced  surface  heating 
in  the  northern  bight,  possibly  since  the  area  is  sheltered  from  the  northwesterly  winds 
by  the  Santa  Cruz  mountains.  Temperatures  at  200  dbar  are  slightly  colder  over  the 
canyon,  especially  towards  the  northern  wall,  as  evidenced  by  the  minimum  temperature 
(7.65 ’C)  at  station  6. 

The  mean  dynamic  topographies  (Figure  28  on  page  49)  at  the  sur¬ 
face  relative  to  50  and  100  dbar,  show  higher  relative  heights  near  the  mouth  of  the  bay, 
roughly  centered  over  the  canyon  axis  with  maximum  values  of  .192  and  .194  dy  m  at 
stations  6  and  7  respectfully  and  implies  anticvclonic  geostrophic  circulation  within  the 
mixed  layer  over  the  Monterey  Bay.  Deeper  circulation  is  cyclonic  with  low  dynamic 
heights  over  the  canyon  axis  at  stations  6  and  7  of  .388  and  .386  dy  m  respectively, 

C.  MEAN  FLOW  FROM  ADCP  DATA 
1.  Method 

As  an  alternative  method  in  determining  a  mean  flow  (and  to  compare  with 
gcostrophically  derived  mean  velocities),  ADCP  data  were  analyzed  for  Line  1.  This  was 
accomplished  by  first  breaking  the  ADCP  data  up  into  two  portions,  the  first  with  the 
ship  on  a  heading  of  180’  T  (i.e.  between  stations  1  and  7),  and  the  second  part,  when 
the  ship  was  on  a  heading  of  roughly  145T  (between  stations  7  and  10).  This  was  nec¬ 
essary  so  the  applied  program  could  compute  velocities  that  would  correspond  to  the 
directions  obtained  from  geostrophic  calculations.  Data  were  averaged  every  4  km  along 
each  portion  of  Line  1  (stations  1-7  and  7-10).  After  the  velocities  for  each  pass  were 
computed,  they  were  again  averaged,  over  all  5  passes.  Horizontal  distance  averaging 
was  done  since  the  3  minute  profiles  were  not  co-!ocatcd  from  pass  to  pass  (see 
Figure  29  on  page  51).  Finally,  mean  files  from  both  portions  (stations  1-7  and  7-10) 
were  combined  to  produce  one  plot  of  ADCP  mean  velocities  that  would  extend  from 
station  1  to  10  so  that  direct  compaiisons  could  be  made  with  the  geostrophic  velocities. 
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Figure  29.  Line  1  ADCP  averaging  sciieme:  Top  figure  shows  averaging  interval 
for  stations  1-7  and  the  bottom  figure  for  stations  7-10.  Horizontal 
averaging  as  well  as  vertical  averaging  of  each  pass  was  done  to  mini¬ 
mize  errors  induced  through  slight  offsets  of  data  profiles  along  each 
pass. 
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(see  Figure  30  on  page  53,  and  Appendix  B  for  numerical  outputs).  Standard  errors 
were  calculated  using  the  mean  ADCP  velocities  and  the  velocities  obtained  during  each 
separate  pass  in  the  same  manner  as  was  done  with  the  geostrophic  data. 

2.  Results 

a.  ADCP  Velocities  Obtained  During  Successive  Passes  of  Line  1 

Contrary  to  velocities  obtained  through  geostrophy  from  individual  passes, 
ADCP  velocities  will  not  be  influenced  by  mass  fluctuations  and  detect  the  total  veloci¬ 
ties  (both  the  geostrophic  and  ageostrophic  components)  present  at  the  time.  Although 
ADCP  velocities  collected  during  each  pass  undoubtably  are  also  influenced  by  internal 
wave  motions,  it  is  felt  that  results  obtained  during  each  pass  warrant  a  brief  description. 

ADCP  velocities  obtained  during  successive  passes  of  Line  1  ranged  from 
light  flow  up  to  20  cm  s_1  in  pass  3.  Numerous  current  reversals  between  successive 
passes  were  noted  and  effectively  demonstrate  the  active  dynamical  processes  at  work 
within  the  bay.  Also,  large  shear  zones,  both  in  the  vertical  and  horizontal,  were  in  ev¬ 
idence  in  virtually  all  passes.  In  general,  ADCP  velocities  were  less  than  the  geostrophic 
velocities  calculated  during  each  pass,  especially  at  depth.  This  is  not  surprising  since 
the  calculated  geostrophic  response  to  the  large  swings  in  isopycnal  slopes,  created  by 
the  large  amplitude  internal  waves,  would  tend  to  produce  unrealistically  large  velocities, 
especially  in  the  thermocline. 

No  definite  relationship  between  the  surface  tidal  heights  and  flow  within 
the  bay  could  be  determined  from  the  ADCP  passes.  During  pass  3  of  Line  1,  the  sur¬ 
face  tide  was  in  the  ebbing  phase  between  stations  4  and  9  (see  Figure  17  on  page  29), 
going  from  Higher  High  Water  (HHW)  at  station  4  to  Higher  Low  Water  (HLW)  at 
station  9.  During  this  period  the  ADCP  data  depicted  virtually  all  flow  going  upcanyon 
(into  the  bay)  with  a  20  cm  sr'  maximum  centered  over  the  canyon  axis  at  a  pressure  of 
approximately  100  dbar.  On  the  next  pass  (pass  4),  the  tide  was  in  a  flooding  stage  be¬ 
tween  station  5  (at  Lower  Low  Water  (LLW))  and  station  10  (at  Higher  High  Water 
(HHW)).  In  this  case  a  total  reversal  in  flow'  was  noted  with  no  evidence  of  the  20  cm 
s~‘  maximum  previously  noted,  Virtually  all  flow'  was  directed  out  of  the  bay  with  a  -18 
cm  s''  maximum  at  the  surface  near  stations  8  and  9.  These  two  passes  tend  to  support 
Gatje  and  Pizingcr's  (1965)  current  meter  observations  which  correlated  downcanyon 
flow  with  surface  flood  tides,  and  upcanyon  flow  with  surface  tidal  heights  in  the  ebbing 
stage.  The  apparent  correlation  between  surface  tidal  heights  and  bay  flooding  flushing 
periods  implied  by  these  two  passes  docs  not  hold  in  the  other  passes  however.  For  in¬ 
stance,  pass  1,  again  in  a  flooding  stage  from  LLW  at  station  1  to  Lower  High  Water 
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Figure  30.  Line  1  mean  ADCF  velocities  and  standard  error:  Top  figure  shows 
mean  ADCP  velocities  along  Line  1.  Negative  values  indicate  flow  out 
of  the  bay,  positive  values  indicate  flow  into  the  bay.  Bottom  figure 
represents  standard  error  values  calculated  for  the  ADCP  velocities. 
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(LHW)  at  station  10,  depicted  outflow  over  the  northern  shelf  and  along  the  north  can¬ 
yon  wall,  but  inflow  over  the  :enter  of  the  canyon  and  over  the  southern  shelf.  Also, 
during  pass  2,  which  went  from  HHW  at  station  4  to  nearly  LLW  at  station  10  (ebbing 
tides),  the  ADCP  velocity  section  depicted  outflow  at  the  surface  but  inflow  at  100  dbar. 
These  passes  tend  to  support  observations  made  by  Shepard,  et  al.  (1979),  and  by 
Broenkow  and  Mckain  (1972),  that  the  baroclinic  and  barotropic  tides  are  generally  not 
in  phase. 

b.  Mean  ADCP  Velocities  and  Standard  Errors 

The  mean  ADCP  velocity  structure  shows  broad  surface  inflow  over  the 
northern  half  of  the  bay  from  station  1  to  station  7  with  the  core  situated  approximately 
17  km  south  of  station  1  having  a  maximum  velocity  of  15  cm  s~'.  A  narrower  region 
of  outflow  is  revealed  along  the  southern  half  of  the  bay,  south  of  station  7,  with  a 
maximum  velocity  of  10.5  cm  s~'  centered  approximately  34.5  km  south  of  station  1 
along  the  Line  1  track.  Weak  outflow  (-2  cm  s'1)  is  indicated  near  the  northern  canyon 
wall  at  182  dbar  with  stronger  inflow  (13.7  cm  s~'  at  150  dbar)  indicated  over  the 
southern  half  of  the  canyon  axis  between  stations  7  and  8.  Standard  errors  ranged  from 
1.7  cm  s_l  for  the  surface  inflow  center  to  4.1  cm  s-‘  for  the  core  of  the  deeper  inflow 
between  stations  7  and  8  (see  Figure  30  on  page  53).  This  represents  a  better  signal  to 
noise  ratio  for  the  ADCP  data  than  for  the  velocities  derived  through  geostrophy  and 
implies  that  ADCP  techniques  may  have  an  advantage  over  geostrophy  when  attempting 
to  determine  mean  flow  patterns  within  the  Monterey  Bay. 

c.  Comparison  of  Mean  ADCP  and  Geostrophic  Velocities 

Table  3  on  page  55  summarizes  the  maximum  velocity  values  and  their  lo¬ 
cations,  obtained  by  using  both  mass  averaged  geostrophy  and  by  averaging  the  ADCP 
velocities.  Considering  the  major  differences  in  the  two  techniques,  the  results  of  both 
methods  compare  favorably  in  the  ..  ation  and  in  the  magnitude  of  the  flow  centers. 
Horizontal  differences  in  the  location  of  maximum  velocities  were  at  most  4  km's  apart. 
This  is  expected  since  the  ADCP  continuously  collects  and  averages  velocity  data,  while 
geostrophic  velocities  can  only  be  calculated  for  the  water  column  centered  between 
CTD  stations.  Maximum  velocities  at  depths  were  also  fairly  close,  differing  by  an  av¬ 
erage  of  50  m.  Velocity  magnitudes  were  comparable  for  the  surface  flows,  differing  by 
an  average  of  4  cm  s  '.  At  depth,  velocities  at  the  inflow  region  agreed  well  on  the 
southern  side  of  the  canyon,  with  differences  of  only  3.6  cm  s~'  between  the  two  meth¬ 
ods.  The  strong  geostrophic  outflow  region  at  depth  along  the  northern  slope  of  the 
earn  on  did  not  fare  as  well.  The  sense  of  flow  was  the  same  but  the  velocities  differed 
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by  almost  8  cm  srK  A  few  possible  explanations  for  this  discrepancy  in  fiow  at  this  lo¬ 
cation  follow: 

•  Anomalously  high  AGC  values  were  observed  during  the  processing  of  the  ADCP 
data  (see  Chapter  II)  over  this  region.  The  "scattering  layer",  cited  as  the  possible 
culprit  for  the  unusual  signal,  could  have  adversely  effected  ADCP  measurements 
below,  rendering  data  collected  in  this  area  "suspect." 

•  Standard  error  calculations  are  smaller  than  mean  values  of  both  geostrophic  and 
ADCP  flow,  with  the  exception  of  the  ADCP  mean  obtained  for  this  deep  outflow 
region  where  the  mean  is  -2.10  cm  s'1  and  the  standard  error  is  3.18  cm  s~l.  This 
is  another  reason  to  label  the  ADCP  data  obtained  in  this  area  as  "suspect". 

•  CTD  station  spacing  was  very  tight  (approximately  2  km  vice  4  for  the  other 
stations)  over  the  northern  canyon  shelf  whereas  ADCP  velocities  were  averaged 
every  4  km  over  the  whole  of  Line  1.  This  larger  averaging  interval  could  have  re¬ 
sulted  in  smoothing  the  center  of  the  outflow  core  and  resulted  in  a  much  smaller 
velocity. 


Table  3.  SUMMARY  OF  LINE  1  MEAN  VELOCITIES:  Negative  velocities  in¬ 
dicate  flow  out  of  the  bay,  positive  values  denote  flow  into  the  bay. _ 


Type 

Location  (km 
from  sta  1) 

Depth(dbar) 

Std  Err 
(cm  .r1) 

GHOSTROPH1C 

13.16 

0 

15.62 

9.95 

ADCP 

17.32 

0 

14.97 

1.72 

GEOSTROPHIC 

33.17 

0 

i 

GO 

4.81 

ADCP 

34.48 

0 

-10.48 

3.49 

GEOSTROPHIC 

16.78 

250 

-10.65 

8.3 

ADCP 

17.32 

182 

-2.10 

3.18 

GEOSTROPHIC 

25.89 

150 

10.15 

9.11 

ADCP 

29.99 

118 

13.71 

4.10 

D.  DESCRIPTION  OF  NOAA  AVHRR  SATELLITE  IMAGERY 

Cloud-free  Satellite  imagery  was  available  for  04,  08,  and  11  May  1988  (Figure  31 
on  page  57  and  Figure  32  on  page  58).  The  04  May  image  reflects  strong  coastal  up- 
welling  near  Pt.  Ano  Nuevo.  35  km  to  the  north  of  the  Monterey  Bay.  This  cold,  up- 
welled  water  (7  to  8’C)  appears  to  extend  into  Monterey  Bay  from  the  north.  Warm 
water  from  inside  the  bay  appears  to  be  advecting  southward  past  Pt.  Sur.  The  implied 
circulation  is  inflow  at  the  north  of  the  bay  and  outflow  at  the  south. 
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The  images  taken  on  08  and  1 1  May  are  nearly  identical  and  only  the  08  May  image 
is  shown.  Skin  temperatures  over  the  bay  showed  a  sharp  rise  in  temperature  of  ap¬ 
proximately  5'C  over  the  four  day  period  from  04  May.  Whether  this  is  a  result  of  warm 
water  advected  into  the  bay,  or  simply  a  product  of  solar  heating  which  ensued  upon  the 
cessation  of  winds  and  upwelling,  cannot  be  determined.  Coastal  upwelling  is  still  in 
evidence  to  the  north  and  south  of  the  bay  on  the  08  May  image.  Colder  waters,  either 
being  upwelled  along  the  coast  or  advected  into  the  northern  part  of  the  bay,  again  ap¬ 
pear  on  the  11  May  imagery. 

E.  METEOROLOGICAL  SYNOPSIS 

1.  The  Period  Prior  to  Data  Acquisition  (25  April  -  08  May  1988) 

On  25  April,  the  central  California  coast  was  dominated  by  a  1026  mb  high 
pressure  center  situated  over  the  center  of  the  state.  This  pattern  resulted  in  a  relatively 
light  pressure  gradient  overall,  allowing  the  mesoscale  land-sea  breeze  processes  to 
dominate  the  flow  pattern.  A  low  pressure  system  was  located  approximately  900  n  mi 
west  of  Washington  state.  By  the  26rt  the  high  began  to  break  down  under  the  influence 
of  this  low  pressure  system  and  by  the  27'*,  this  low,  with  its  associated  cold  front,  passed 
over  the  Monterey  Bay.  While  no  significant  weather  or  winds  were  produced  by  this 
system,  it  did  clear  the  path  for  a  second  cold  frontal  system  which  was  approximately 
600  n  mi  west  of  Monterey  Bay  at  this  time.  By  28  April,  the  pressure  field  in  the  vicinity 
of  the  Monterey  Bay  began  feeling  the  influence  of  this  system  generating  a  light  (3  to 
5  m  s“‘)  southwesterly  pre-frontal  wind  region  regime.  This  strong  frontal  system  passed 
over  the  Monterey  Bay  early  in  the  morning  of  29  April  causing  the  winds  to  rapidly  veer 
to  the  northwest  and  increase  in  magnitude  to  18  m  .r1,-  As  this  front  continued  east¬ 
ward  across  the  state,  strong  high  pressure  ridging  behind  this  front  moved  over  the 
Monterey  Bay  and  maintained  a  strong  pressure  which  continued  generating  13  to  18 
m  s~'  winds  over  the  area  through  the  morning  of  01  May.  From  the  afternoon  of  01 
May  through  07  May,  the  winds  remained  out  of  the  northwest  at  moderate  velocities 
(5-8  m  s’1)  increasing  slightly  the  morning  of  the  4:*  when  a  third  front,  weaker  than  the 
previous  one,  passed  over  the  Monterey  Bay.  In  summary,  this  period  was  very  "unset¬ 
tled"  for  May,  with  three  separate  frontal  systems  passing  over  the  bay.  Winds  were 
predominantly  out  of  the  northwest  to  north,  at  times  very  strong,  for  a  period  of  ap¬ 
proximately  nine  days,  from  29  April  to  07  May. 
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Figure  31. 


NOAA  AVHRR  satellite  imagery  from  2336  Z  04  May  1988:  Darker 
shades  along  the  coast  to  the  north  and  south  of  the  Monterey  Bay 
depict  colder  water.  Lighter  shades  over  the  extreme  eastern  portion 
of  the  bay  and  within  the  northern  and  southern  bights  of  the  bay  de¬ 
pict  wanner  water.  The  light  streak  extending  southwest  from  the 
Monterey  Peninsula  represents  warmer  waters.  The  large  areas  of  black 
over  the  ocean  are  cloud  cover.  Latitude  and  longitude  lines  are  in 
one-half  degree  increments.  Rectangle  at  left  is  SST  scale  meant  for 
use  with  color  copy  and  is  not  usable  with  this  black  and  white  re¬ 
production. 
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Figure  32.  NOAA  AVHRR  satellite  imagery  from  2247  Z  08  May  1988:  Darker 
shades  within  the  Monterey  Bay  and  to  the  west  of  the  bay  depict 
warmer  water.  Darker  areas  along  the  coast,  south  of  the  Monterey 
Peninsula,  depict  cooler  upwelled  water.  Black  areas  over  the  ocean  are 
clouds.  Latitude  and  longitude  lines  are  the  same  as  Figure  31.  Rec¬ 
tangle  at  bottom  is  SST  scale  meant  for  use  with  color  copy  and  is  not 
usable  with  this  biack  and  white  reproduction. 
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2.  Synoptic  Situation  During  the  Period  08  -  11  May  1988 

The  synoptic  meteorological  situation  during  the  period  of  data  acquisition  was 
high  pressure,  centered  to  the  southwest  of  the  bay,  which  was  strong  enough  during  this 
period  to  block  the  advance  of  any  atmospheric  low  pressure  systems  or  cold  fronts. 
The  pressure  gradient  was  not  strong  enough  to  influence  the  mesoscale  land-sea  breezes 
which  dominated  this  entire  period  (see  Figure  21  on  page  37).  Calm  conditions  pre¬ 
vailed.  with  very  low  (less  than  1  m)  sea  states  during  the  entire  survey  cruise. 

3.  The  Monterey  Bay  Aquarium  Wind  Field 

Review  of  the  MBA  wind  data  taken  over  a  3  month  period,  from  the  last  week 
of  April  to  the  last  week  in  June  1988,  shows  the  predominant  diurnal  land-sea  breeze 
fluctuations  evident  during  periods  of  light  winds.  The  wind  events  previously  men¬ 
tioned  during  the  last  part  of  April  and  continuing  to  approximately  07  May  is  evident 
in  this  time  series  (Figure  21  on  page  37).  During  this  period  winds  were  predominantly 
out  of  the  northwest  and  attained  velocities  of  up  to  15  m  s~‘.  The  overall  mean  wind 
speed  for  the  three  month  period  was  approximately  3.5  m  ,r>  from  2607'.  Calculated 
mean  wind  speed  and  direction  for  the  period  of  the  May  19SS  cruise  was  2.3  m  from 
260T  .  Also  of  note  is  the  predominantly  stronger  winds  during  afternoon  seabreeze 
events,  with  winds  dying  off  rapidly  during  the  evening  hours. 

F.  SUMMARY  OF  MEAN  CURRENT  FLOW 

Surface  How  implied  through  satellite  imagery,  geostrophic  flow,  ADCP  flow,  and 
dynamic  height  analyses  were  consistent.  They  show  mean  surface  inflow  over  the 
northern  half  of  the  bay  and  stronger,  more  concentrated  outflow  over  the  southern  half 
of  the  bay.  At  depth,  geostrophic  and  ADCF  data  were  consistent  on  the  southern  side 
of  the  canyon  and  show  inflow  on  the  order  of  10  cm  s-1.  Because  of  the  possible  prob¬ 
lems  noted  in  the  acquisition  of  ADCP  data  along  the  northern  wall  of  the  canyon,  the 
geostrophic  mean  is  believed  to  be  more  reliable  over  this  area.  For  this  reason,  and 
because  the  geostrophic  and  ADCP  mean  data  were  similar  in  the  other  areas  of  the  bay, 
the  geostrophic  mean  field  was  used  to  force  the  numerical  ocean  model  of  the  Monterey 
Bay. 
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G.  NUMERICAL  OCEAN  MODELING  OF  THE  MONTEREY  BAY 
CIRCULATION 


1.  Model  Description 
a.  Dynamics 

The  model  used  in  this  study  is  a  two  layer,  semi-implicit,  primitive  equation 
model  which  uses  a  63  by  37  element  C-grid  for  finite  differencing.  The  numerical 
scheme  utilized  in  this  model  was  initially  derived  by  Hurlburt  (1974)  to  study  ocean 
mesoscale  circulation  features  (Bruner,  1988)  and  later  adapted  for  use  in  studying  cir¬ 
culation  processes  within  the  Monterey  Bay  by  Prof.  David  C.  Smith  IV  of  the  Naval 
Postgraduate  School,  Monterey,  CA.  Linear  test  cases  with  analytic  solutions  have  been 
run  in  tandem  by  Smith  and  Reid  (1982)  to  validate  the  model  (Bruner,  1988).  The 
equations  used  in  the  model  are  the  vertically  integrated  momentum  equation: 


ci 


+  (V.  v,+  =  +  6pn„  (12) 


and  a  continuity  equation: 
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for  layer  (i=  1  upper  and  i  =  2  lower)  thickness  h„  transports  l]  and  velocities  v,.  The 
relation: 


P\=g{h\+h2+d)  (14) 

was  used  to  calculate  pressure  in  the  upper  layer,  and: 

r2  =  r\-g'h\  (i5) 

was  used  for  lower  layer  pressure  calculations  where  the  reduced  gravity  ( g ')  is  calculated 
using: 


g  =g- 


(P2-P1) 


Pi 


(16) 


The  Coriolis  parameter  if)  is  equal  to  2£2  sin  </>  where  Q  is  the  angular  velocity  of  the 
earth  in  rad  s  l  and  <b  is  the  latitude.  The  Kronccker  Delta  (()„)  is  equal  to  one  when  i  =  j 
and  zero  when  i  ¥=  j.  Subgrid  scale  dissipation  is  modeled  using  biharmonic  (Va)  friction 
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with  the  diffusivity  constant  ( B„ )  equal  to  -  1  x  107  m4  s_1.  Wind  stress  between  the  air 
and  water  (tw)  is  included  as  a  forcing  parameter  along  with  bottom  friction  which  em¬ 
ploys  a  linear  bottom  drag  coefficient  (r)  equal  to  5  x  10~4ms_l  (Winant  and  Beardsley, 
1979)  in  the  relation  -rv2  (Shaw  and  Csanady,  1983).  The  reduced  gravity  ( g ')  value  used 
was  2  x  10‘2  mrl 

b.  Boundary  Conditions  and  Domain 

The  north,  east,  and  south  boundaries  are  closed  (no  slip).  The  west 
boundary  is  open  (following  Camerlengo  and  O'Brien,  1980)  except  where  forcing  (in¬ 
flow  and  outflow)  is  specified.  Forcing  along  the  western  boundary  can  be  applied  sep¬ 
arately  to  the  upper  (w'hich  encompasses  the  free  surface  to  50  m)  and  lower  (50  m  to 
the  bottom)  layers,  allowing  both  barotropic  and  baroclinic  modes. 

The  63  by  37  grid  is  spaced  at  500  m  intervals  resulting  in  a  model  domain 
encompassing  a  32  by  18  km  rectangle.  The  domain  is  orientated  roughly  north- 
northwest  to  south-southeast  (see  Figure  33  on  page  62)  to  encompass  as  much  area 
of  the  Monterey  Bay  as  possible.  Realistic  topography  (Figure  33  on  page  62)  was  in¬ 
corporated  by  applying  an  interpolated  field  of  gridded  bathymetry.  Because  of  numer¬ 
ical  considerations,  the  layer  interface  cannot  intersect  either  the  free  surface  or  bottom 
topography.  For  this  reason,  bottom  topography  shallower  than  100  m  was  not  allowed 
so  unrestricted  vertical  displacement  of  the  interface  could  occur  (Bruner,  1988). 

2.  Initialization  Methods 

a .  Methods  of  Model  Forcing 

The  first  problem  encountered  was  whether  to  use  the  mean  flow  obtained 
through  geostrophy  or  through  ADCP  measurements.  Since  the  sense  and  magnitude 
of  the  two  means  wrere  similar  and,  because  of  the  doubtful  reliability  of  the  ADCP  data 
at  depth  near  the  northern  wall  of  the  canyon  (discussed  at  length  in  the  ADCP  section), 
it  was  decided  to  use  the  mean  obtained  through  geostrophy.  Aside  from  that  one  re¬ 
gion,  they  were  not  significantly  different.  Since  the  model  is  a  two  layer  system,  the 
mean  geostrophic  flow  was  averaged  over  the  top  50  m  to  obtain  v,  values,  then  averaged 
between  50  m  and  the  bottom  depth  between  each  station  (Figure  34  on  page  63)  to 
obtain  values  for  v2.  This  averaging  was  performed  in  order  to  conserve  the  mean  mass 
transport  of  the  mean  geostrophic  flow  field.  This  procedure  left  the  surface  velocities 
relatively  unscathed,  but  resulted  in  lowered  velocities  at  depth  since  most  of  the  mean 
deep  flow  lay  above  250  m.  The  next  obstacle  was  how  to  force  the  model  with  known 
data  at  inconsistent  intervals  and  not  in  coincidence  with  the  western  boundary  of  the 
model  (see  Figure  33  on  page  62).  Since  both  Line  1  and  the  western  boundary  of  the 
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Figure  33.  Numerical  model  bathymetry  (top)  and  spatial  domain  (bottom):  Top 
figure  shows  the  bathymetry  used  in  the  Monterey  Bay  numerical 
model,  bottom  figure  shows  the  domain  (box)  and  its  orientation  within 
the  bay.  Line  1  (stations  1-10;  and  Line  2  (stations  10-17)  are  also  de¬ 
picted  in  relation  to  the  model  domain.  Arrows  roughly  approximate 
the  method  used  in  projecting  the  mean  geostrophic  flow  onto  the 
model's  western  (open)  boundary. 
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Figure  34.  Cross  sections  of  model  inputs:  Top  (left)  figure  shows  the  result  of 
vertically  averaging  the  mean  geostrophic  How  along  Line  1.  Top 
(right)  figure  depicts  the  projection  of  these  averaged  flows  onto  the 
western  boundary  of  the  model.  The  bottom  figure  represents  the  mass 
balanced  flow  input  into  the  model  during  subsequent  experiments. 
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model  are  orientated  roughly  perpendicular  to  bottom  slope,  the  simplest  solution  was 
to  project  the  mean  flow  at  known  points  onto  the  western  boundary  of  the  model  grid 
(see  Figure  33  on  page  62)  along  isobaths.  This  of  course  is  a  realistic  assumption  if 
vorticity  is  conserved  over  the  mouth  of  the  bay  and  if  the  predominant  flow  was  per¬ 
pendicular  to  Line  1.  Once  the  velocities  were  identified  on  the  model  boundary,  the 
data  were  run  through  a  linear  interpolation  scheme  to  define  velocity  values  at  each  grid 
point  along  the  boundary  (see  Figure  34  on  page  63  and  Appendix  C).  Runs  were  then 
made  using  this  farcing  scheme,  both  with  and  without  the  mean  wind  and  bottom  fric¬ 
tion  (Table  4  on  page  65).  These  model  runs  ran  approximately  8  days  before  the 
interface  intersected  the  bottom.  When  this  occurred  the  runs  terminated  because  model 
assumptions  were  violated.  This  was  not  surprising  since  the  original  mean  geostrophic 
flow  was  not  exactly  mass  balanced  (positive  net  transport,  or  flow  into  the  bay),  thus 
surplus  water  was  continually  being  pumped  into  the  bay  at  each  time  step  with  no 
outlet  available.  In  an  attempt  to  stabilize  the  model,  the  input  mass  field  (wiJ  s~l)  ori¬ 
ginally  derived  from  the  geostrophic  mean  flow  was  adjusted  so  that  net  transport  in 
both  layers  added  up  to  zero.  This  was  accomplished  by  applying  linear  correction  co¬ 
efficients  to  both  layers  (Appendix  C).  This  adjustment  of  mass  balance  is  considered 
justifiable  and  realistic  for  the  following  reasons: 

•  The  "true"  mean  flow  would  be  mass  balanced. 

•  The  amount  of  adjustment  applied  to  the  boundary  forcing  is  probably  less  than 
the  geostrophic  resolution. 

•  As  a  consequence  of  the  model's  two  layered  system,  averaging  of  the  mean  flow 
was  only  conducted  down  to  50  m  for  the  upper  layer.  Inspection  of  the  mean 
geostrophic  flow'  field  (Figure  24  on  page  42)  shows  that,  while  it  definitely  appears 
to  be  an  uncoupled,  two  layer  system,  the  top  layer  most  likely  extends  deeper  than 
50  m.  Therefore,  much  of  the  flow  in  the  low'er  layer  wras  artificially  introduced 
during  the  partitioning  process. 

Runs  with  the  balanced  flow'  field  appeared  more  stable  numerically,  and  ran  out  to  an 
average  of  19  days. 

To  see  what  the  model  response  would  be  to  wind  forcing  alone,  a  run  was 
made  simulating  the  actual  winds  over  the  bay  during  a  strong  frontal  passage  (similar 
to  what  was  described  in  the  Meteorological  Synoptic  Situation  section).  Input  winds 
were  first  "ramped"  (i.e.  gradually  increased  each  time  step  to  avoid  numerical  instabil¬ 
ity)  from  0  to  2  m  r1  from  the  south  to  southeast.  Then,  from  da>  3  to  day  7,  the  winds 
were  ramped  from  southerly  to  southwesterly,  and  the  winds  were  ramped  up  to  8.5  m 
s~l  to  simulate  pre-frontal  winds.  From  day  7  to  day  8,  winds  were  kept  strong,  but  di- 
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rection  was  shifted  each  time  step  from  southwest  to  northwest  to  represent  frontal 
passage.  Winds  were  gradually  decreased  from  8.0  mr1  to  1.0  m  s~\  while  direction 
was  backed  to  westerly,  over  the  next  two  days  to  day  10.  Past  day  10  the  winds  were 
turned  off  so  the  "spun  down"  state  could  be  observed.  Both  layers  were  initially  at  rest 
and  the  entire  western  boundary  was  left  open  (allowing  inflow  and  outflow  anywhere). 

All  model  runs  were  made  with  bottom  topography  included  (for 
topographic  effects  see  Bruner,  1988).  Open  boundary  conditions  were  specified  wher¬ 
ever  mean  forcing  was  determined  to  be  zero  along  the  boundary  to  allow  inflow  or 
outflow  at  that  point.  As  previously  mentioned,  the  entire  boundary  was  open  for  the 
wind  forcing  experiments.  The  mean  MBA  wind  speed  and  direction  was  used  in  those 
mean  forcing  runs  which  applied  a  constant  wind  stress  throughout.  A  summary  of  the 
five  runs  conducted  is  provided  in  Table  4  (outputs  appear  in  Figure  35  on  page  66  to 
Figure  39  on  page  70). 


Table  4.  SUMMARY  OF  MONTEREY  BAY  NUMERICAL  MODEL  SIMU¬ 
LATIONS 


Run  No. 

Boundary'  Forcing 

Bottom  Friction 

Wind  Stress 

1 

Geostrophic  Mean 

No 

No 

2 

Geostrophic  Mean 

Yes 

260T  at  2.3  m  r1 

3 

Mass  Balanced 
Geostrophic  Mean 

No 

No 

4 

Mass  Balanced 
Geostrophic  Mean 

Yes 

mam 

5 

None 

No 

Variable  (see  text) 

3.  Results  of  Experiments 

Results  of  run  number  one  (Figure  35  on  page  66)  shows  anticyclonic  rotation 
over  the  southern  half  of  the  bay  with  a  closed  cyclonic  eddy  situated  in  the  northern 
bight  region  in  the  surface  layer.  At  depth,  velocities  are  weakly  cyclonic  (less  than  2.5 
cm  i_1)  over  the  mouth  of  the  canyon  with  weak  anticyclonic  curvature  in  the  northern 
and  southern  bights.  The  addition  of  mean  winds  and  bottom  friction  (run  two)  did  not 
seem  to  have  a  marked  effect  on  the  overall  flow  (Figure  36  on  page  67).  Flow  patterns 
w'ere  similar  in  both  cases,  with  the  exception  that,  as  expected,  velocities  w'ere  decreased 
in  run  two  by  the  effect  of  bottom  friction.  To  gain  a  quantitative  comparison  with  the 
mean  flow'  calculated  across  Line  2,  geostrophic  velocities  w’ere  averaged  down  to  50  m 
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Figure  35.  Monterey  Bay  model  output  from  run  No.  1:  North  is  to  the  left  in 
each  figure  (i.e.  figures  are  orientated  as  if  you  are  looking  into  the 
bay).  Arrows  are  velocity  vectors,  solid  lines  are  isotachs  (cm  .r1).  Top 
figure  depicts  upper  level  flow  and  bottom  figure  shows  deep  flow,  both 
at  day  7  of  the  run.  Mean  gcostrophic  flow  was  used  to  force  the  model 
along  the  western  (bottom  of  plots)  houndary.  The  effects  of  wind  and 
bottom  friction  are  neglected. 
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Figure  36.  Monterey  Bay  model  output  from  run  No.  2:  North  is  to  the  left  in 
each  figure  (i.c.  figures  are  orientated  as  if  you  are  looking  into  the 
bay).  Arrows  arc  velocity  vectors,  solid  lines  are  isotachs  (cm  s-1).  Top 
figure  depicts  upper  level  flow  and  bottom  figure  shows  deep  flow,  both 
at  day  7.25  of  the  run.  Mean  geostrophic  flow  was  used  to  force  the 
model  along  the  western  (bottom  of  plots)  boundary.  The  effects  of 
wind  and  bottom  friction  arc  included. 
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Figure  37.  Monterey  Bay  model  output  from  run  No.  3:  North  is  to  the  left  in 
each  figure  (i.c.  figures  are  orientated  as  if  you  are  looking  into  the 
bay).  Arrows  are  velocity  vectors,  solid  lines  are  isotachs  (cm  .r1).  Top 
figure  depicts  upper  level  flow  and  bottom  figure  shows  deep  flow,  both 
at  day  15.8  of  the  run.  Mass  balanced  mean  gcostrophic  flow  was  used 
to  force  the  mode!  along  the  western  (bottom  of  plots)  boundary.  The 
effects  of  wind  and  bottom  friction  arc  neglected. 
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Figure  38.  Monterey  Bay  model  output  from  run  Nr.  4:  North  is  to  the  left  in 
each  figure  (i.c.  figures  arc  orientated  as  if  you  are  looking  into  the 
bav).  Arrows  arc  velocity  vectors,  solid  lines  are  isotachs  (cm  Top 
figure  depicts  upper  level  flow  and  bottom  figure  shows  deep  flow,  both 
at  day  15.8  of  the  run.  Mass  balanced  mean  gcostrophic  flow  was  used 
to  force  the  model  along  the  western  (bottom  of  plots)  boundary.  I  he 
effects  of  wind  and  bottom  faction  arc  included. 
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Figure  39.  Monterey  Bay  model  output  from  run  No.  5:  North  is  to  the  left  in 
each  figure  (i.e.  figures  ate  orientated  as  if  you  arc  looking  into  the 
bay).  Arrows  are  velocity  vectors,  solid  lines  are  isotachs  (cm  .r1).  Top 
figure  depicts  upper  level  flow  and  bottom  figure  shows  deep  flow,  both 
at  day  13.6  of  the  run.  Wind  forcing  alone  was  used  to  force  flow  in 
the  model.  Both  layers  were  initially  at  rest.  The  entire  western 
boundary'  (bottom  of  plots)  was  open  to  allow  flow  to  develop  at  will. 
The  effect  of  bottom  friction  was  neglected. 


and  from  50  m  to  the  bottom  and  compared  with  model  outputs  at  these  grid  points 
(Table  5  on  page  71,  see  also  Figure  33  on  page  62  for  station  positions).  The  mean 
Line  2  data  were  not  used  to  force  the  model  in  any  way.  The  mean  geostrophic  flow 
between  stations  17  and  16  does  show  light  outflow  in  the  extreme  northern  part  of  the 
bay,  although  not  to  the  extent  the  model  shows.  Outflow  is  not  indicated  in  the  mean 
Line  1  geostrophic  velocities  implying  that  some  sort  of  weak,  closed  circulation  may, 
in  fact,  exist  in  the  northern  bight  area.  Model  output  across  Line  2  also  shows  strong 
inflow  in  the  lower  layer  over  the  northern  part  of  the  bay  and  outflow  in  the  southern 
part.  Line  2  mean  velocities  show  no  such  flow  in  these  areas.  In  fact,  the  mean  Line  2 
geostrophic  flow  shows  very  little  flow  in  the  lower  layer  at  all,  with  1.94  cm  s~l  inflow 
between  stations  15  and  14,  and  -0.99  cm  s_1  outflow  between  stations  14  and  13.  In 
general,  the  overall  flow  patterns  predicted  by  the  bay  model  are  much  stronger,  in  both 
upper  and  lower  layers,  than  those  found  in  the  observed  mean  flow  field.  Circulation 
in  the  northern  bight  which  was  generated  in  runs  1  and  2  can  be  explained  in  terms  of 
compression  and  expansion  of  vortex  tubes.  As  water  is  pumped  into  the  northern  part 
of  the  bay  (too  much  water)  in  the  upper  layer,  it  travels  into  the  northern  bight  region 
where  it  forces  the  interface  down  producing  vortex  stretching  in  the  upper  layer  and 
compaction  of  tubes  in  the  lower  layer  resulting  in  cyclonic  circulation  in  layer  one,  and 
anticyclonic  circulation  in  layer  two. 


Table  5.  COMPARISON  OF  MODEL  OUTPUT  AND  LINE  2  GEOSTROPHIC 


MEAN:  Negative  velocities  indicate  flow  out  of  the  bay,  positive  values 
denote  flow  into  the  bay.  Velocities  are  in  cm  s~'. _  _ 


Station 

Observed  Mean 

Run  2 

Run  4 

Run  5 

v'i 

V*2 

0 

v. 

Vj 

17-16 

-0.7 

0 

-5 

9 

1 

0 

4 

10 

16-15 

-0.3 

0 

2 

3 

3 

0 

3 

7 

15-14 

4 

2 

10 

0 

4 

0 

1 

-.5 

14-13 

4 

-1 

6 

0 

5 

-.5 

-.6 

-.5 

13-12 

4 

0 

5 

1 

3 

1 

-.5 

-2 

12-11 

2 

0 

4 

1 

1 

-.5 

-1 

12 

11-10 

-9 

0 

-1 

0 

-3 

-.5 

-1 

2 

The  effect  of  balancing  the  mass  transport  along  the  western  boundary 
produced  major  changes  in  the  flow  pattern  (Figure  37  on  page  68).  The  circulation  in 
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the  northern  bight  produced  by  run  one  and  run  two  is  no  longer  evident  in  runs  three 
and  four  (Figure  37  on  page  68  and  Figure  38  on  page  69).  The  surface  layer  consists 
of  a  broad  anticyclone  encompassing  practically  the  entire  bay,  and  centered  over  the 
canyon  axis.  In  run  three,  flow  in  the  lower  layer  generally  follows  topography  with 
weak  cyclonic  circulation  over  the  canyon  mouth,  somewhat  stronger  anticylonic  rota¬ 
tion  over  the  northern  shelf,  and  anticyclonic  circulation  contained  in  the  southern  bight 
region.  Flow  appears  to  diverge  from  the  cyclone  in  the  canyon's  mouth,  travel  deeper 
into  the  bay  along  the  canyon  axis  (east),  then  enters  the  southern  bight  via  the  canyon 
head  where  it  ties  into  the  previously  mentioned  lower  layer  anticylone  located  in  this 
area.  The  addition  of  bottom  friction  effectively  keeps  the  flow  trapped  in  the  canyon 
(Figure  38  on  page  69).  Magnitude  and  direction  of  flow  in  both  layers  from  run  four 
compares  favorably  with  the  observed  mean  geostrophic  flow  (Table  5  on  page  71)  with 
the  exception  of  not  predicting  weak  outflow  in  layer  one  over  the  extreme  northern  part 
of  the  bay.  The  model  also  appears  to  under  predict  the  outflow  in  the  southern  part 
of  the  bay  in  layer  one  (-3  vice  -9  cm  s~'  which  was  observed).  A  Silicon  Graphics 
"movie"  of  this  run,  made  at  1000  s  intervals,  showed  that  the  model  was  highly  stable 
using  this  balanced  forcing.  The  predicted  high  dynamic  height  center  remained  quasi¬ 
stationary  over  the  canyon  mouth  throughout  the  19  day  animation,  agreeing  well  with 
the  inferred  flow  from  the  actual  mean  dynamic  height  analyses  (Figure  28  on  page 
49).  Again,  the  model  realistically  predicts  the  observed  Line  2  mean  flow,  especially  in 
the  lower  layer,  which  supports  the  hypothesis  that  much  of  the  cyclonic  circulation 
within  the  canyon  at  depth  does  not  penetrate  into  the  canyon  to  any  appreciable  extent, 
but  rather  recirculates  near  the  mouth  of  the  canyon. 

The  final  model  simulation  (Figure  39  on  page  70),  which  used  variable  wind 
stress  as  the  sole  forcing  parameter,  showed  that  strong  atmospheric  frontal  passages, 
such  as  the  type  previously  described,  can  set  up  circulation  patterns  which  are  similar 
in  nature,  at  least  in  the  upper  layer,  to  those  observed  in  the  mean  geostrophic  flow 
(Figure  39  on  page  70).  The  upper  layer  outflow  aiong  the  southern  portion  of  the  open 
boundary  is  forced  to  the  southwest,  similar  to  what  was  depicted  in  the  04  May  satellite 
imagery.  Reproducing  flow  similar  to  that  observed  in  the  lower  layer  mean  geostrophic 
flow  using  this  method  was  not  as  successful.  The  model  output  shows  anticyclonic 
turning  in  the  northern  shelf  area  and  cyclonic  turning  in  the  lower  layer  over  the 
southern  shelf  break.  These  flows  then  converge  near  the  canyon  head  and  flow  out  to 
sea  along  the  axis.  One  major  shortcoming  of  using  this  forcing  method  is  the  fact  that 
theie  is  no  model  topography  above  sea  level  to  produce  lee  effects  such  as  in  the 
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northern  and  southern  bight  regions  which,  in  reality,  are  fairly  protected  from 
northwesterly  winds  (Breaker  and  Broenkow,  1989).  Animation  of  this  simulation 
showed  height  anomalies  to  vary  widely  in  location  and  intensity.  This  run  never  quite 
reached  a  stable  state,  even  after  day  10,  when  wind  stress  was  turned  off  implying  wind 
alone  is  not  the  only  factor. 


IV.  DISCUSSION 


A.  COASTAL  PROCESSES 

Strong  northwesterly  winds,  produced  by  post-frontal  high  pressure  ridging  across 
California,  prevailed  over  the  central  California  coast  from  the  last  week  of  April 
through  the  first  week  of  May  1988.  The  oceanic  response  to  this  wind  stress  was 
coastal  upwelling,  produced  by  Ekman  transport,  and  the  appearance  of  a  strong 
equatorward  geostrophic  jet  along  the  central  coast.  NOAA  satellite  imagery  from  04 
May  1988  clearly  shows  this  upwelling  and  the  presence  of  equatorward  flow,  entering 
the  Monterey  Bay  from  the  north,  and  advecting  warm  Bay  waters  south  past  Pt.  Sur, 
37  km  to  the  south  of  the  Monterey  Bay.  This  pattern  is  similar  to  observations  made 
by  Breaker  and  Broenkow  (1989),  who  noted  that,  during  upwelling  events,  cold,  up- 
welled  water  was  advected  across  the  bay,  producing  a  strong  thermal  front  near  Pt. 
Pinos.  A  geostrophic  velocity  cross  section,  produced  from  hydrographic  data  collected 
along  a  line  perpendicular  to  Pt.  Sur  on  04  May  1988  depicts  a  strong  25  cm  s~l 
equatorward  flow  extending  from  the  coast  to  50  km  seaward  and  limited  to  pressures 
of  less  than  100  dbar,  typical  of  coastal  upwelling  geostrophic  jets  (Figure  40  on  page 
75).  These  descriptions  of  the  winds,  upwelling,  and  coastal  geostrophic  flow  are  all 
consistent  with  the  findings  and  observations  of  Halliwell  and  Allen  (1987),  Buyer 
(1983),  and  Hickey  (1979).  Of  final  note,  the  CUC  was  not  evident  on  the  04  May 
geostrophic  cross  section  at  Pt.  Sur,  consistent  with  Chclton's  (1984)  observations  that 
this  poleward,  deep  current  was  weak  and  equatorward  during  April  and  May  at 
Monterey  and  Pt.  Sur. 

The  08-11  May  1988  Monterey  Bay  hydrographic  survey  commenced  at  the  cessa¬ 
tion  of  strong  northwesterly  winds.  Winds  were  light  and  seas  were  calm  throughout  the 
period.  As  a  result  of  this  period  of  calm  winds,  satellite  imagery  taken  on  08  and  1 1 
May  1988  showed  waters  within  the  bay  to  have  w'armed  significantly,  either  by 
advection  of  warmer  oceanic  waters  into  tne  bay  (Breaker  and  Broenkow,  1989),  or 
simply  through  solar  radiative  processes.  While  no  quantitative  evidence  exists  to  sup¬ 
port  the  continued  presence  of  the  equatorward  coastal  current  during  this  period,  weak 
upwelling  was  still  visible  to  the  north  and  south  of  the  Monterey  Bay. 
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Figure  40.  Pt.  Sur  geostropliic  velocity  cross-section:  Geostrophic  Velocities  were 
derived  from  C  I'D  casts  collected  on  04  and  05  May  1988.  The  section 


runs  roughly  east-west,  perpendicular  to  the  coast,  from  Pt.  Sur 
seaward.  Pt.  Sur  is  located  roughly  20  n  mi  (37  km)  south  of  the 
Monterey  Peninsula.  Positive  values  denote  poleward  flow,  negative 
values  denote  cquatorward  How. 


B.  INTERNAL  WAVES 

The  mean  flow  signal,  especial!}  over  the  shelf  break  and  the  Monterey  Submarine 
Cans  on,  was  highly  aliased  by  the  internai  wave  field  throughout  the  data  acquisition 
period  (Figure  19  on  page  34).  'lire  process  of  averaging  the  mass  field  and  ADCP  ve¬ 
locities  appears  to  have  effectively  removed  the  majority  of  tins  'noise.  Examination 
of  the  ADCP  passes  along  Line  1  showed  that  surface  tidal  heights  and  the  flow  within 
the  canyon  were  not  in  phase.  The  internal  wave  amplitudes  of  30  m,  observed  in  Line 
1,  aic  smaller  than  the  50-120  m  amplitudes  observed  by  Shea  and  Biocnkow  (1982)  near 


the  canyon  head.  This  difference  can  be  explained  by  their  hypothesis  that  internal  wave 
energy  is  focused  near  the  head  as  they  shoal,  thus  producing  the  somewhat  higher  am¬ 
plitudes  they  observed.  This  is  also  consistent  with  Broenkow  and  McKain's  (1972) 
observations  that  internal  tidal  waves  generally  decreased  seaward.  Why  observed 
internal  wave  amplitudes  observed  in  Line  2  data,  which  is  inshore  of  Line  1,  were 
smaller  than  those  observed  in  Line  1  is  unknown. 

C.  CIRCULATION 

The  mean  flow  during  the  data  acquisition  period  is  an  uncoupled  two  layer  system, 
consisting  of  a  warm  core  anticyclone  at  the  surface  over  a  cold  core,  cyclone  which  is 
"trapped"  within  the  confines  of  the  Monterey  Submarine  Canyon  under  the 
thermocline.  Inflow  at  the  surface  is  broad,  extending  from  the  northern  part  of  the  bay 
south  to  the  center  of  the  canyon,  while  outflow  is  relatively  stronger  but  narrower  and 
is  confined  to  the  area  between  the  Monterey  Peninsula  and  the  center  of  the  canyon 
(station  7).  Mean  flow  at  depths  of  100  to  200  m  consists  of  two  cores  with  outflow 
along  the  northern  wall  of  the  canyon  and  inflow  along  the  southern  half  of  the  canyon 
axis.  Both  ADCP  and  geostrophic  means  (Figure  30  on  page  53  and  Figure  24  on  page 
42),  along  with  dynamic  topographies  (Figure  28  on  page  49),  support  this  description. 
The  only  major  discrepancy  noted  between  the  ADCP  and  the  geostrophic  means  was 
in  the  outflow  region  along  the  north  canyon  wall. 

Evaluation  of  the  mean  flow  pattern  in  Line  2  also  shows  surface  inflow  over  the 
center  of  the  bay  with  outflow  along  the  southern  portion.  Surface  velocities  are  gen¬ 
erally  weaker  in  Line  2  and  no  evidence  of  cyclonic  circulation  at  depth  is  evident,  just 
a  weak  inflow  near  the  center  of  the  canyon.  This  leads  to  the  conclusion  that  the  sur¬ 
face  flow  is  turning  over  Line  2,  thus  accounting  for  the  small  east-west  velocity  com¬ 
ponents,  and  that  the  cyclonic  flow  noted  in  the  canyon  itself  is  limited  to  the  deeper 
portions  of  the  canyon,  near  the  mouth. 

As  seen  in  the  cross-sections  of  mean  temperature  and  density,  as  well  as  in  the 
mean  10'C  isotherm  plots,  the  isotherms  and  isopycnals  had  a  much  more  pronounced 
slope  from  the  canyon  axis  to  the  Monterey  Peninsula  explaining  the  higher  geostrophic 
and  ADCP  velocities  found  there.  Map  contours  of  mean  temperature  and  dynamic 
heights  support  these  flow  patterns  at  the  surface  and  at  depth.  The  tighter  dynamic 
height  gradient  over  the  southern  portion  of  the  bay  is  also  indicative  of  stronger 
geostrophic  flow  in  this  area.  Satellite  imagery,  especially  the  04  May  1988  pass,  sup 
ports  anticyclonic  surface  flow.  The  surface  flow  pattern  agrees  with  the  findings  of 
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Bigelow  and  Leslie  (1930),  Skogsberg  (1936),  Lammcrs  (1971),  and  with  Pirie  and 
Stellar's  (1974)  inferred  southward  flow  from  March  through  July.  The  noted  reversal 
in  flow  confirms  Skogsberg's  (1936)  suspicions  that  flows  within  the  bay  at  different 
depths  may  be  in  opposition  to  one  another.  Lasley  (1977),  Smethie  (1973),  Broenkow 
and  Smethie  (1978),  Reise  (1973),  Moomy  (1973),  and  Breaker  and  Broenkow  (1989) 
determined  flow  to  be  predominantly  cyclonic  within  the  bay  but  most  acknowledged 
that  periods  of  flow  reversals  do  occur  within  the  bay.  During  the  period  of  this  study, 
one  of  these  "reversals"  was  obviously  taking  place. 

As  previously  mentioned  the  CUC  was  not  in  evidence,  at  least  at  Pt.  Sur,  during 
the  period  so  momentum  transfer  from  this  current  does  not  seem  a  viable  flow  gener¬ 
ation  mechanism  (Breaker  and  Broenkow,  1989)  during  the  May  1988  cruise.  The  best 
explanation  for  the  observed  flow  patterns  stems  from  the  results  of  Klinck's  (1989) 
three-layer  geostrophic  adjustment  model  of  flow  over  a  canyon.  Despite  the  obvious 
differences  between  modeled  shelf  and  canyon  geometry',  and  the  actual  bathymetry  of 
the  Monterey  Canyon  and  vicinity,  there  are  similarities  between  his  results  and  the  flow 
patterns  observed  in  the  Monterey  Bay.  His  model  results  showed  that,  for  a  narrow 
canyon  (shown  to  be  the  case  here  in  Chapter  III),  anticyclonic  surface  flow  over 
cyclonic,  canyon  trapped  flow  would  be  produced  by  the  geostrophic  adjustment  process 
that  takes  place  when  a  surface  current  encounters  a  submarine  canyon.  Similarities 
between  the  50  km  wide  current  used  in  forcing  his  model,  and  the  50  km  wide, 
equatorvvard  surface  current  depicted  in  the  Pt.  Sur  geostrophic  velocity  cross-section 
(Figure  40  on  page  75)  is  also  noted.  Klinck's  results  also  indicated  that  once  a  canyon 
width  narrows  to  less  than  half  of  the  internal  Rossby  radius,  it  has  little  perturbing  ef¬ 
fect  on  the  cross  canyon  flow.  The  calculated  internal  Rossby  radius  at  station  13  was 
1 1.4  km  and  the  canyon  width  in  this  area  is  approximately  6  to  7  km.  This  may,  in  part, 
explain  the  absence  of  any  significant  deep  flow  at  depth  in  Line  2,  where  the  Monterey 
Submarine  Canyon  begins  to  narrow  rapidly. 

Whether  or  not  flow  reversals  within  the  Monterey  Bay  are  correlated  to  strong 
upwelling  events  cannot  be  determined  from  the  results  of  this  study  alone.  The  findings 
of  Broenkow  and  Smethie  (1978)  implied  predominantly  northward  flow  within  the  bay 
even  during  the  spring  and  summer  upwelling  period.  This  may  suggest  that  the  observed 
mean  flow  pattern  may  only  develop  during  periods  of  strong  equatorward,  coastal 
flows.  The  minimum  velocity  required  to  drive  such  flows,  the  amount  of  time  this 
coastal  current  requires  to  "spin-up"  the  observed  flow,  or  the  length  of  time  this  flow' 
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pattern  remains  after  cessation  of  this  current,  remain  unknown  quantities  and  are  left 
as  topics  of  future  research. 

D.  MEAN  FLOW  MODELING 

The  general  flow  pattern  described  above  was  also  reflected  in  numerical  model 
simulations,  particularly  when  forced  by  a  balanced  flow  (runs  3  and  4).  The  differences 
between  the  unbalanced  (Figure  36  on  page  67)  and  balanced  (Figure  38  on  page  69) 
runs  can  be  seen  by  analyzing  model  outputs  of  height  anomalies  (Figure  41  on  page 
80).  Run  2  height  anomalies  showed  surface  heights  of  approximately  7  cm  over  the 
mouth  of  the  canyon  with  a  relative  minimum  of  5.6  cm  over  the  northern  bight  resulting 
in  anticyclonic  and  cyclonic  circulations  respectively  in  the  upper  layer.  The  interface 
from  run  2  showed  a  depression  of  28  m  over  the  northern  bight  but  38  m  over  the 
canyon  mouth  creating  relative  highs  and  lows  respectively.  Analysis  of  Run  4  height 
anomalies  only  showed  one  high  in  the  upper  layer  and  one  low  in  the  lower  layer,  both 
vertically  stacked  and  situated  roughly  over  the  mouth  of  the  Monterey  Submarine 
Canyon.  These  balanced  runs  seem  more  realistic  than  those  produced  by  the  unbal¬ 
anced  mean  forcing  (runs  1  and  2)  and  depict  broad  anticyclonic  flow  in  the  upper  layer 
which  effectively  fills  the  bay.  These  runs  also  appear  to  agree  better  with  observed  mean 
flows  derived  from  Line  2,  and  also  with  geostrophic  flows  implied  from  dynamic  to¬ 
pographies  (Figure  28  on  page  49).  The  strong  southwesterly  outflow  near  Pt.  Pinos, 
noted  in  the  observed  mean  flow  and  in  satellite  imagery,  appears  to  have  been  handled 
well  by  the  model.  When  the  effects  of  bottom  friction  were  added,  the  majority  of  the 
deep  flow  was  confined  to  the  canyon  topography,  also  supported  by  the  observed  weak 
flow  at  depth  in  the  Line  2  mean.  The  effects  of  friction  are  similar  to  those  noted  by 
Garcia  (1971)  who  observed  that  in  deep  areas  of  the  bay,  friction  played  a  small  role 
but  in  shallower  areas,  friction  significantly  altered  the  flow.  Some  evidence  exists  in  the 
mean  Line  2  data  to  support  outflow  (or  return  flow)  along  the  northern  coast  of  the 
bay  where  the  inflow  diverges  however  not  to  the  extent  of  that  predicted  by  the  model 
from  the  unbalanced  mean  forcing  (runs  1  and  2). 

The  major  limitations  of  these  model  results  include: 

•  The  assumptions  used  in  projecting  the  mean  flow  field  onto  the  model's  open 
boundary.  The  mean  flow  inferred  by  Line  1  data  (west  of  the  open  boundary) 
may  not,  in  reality,  have  extended  all  the  way  into  the  model's  western  boundary, 
as  was  assumed  when  forcing  was  applied.  However,  in  view  of  satellite  imagery 
and  dynamic  height  analyses,  it  is  felt  that  model  flow  was  fairly  representative  of 
the  actual  flow  encountered  over  the  period  of  data  acquisition. 
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•  The  upper  layer  actually  extends  to  the  bottom  over  the  shelf  regions  of  the 
Monterey  Bay  however,  due  to  model  assumptions,  these  areas  must  be  treated  as 
two  layers  to  avoid  surfacing  or  grounding  of  the  layer  interface. 

•  The  observed  flows  indicate  that  a  three-layer  model  is  actually  required,  i.e.  the 
top  layer  to  model  flow  over  the  shelf  area,  a  second  layer  to  model  flow  at  inter¬ 
mediate  depths  of  100  to  500  or  600  m,  and  a  bottom  layer  to  model  flow  in  the 
quiescent  regions  of  the  Monterey  Submarine  Canyon  below  about  600  m. 

Results  of  these  runs  are  in  general  agreement  with  Bruner's  (1988)  observations  that, 
when  lower  layer  velocities  are  weak  and  vertical  shear  is  large,  the  model  effectively 
"decouples"  the  two  layers,  resulting  in  independent  flow  patterns.  Bruner  noted  that, 
for  narrow  inflow  to  the  north,  and  broadened  outflow  to  the  south,  anticyclonic  circu¬ 
lation  ensued  over  the  entire  bay.  In  contrast,  anticyclonic  flow  also  resulted  in  the  runs 
produced  in  this  study  for  broad  inflow  over  the  northern  bay,  and  narrow  outflow  over 
the  southern  bay. 

Finally,  model  runs  which  were  forced  solely  by  variable  wind  stress,  while  obtaining 
similar  mean  flow  patterns  in  the  upper  level,  did  not  correspond  to  observed  mean  flow 
at  depth.  The  wind  forcing  did,  however,  show  decoupling  of  the  two  layers.  Obviously, 
other  factors  besides  wind  stress  were  at  work  here  but  these  experiments  do  imply  that 
wind  forcing  may  play  some  part  in  shaping  mean  flow  patterns,  at  least  at  the  surface, 
during  strong  wind  events  lasting  on  the  order  of  a  week.  Another  interesting  aspect  of 
this  model  result  is  the  fact  that  approximately  4  days  after  the  wind  forcing  was  turned 
off  (Figure  39  on  page  70),  flow  remained  relatively  strong  in  both  layers  implying  th 
"signal"  of  these  wind  events  may  be  seen  in  the  mean  flow  for  a  reasonable  amount  of 
time  after  the  winds  have  abated.  Modeled  current  flow'  appeared  to  respond  rapidly  to 
changes  in  both  wind  direction  and  magnitude,  agreeing  with  Hickey's  (1979)  observa¬ 
tions  that  oceanic  responses  to  changes  in  wind  stress  are  rapid  (O(days)).  However,  as 
noted  in  Chapter  III,  model  runs  using  only  wind  stress  forcing  did  not  achieve  any  ap¬ 
preciable  level  of  stability  and  displayed  highly  variable  height  anomaly  patterns. 

From  modeling  results  it  can  be  seen  that  the  combined  effects  of  winds,  bottom 
friction,  and,  of  course,  boundary  forcing  all  play  key  roles  in  producing  mean  circu¬ 
lation  patterns  within  the  Monterey  Bay.  Strong  winds  appear  to  have  a  significant  ef¬ 
fect  on  surface  flow,  while  light  mcsoscale  land  and  sea  breezes  typically  encountered 
under  weak  high  atmospheric  pressure  gradients  have  little  effect  on  surface  flow.  Bot¬ 
tom  friction  played  a  key  role  in  suppressing  lower  level  flow  out  of  the  Mon^rey  Sub¬ 
marine  Canyon.  It  was  also  shown  that  the  model  was  highly  sensitive  to  mass  balanced 
flow,  yielding  strikingly  different  flow  patterns  from  balanced  and  unbalanced  flows. 
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Figure  41.  Run  2  (left)  and  run  4  (right)  height  anomalies.:  Left  figuics  depict 
Run  2  height  anomalies  (m)  for  the  upper  layer  (top)  and  bottom  layer 
(bottom).  Right  figures  show  Run  4  height  anomalies  for  upper  layer 
(top)  and  bottom  layer  (bottom). 


Finally,  while  the  model  used  in  this  study  could  not  be  forced  by  an  alongshore  current, 
this  component  was  probably  at  woik  during  the  period  of  this  study.  Actual  flow  was 
undoubtably  influenced  by  an  alongshore  geostrophic  jet  and,  as  previously  stated,  sim¬ 
ilarities  to  results  obtained  by  Klinck  (1989)  arc  noted.  This  near  coastal  jet  was  prob¬ 
ably  the  major  generation  mechanism  of  the  observed  mean  flow  patterns. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

Both  prior  to  and  during  the  period  studied  in  this  report,  upwelling  was  occurring 
along  the  entire  central  California  coast.  Both  satellite  imagery  and  the  04  May  1988 
geostrophic  velocity  cross  section  at  Pt.  Sur  supports  this  conclusion.  Within  the 
Monterey  Bay  itself,  the  observed  internal  wave  field  was  quite  active  and  effectively 
masked  the  mean  signal  throughout  the  period.  Previous  studies  suggest  these  internal 
waves  are  quite  common  within  the  bay  and  not  an  anomalous  situation  during  this 
study.  Averaging  the  mass  field  and  measured  ADCP  flows  effectively  filtered  out  this 
"noise."  Because  of  this  high  noise  level,  any  type  of  hydrographic  data  acquired  over 
a  short  time  period  from  the  bay  is  of  little  use.  Only  long  term  averaged  data  can  be 
used  to  infer  any  type  of  semi-permanent  flow  pattern  w'ithin  the  bay.  Also,  because  of 
the  large  amplitudes  of  these  internal  waves,  geostrophic  velocities  derived  from  single 
section  mass  field  observations  will  typically  yield  unrealistically  high  values,  especially 
within  the  thermoclinc. 

During  periods  of  strong  surface  inflow  in  the  northern  part  of  the  bay,  the  flow 
within  the  Monterey  Bay  becomes  a  two-layered  uncoupled  system  consisting  of 
anticyclonic  How  at  the  surface  and  weak,  canyon  trapped,  cyclonic  flow  below  the 
thermocline,  limited  mainly  to  the  deep  waters  over  the  canyon  mouth.  The  main 
mechanism  responsible  for  this  flow  pattern  during  May  1988  was  the  alongshore, 
geostrophic,  equatorward  flow  turning  into  the  northern  half  of  the  bay.  As  this 
equatorward  flow  encountered  the  Monterey  Submarine  Canyon,  the  geostrophic  ad¬ 
justment  process  caused  upwelling  within  the  canyon  which  resulted  in  vortex  tube 
compaction  (anticyclonic  flow)  in  the  upper  layers  and  stretching  of  vortex  tubes  at 
depth  (cyclonic  flow)  below  the  thermocline  (Klinck,  1989).  Numerical  simulations  of 
the  mean  flow  pattern  suggest  that  the  effects  of  bottom  friction  may  act  to  confine 
lower  layer  flows  to  the  canyon  itself.  Other  numerical  simulations  suggested  that  wind 
stress  may  induce  anticyclonic  surface  flow  when  pre-frontal  and  post-frontal  winds  were 
used  to  force  the  model  although  these  experiments  were  inconclusive  due  to  lack  of 
evolution  to  a  steady  state. 

If  surface  flow  within  the  Monterey  Bay  is  predominantly  cyclonic  as  previous  au¬ 
thors  suggest  (Breaker  and  Brocnkow,  1989),  upwelling  periods  nay  coincide  with  p  i- 
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ods  of  flow  reversals  that  have  been  observed  in  previous  studies  (Breaker  and 
Broenkow,  1989).  The  amount  of  time  this  flow  reversal  lasts  after  the  cessation  of  up- 
welling  and  the  alongshore  geostrophic  jet  is  unknown.  The  suggested  scenario  leading 
up  to  the  mean  flow  pattern  encountered  in  this  study  involves: 

•  A  7  to  10  day  period  of  strong  northwesterly  winds  produced  by  high  pressure 
ridging  across  California  behind  a  series  of  atmospheric  cold  fronts. 

•  Wind  stress  induced  Ekman  transport  resulting  in  strong  upwelling  along  the  cen¬ 
tral  California  coast. 

•  Development  of  a  strong  geostrophic,  equatorward,  coastal  jet  along  the  central 
coast. 

•  A  geostrophic  adjustment  process,  such  as  that  described  by  Klinck  (1989),  as  this 
coastal  jet  encounters  the  Monterey  Submarine  Canyon,  resulting  in  the  flow  pat¬ 
terns  described  above. 

Since  winds  are  upwelling  favorable  during  most  of  the  summer  season,  anticyclonic 
circulation  within  the  Monterey  Bay  may  be  more  common  than  previously  believed. 

B.  RECOMMENDATIONS  FOR  FUTURE  STUDIES 
1.  Currents 

The  internal  wave  field  within  the  Monterey  Bay  needs  to  be  resolved.  Once  a 
solid  grasp  of  the  frequency  and  wavenumber  spectrum  is  known,  different  sampling 
schemes  can  more  effectively  cancel  out  their  influence  on  the  mean  field.  In  this  way 
mean  flow  can  be  resolved  with  a  minimum  number  of  samples,  thereby  saving  expensive 
ship  time.  To  accomplish  this,  an  array  of  fast  sampling  thermistor  chains  should  be 
deployed  at  strategic  locations  throughout  the  bay  such  as  on  the  shelf  breaks.  An  al¬ 
ternative  approach  to  averaging  out  the  internal  wave  aliasing  is  simply  through  statis¬ 
tical  averaging.  A  higher  sampling  population  should  result  in  smaller  standard  errors 
and  a  better  resolution  of  the  mean  flow,  although  cost  constraints  could  be  more  of  a 
consideration  in  this  case. 

Because  of  non-linear  processes  at  work  within  the  Monterey  Bay,  current  me¬ 
ters  and  ADCP  velocities  are  probably  of  somewhat  more  value  than  geostrophy  in  de¬ 
termining  flow  within  the  bay.  The  better  resolution  of  the  mean  flow  by  ADCP 
techniques  was  apparent  in  the  smaller  calculated  standard  errors  of  this  data  set. 
Again,  the  longer  the  time  series,  the  better.  Bottom  moored  ADCP  data,  especially 
along  the  shelf  break,  would  yield  valuable  information  concerning  the  cross-canyon 
velocities  and  upwelling,  if  any,  that  is  taking  place  from  canyon  depths.  Seasonal 
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spectral  analyses  of  the  MBA  time  series  data  should  continue  to  reveal  much  about  the 
shelf  processes  that  are  occurring  within  the  Monterey  Bay. 

Finally,  as  with  any  scientific  hypothesis,  results  should  be  reproducible. 
Therefore,  a  similar  study,  conducted  during  or  shortly  after  a  comparable  upwelling 
event,  should  be  undertaken  to  see  if  flow  patterns  similar  to  those  found  in  this  study 
exist.  In  this  way,  flow  "reversals"  within  the  Monterey  Bay  could  perhaps  be  correlated 
with  strong  upwelling  events.  Likewise,  a  study  conducted  during  a  strong  poleward 
flow  event  would  also  be  enlightening.  Current  meter  moorings  placed  on  the  shelf  re¬ 
gions  to  the  north  and  the  south  of  the  Monterey  Bay  during  the  period  of  these  studies 
could  verify  correlations,  if  any,  between  these  alongshore  flows  and  circulation  within 
the  bay. 

2.  Numerical  Studies 

While  many  insights  into  the  behavior  of  flow  within  the  Monterey  Bay  were 
gained  from  the  numerical  model  used  in  this  study,  a  few  added  features  would  greatly 
enhance  its  meaningfulness,  namely: 

•  Future  studies  which  use  actual  flow  to  force  this  model  should  be  from  data  along 
a  line  coincident  with  the  model's  open  boundary.  This  would  negate  the  necessary 
assumptions  made  in  this  study  in  projecting  the  flow  field  onto  the  model  grid. 

•  As  mentioned  earlier,  it  would  appear  that  a  three-level  model  would  be  more  ap¬ 
propriate  in  modeling  the  Monterey  Bay  currents.  This  would  diminish  the  effects 
of  vertical  averaging  and  keep  flows  closer  to  their  "true"  velocities. 

•  The  addition  of  a  realistic  coastline  and  surface  topography  would  enhance  the  ef¬ 
fectiveness  of  wind  forced  flow  experiments  within  the  bay. 

•  Since  local  heating  may  play  a  role  in  circulation  processes  within  the  bay  (Bigelow 
and  Leslie,  1930),  the  addition  of  thermodynamic  equations  to  model  air-sea  heat 
exchanges  would  most  likely  yield  interesting  results. 

•  Expanding  the  model  domain  used  in  this  study  westward,  to  encompass  the  deeper 
seaward  portions  of  the  Monterey  Bay  Submarine  Canvon,  and  northward  and 
southward,  to  model  the  shelf  regions  in  these  areas,  would  provide  more  realistic 
results.  The  additional  capability  of  forcing  the  model  with  northward  or  south¬ 
ward  alongshore  currents,  similar  to  Klinck's  (1989),  would  also  be  of  great  value. 
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APPENDIX  A.  MEAN  GEOSTROPHIC  VELOCITIES 


The  following  pages  contain  the  mean  Gcostrophic  velocities  for  the  section  across 
the  mouth  of  the  Bay  (stations  1  to  10).  Velocities  were  cleaved  by  calculating  the  mean 
mass  field  (mean  temperatures  and  salinities)  from  data  collected  during  the  08-11  May 
19S8  NTS  Student  Cruise,  and  then  using  the  thermal  wind  relation  to  calculate  the  ve¬ 
locities. 
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--  GEOSTRPHIC 

4)  AV  DVNC 

VELOCITY 

5) 

BETWEEN  3 

SUM  DYNC 

TAT  IONS  A  AND  5 

A)  SUM  DYNC  5) 

VELOCITYCCM/S) 

TRANSPORT (M5/S ) 

SUM  OF  TRAN* 

0. 

0.0572 

0.0578 

0.30500 

0.30802 

IA.J7 

0 . 784  77E*04 

0.29556E-05 

25. 

0 , 06*6 

0.06*6 

0.24776 

0.25021 

11  .t* 

0.6701°E*04 

0.21 709E*05 

5C . 

0.0466 

0.046, 

0. 1983S 

0.20058 

10.5? 

0. 6372°E*04 

0.1S007E*05 

75 

0.0407 

0.0617 

0.15371 

0. 155°6 

10.5° 

0  6*S52E*04 

0.86337E404 

100. 

0.0586 

o.o:°3 

0.112°° 

0.11422 

5.85 

0 .25°28E  *04 

0 . 36805E*04 

125. 

o  c:?4 

0  037° 

0  07436 

0.074*4 

2.76 

0 . °588 1E*03 

0. 1086’E*06 

ISO. 

0.0570 

0.0371 

0 . 036°  7 

0.03706 

0.42 

0 . 127°3E*03 

0 . 12 7°5E*05 

175. 

0.0188 

0.0364 

0  00000 

0.00000 

0  00 

o  oooooe*oo 

0 . OOOOOE*QO 

O  TOTAL  TRANSPORT  »  0.29556E*05 

0  VELOCITY  CALCULATED  ASSUMING  A  LEVEL  Oc  NO  MOTION  OF  175.  D8 


DISTANCE  BETWEEN  STATIONS  :  2. <0®  KM 


84 


0 

******  C*EC 

iGTRPMIC  VELOCITY  BU-EEN  C 

TAT I0N3  5  AMD 

6 

0 

:<cb ) 

AV  DVHC  5) 

AV  D»N<  6) 

SUM  OAM( 

5)  GUM  DYMC 

o> 

VELGCI IVICM/3 

)  TRAN3PORT(M3/$) 

SUM  OF  TRANS 

0. 

0  C573 

0.053° 

0  c 1508 

0  °lt.66 

b .  5 1 

0.:3187E‘06 

-0 . 70078E*05 

:s. 

0  0«*°b 

0  0SO7 

0.8572b 

0.85  75: 

1.1° 

-0.80090E*02 

-0.72397E*05 

50 

0  C  16  o 

0.06SC 

0.S07t>3 

0.80-8: 

-3.85 

-0 . C  7°b6E  *06 

-0.71596E*05 

rs. 

C  061? 

0  C6l*» 

0. 

0 . 7© 164 

-5  66 

-0 . 3685  IE*  06 

-0 . 68800E  *05 

ICO. 

o  o:°2 

o.o;«& 

0. 7:1:7 

C  .  71  °°8 

-6.16 

-0.6 1 C78E  *06 

-O.65315E*0S 

i:s . 

O.C37° 

0.0383 

0 . *8 1 °° 

C.c806: 

-7.50 

-O.E072?£*06 

-0.61207E*05 

150. 

c  q;?i 

o.c:?:' 

0  66611 

0.66:15 

-°.26 

-O.588°0E*04 

-0.56133E*05 

175 

o  o:^6 

0 . C3«5 

0.607C5 

0.606°0 

-10. :i 

-0.6:0b5E*06 

-0.50:66E*05 

too. 

0.0*58 

e.0358 

0  5?0*1 

0  5686: 

-10.60 

-0.62b°5£*06 

-0.66037£*OS 

::s . 

0.C3S2 

0.035: 

C. 53681 

o.s::*: 

-lo.c: 

-0  6365 7E *06 

-0.37768E*05 

250 

0.024* 

0.034b 

C.4°°b6 

o.6«?6: 

-10  65 

-0.62716E*06 

-o.316::e*o5 

275 . 

o  o:<i 

0.032° 

0.4650b 

0 . 6t,:s5 

-10.51 

-0.606°CE*06 

-0.:5051E*05 

rcc . 

o.o::6 

0-0333 

C.630°6 

o.6:e°s 

-°  65 

-0.5626CE*04 

-0. J8°82E*05 

::s 

o.c::° 

0-0238 

0.2°74° 

0 . 3°55° 

-° .  03 

-0  5:°01E*06 

-0 . 1 225bE*05 

:so. 

0.02:3 

0  02::' 

0. So458 

0.2b27« 

-8.52 

-0.S0«2oE*04 

-0.80656E*04 

■  ;c 

0  c:ic 

O.C3I7 

0.22:21 

0.33055 

-8.38 

-C  .508 75E*C6 

-0.:°71°E*06 

600 

c.021 1 

C.C307 

0 . 200 t8 

o.:°c*8° 

-8.51 

-C.65:81E*06 

0 . : 1 1 55E*06 

6:5. 

0  o:c6 

0  c:°8 

0  2b°56 

0.2b01° 

-b  .53 

-0.310°:E*C6 

0.66637E*04 

650. 

0  o:«7 

o.o:«i 

o.:2°is 

O.C3825 

-2  80 

-C.16:;ie«06 

0.«75:8E*06 

675. 

c.c:&° 

0  o:es 

0  :o°<i 

o.:o°:i 

-0 .  °6 

0.60062E*0: 

0.11 1801*05 

500. 

O.OC81 

0.0:79 

0.1806° 

C  18073 

1.16 

0. io?:6£*o6 

o.jm:ce*o5 

s:s 

c.c:r:' 

0  0:5° 

0.15:66 

0  is:°o 

:  6  ° 

0. 18°61E*04 

0.100:8E*05 

550 

0.02*4 

0.02*3 

0.  i:s:7 

0 . i:o07 

3.81 

o.:s:o:e*06 

0  8 I 5 15E*06 

575. 

0.0:55 

0.0:57 

0 . 0°885 

0.0°°91 

6  56 

0  :6561E‘C6 

0.561 i:E*06 

cCO. 

0.0:4° 

o.o:s: 

0.07326 

0.0760° 

3  5° 

0. 17l°7E*06 

0 . 3155 1E*06 

0.0:65 

o.c:6? 

0.068-.5 

0.0688° 

:.i: 

0. 1036:E*06 

0.16I56E*06 

t50 

0.022° 

0  c:6: 

c .  c::°6 

0.0:6:: 

1.23 

O.2°°i:E*03 

o.3°°i:e*o3 

6  75. 

O.00C3 

0.0:36 

o.oooco 

O.OGCCO 

0.00 

0.0000CE*00 

0.00000E*00 

0 

TOTAL 

TRANSPORT  *  -0. 

70078E* v5 

c 

VELOCITY  CALCULA; 

XP  ASSUMING 

A  LEVEL  OF  NO 

MOTION  Op  b75 . 

CB 

d: 

1ST/KCC  BETWEEN 

STATIONS  * 

2.60°  KM 

c 

0  :<cp) 

AV  DY»i< 

**  GEOSTRPHIC 

6)  AV  DVH( 

VELOCITY 

7) 

BETWEEN  : 

SUM  DY*K 

>TAT  I< 

6’ 

3N3  6  AND  7 

SUM  pvM  7) 

VELOCITY(CM/S) 

TRANSPORT (M2/S) 

SUM  OF  TRANS 

c 

C .  CS8° 

O-CbOO 

1.3152° 

1.33080 

10.67 

0  12°4°E*C5 

0.74730E*05 

:s 

0.C5C7 

0  0517 

1  .  15b6b 

1  16077 

8  24 

0 . 10a27E*  05 

0.«077IE>0S 

50. 

0  C6$0 

0.0653 

1  .105  76 

1 . 10°C7 

6.-C 

0 . °0085E  *04 

0 . 49844E  *  05 

75 

0  061° 

o.c*:: 

1  .0607° 

l  0637b 

5.75 

0.80°26E*C4 

O.<0aJ6E*05 

100. 

0 . 02  °6 

0  c4o: 

1.018°: 

1  .0315° 

5.16 

0  67187E*04 

o.s:74:e*o5 

1:5 

0  0283 

0  02°0 

0  °?°26 

0 .  °8 1 3  7 

2 .  °0 

0 . 4  68  b6£  *04 

0.36033E*05 

150 

O.O:;:' 

0.027° 

0  .  °4 1 0° 

0  °6334 

3.43 

0.371 17E*04 

0.31227E*05 

175 

0.02*5 

0  02  71 

0 .  °C385 

0  °C44° 

1  34 

C  97 184E*03 

0. 18635E*05 

:oo. 

0  0258 

0  02b3 

0  8b7 37 

0  8b/4 1 

0.03 

-0 .5  7350E*02 

0. 17654E*C5 

:cs 

0.0  25:' 

0  0356 

0  8  315  7, 

0  83112 

-0  C 5 

-0. 17°14E*04 

0. 18336E*05 

:sc 

0  0266 

0  0350 

0  7  °  6  3  7, 

0  7°55o 

-I  5? 

-0  38° l 7E  *04 

0.2001 7E*05 

275 

0  032° 

0.0261 

0  7617° 

0.7605° 

-3  22 

-0.2?231E*04 

o.::°o°e*c5 

200. 

0.0233 

0  0  2  34 

0. 73786 

0 . 73666 

-3  70 

-  0 . 608  7  °E  *  04 

0 .2664  )£  *05 

335. 

0  0238 

0.0235 

0.69652 

0.69309 

-3.81 

-0  28007E*04 

0 . 30729E*05 

25  0 

0  02:: 

0.0216 

0.66172 

0  66052 

-3.21 

-  0 . 35  6  68E  *  04 

0.34520E*05 

2  75 

0.0217 

0.0210 

0  63°6° 

0 . 638°0 

-1.15 

-0 . 68  2 1 2E • 02 

0 . 27097E*05 

6&C 

0.0307 

0  0:01 

0  5° 73 2 

0 . 5°7°5 

0.32 

0. 1 1C30E*04 

0. J7780E*05 

635 

o.o:°a 

o.o:°s 

0.56716 

0.56780 

1  3° 

C.33238E*04 

0 . 36657E*C5 

650. 

0  o:°i 

0  c:°o 

0.5273° 

0  529:b 

1  ftb 

0  3°°44£  *04 

0.26222E*05 

675 

0.0385 

0  0383 

0  5C816 

0 . 50°38 

3  18 

0 . 26325E  *04 

0.21 2C°E*05 

500. 

0.c:78 

C  03  7  7 

0  4 7°c8 

0.681C8 

3.71 

0.41080£*C4 

0.:7705E»05 

5C5. 

0.0:6° 

0.0371 

0 . 4  5 1  °  1 

0.65227 

3  83 

0 . 284 1 S£*04 

0.23S°7E*05 

550. 

0.036  2 

0  C3oS 

0  43501 

0.63b32 

3  35 

0.3101 3E*04 

0. 19?56E*05 

5  75 

0 .  c:s  7 

0  03  bO 

0 . 2°8  75 

0 . 2°?70 

1  84 

0 . 325b  7E  *04 

0.16654E*05 

600 . 

0  o:s: 

0.0356 

0  27204 

0.2727J 

1  26 

0  !?696£*04 

0.14:°CE*C5 

635 

0  C36? 

0.0368 

0.24784 

0.26828 

1.04 

0. 12887E*04 

0. 13529E*05 

650. 

o.c:-.: 

0 . C34J 

0  23216 

0  3335° 

0.82 

0 . 1 4 1 7  °E  *  04 

0.11 140E  *05 

6  75 

0  c: 2b 

0  0325 

0  2°°°4 

0 . 3°°50 

t  .09 

C  1 682 1£  *04 

o.°7::6£»c4 

7C0 

c  o:;i 

0  c::i 

0  37528 

C  “75°° 

1.1° 

0  1 7  b  03E  *  04 

0.802°2E*04 

735 

0  033  7 

0  033  7 

C. 35337 

0  35388 

1  .1° 

0. 1 b  7  ° 1  £ • 04 

0.627°1E*04 

S5 


7S0. 

0.0223 

0.0224 

0.22°58 

0.22016 

1.0b 

0. 14&S8E*C4 

0 .66000E  *  06 

r;s. 

0.021* 

0.0:21 

0.20728 

0.20776 

0 .  *2 

0. 1 10746*04 

0 . 2 11 12E  *  06 

800. 

0 .0217 

0 . C2 1  ° 

0.18525 

0.18566 

0.57 

0.651 QSE  *02 

0.200286*04 

b:s. 

0.0C15 

0.0215 

0.16262 

0.1627* 

0.31 

7.3*4846*02 

0 , 1 25  1  °E*06 

8S0 . 

0.0211 

0.0212 

0. 1*212 

0.16226 

0.22 

0  lb262E*02 

0 .  *5702E*03 

875. 

0.020° 

0.020° 

0.120*° 

0.1209<- 

0,00 

0.245266*02 

0.7924)6*03 

“CO. 

0.0206 

0.0205 

0*10011 

0.10013 

0.03 

0 .  °8°1QE*02 

0. 76888E*03 

*25. 

0.0203 

0.0202 

0.07*56 

0.07*5* 

O.I'T 

0.20227E*03 

0.669*76*03 

*50 . 

0.0200 

0.0200 

0.05*28 

0.C5927 

0.17 

0.222476*03 

0 ,46760E*02 

0  75 . 

0.01*8 

0.01*8 

0.02*27 

0.03*33 

0.13 

0.  I  707  1E*02 

0 . 266 1 5E*03 

1000 

0.01°5 

0.01*5 

0.01*69 

0.01°56 

0.10 

0 . 736 1 *E  *  02 

0 . 736  1  *E*02 

ic:s. 

0.0068 

0.01*2 

0.00000 

0.00000 

0.00 

0 . 00C00E*O0 

0 . OOOOOE*00 

0  TOTAL  TRANSPORT  s  0.74?:0E‘C5 

0  VELOCITY  CALCULATED  ASSUMING  A  LEVEL  Of  NO  MOTION  Of  IOCS.  DB 
DISTANCE  BETWEEN  STATIONS  :  5.932  KM 


0  GEOSTRPMIC  VELOCITY  BETWEEN  STATIONS  7  AND  8 


0  :(DB) 

AV  DYN(  7  ) 

AV  DYN(  8) 

SUM  DYNC  7) 

SUM  DYNt  8) 

VELOC I TY(CM/S ) 

TRANSPORT (M2/S ) 

SUM  OF  TRANS 

0. 

0.0600 

0.0584 

1.0°°81 

l.0*S°8 

-11.12 

-0.8723*6*04 

0 . 1 0 1 S7E  *  06 

25. 

0.0517 

0.06*0 

1.03*78 

1.02752 

-6.56 

-0 . 2*44*E*04 

0. 1 1011E+06 

50. 

0.0652 

0.0627 

0 . *8808 

0 .  *884* 

1  .1* 

0. 35282E  *06 

0. 11275E*06 

75. 

0.0422 

0.0612 

0 . *4277 

0 .  *6482 

5. *5 

0 . 7320CE*06 

0. 10922E*06 

100. 

0.0402 

0.0400 

C.90061 

0.*C2©6 

8.8* 

0.914316*04 

0.10 1 *0E  *06 

125 

o.o:«o 

0.039* 

0.86029 

0.86271 

*.64 

0.*7785E*04 

0.92761E*05 

15C . 

0.0;?* 

0.0290 

0.82125 

0.82485 

10.15 

0 .  *887*6*04 

0.82°82E*05 

175. 

0.0271 

0.0373 

0.78250 

0.7868* 

9.86 

0.92898E*04 

0.73094E*05 

20C  . 

C  .0263 

0.0 2*7 

0.74642 

0.74*56 

9.14 

0.860246*04 

0 . 63705E*05 

225. 

0.0256 

0.02S8 

0.71014 

0.712?5 

7.86 

0.7484*E*04 

0.55302E*05 

25C. 

0.0250 

0.026* 

0.67457 

0.67708 

7.28 

O.72522E*04 

0.47817E*05 

275. 

0.0241 

0.0343 

0.63*60 

0.64215 

7.3* 

0 . 701 1 7E*04 

0.405646*05 

:oc. 

0.0224 

0.0227 

0.60547 

0.60781 

6.7* 

0.622346*04 

0.335536*05 

225. 

0.0225 

O.C22? 

0.5720* 

0.5741b 

6.00 

0.56917E*04 

0.2722*E*05 

350. 

O.C216 

0.021* 

0.52954 

0.54144 

5.52 

0.512286*04 

0.21527E*05 

37S. 

0.0310 

0.0211 

0.507*1 

0.50*58 

4.85 

0.45097E*04 

0.16415E*05 

6CC. 

0.0301 

0.0204 

0.476*6 

0.67843 

6.28 

0.3°092E*04 

0 . 1  1905E*05 

625. 

C .  02*5 

0.02*7 

0.64681 

0.64806 

3.62 

0.2S025E*04 

0 .7**676*04 

660. 

0.02*0 

0.02*1 

0.41727 

0.61822 

3.05 

0.281 08E*04 

0.469426*04 

475. 

0.C282 

0.C2S5 

0.2892* 

0.38*20 

2.64 

0  21792E*04 

0.188346*04 

SCO. 

0  0277 

0.027* 

0.2600* 

0 . 3t>C?0 

1.77 

0. 14765E*04 

-0.295856*03 

525* 

0.0271 

C.C274 

0.23228 

0.22280 

1.22 

0 . ®0058E*03 

-0.17?24E*04 

550. 

0 . 02©S 

0.02bb 

0.30S26 

0.20S4S 

0.61 

0  47046E«03 

-0.26?2«E*04 

575. 

0 .  C2o0 

0.0260 

0.27871 

0.27882 

0.25 

0.357666*03 

-0.31434E*04 

600. 

0.0256 

0 .0252 

0.25274 

0.25287 

0.28 

0 . 38705E*03 

-0.3$01QE*04 

625 . 

0.0248 

C . 0263 

0.2273* 

0.22752 

0  41 

0.3346*E*03 

-0 . 38881E*04 

650. 

0.0241 

0.0262 

0.202*0 

0 .20270 

0.27 

-0.8*030E*02 

-0.42228E*04 

6  75 

0.C235 

0.0227 

0.17851 

0.1782b 

-0  4S 

-0.767806*03 

-0.41 367E*04 

700. 

0.0231 

0.0231 

0.15501 

0.15466 

-1  v7 

-0. 10707E*06 

-0.33890E*04 

725. 

0.C227 

0.0226 

0.1219* 

0.13151 

-1.10 

-0.*2449E*03 

-0.231836*04 

750. 

0.0224 

0.0222 

0.10*15 

0.1098* 

-0  77 

-0 . 65159E*03 

-0. 139385*04 

775. 

0.0221 

0.0220 

0.08677 

0.09658 

-0.55 

-0. 388S0E*03 

-0.742216*03 

800. 

0.021* 

0.0218 

0.064*5 

0.06657 

-0.24 

-0.1*0666*03 

-0.35391E*03 

825. 

0.0215 

0.0215 

0.04280 

0.04275 

-0  15 

-0.1 I804E*03 

-0.163276*03 

850. 

0.0212 

0.0212 

0.02125 

0.02122 

-0  0* 

-0 . 45235E*02 

-0.452356*02 

875. 

0.020* 

0 

0 

r> 

0 

CO 

0.00000 

0.00000 

0  00 

0 . 00000E*00 

0.000006*00 

TOTAL  TRANSPORT  *  0.101276*06 


0  VELOCITY  CALCULATED  ASSUMING  A  LEVEL  OF  NO  MOTION  OF  87S.  DB 
DISTANCE  BETWEEN  ST  AT  10.. S  :  3. *52  KM 


0 

GEOSTRPMIC 

VELOCITY 

BETWEEN  : 

STATIONS  8  AND  9 

0 

2(DP) 

AV  OYN(  8) 

AV  OVN< 

9) 

SUM  DVN( 

8)  SUM  DYN(  *) 

VELOC I TYf CM/S) 

TRANSPORT (M3/S) 

SUM  of  tran: 

0. 

0  0584 

0.0S69 

0.81-MS 

0  .61424 

-  12.06 

-0  585516*04 

0.9S366E*04 

25. 

0.04*0 

0.0482 

0.75970 

0.75743 

-5.65 

-0.246656*04 

0 . 15 1 92E*05 

50. 

0.0427 

0  04  21 

0.70*66 

0.70*22 

-1  86 

-0.3*201E*CJ 

0.m!8E*05 

7$. 

0.C412 

0.04C9 

0.665** 

0  665 15 

0  63 

0.10248E»04 

0. 180S0E*05 

80 


100. 

O.CAOO 

0.03*8 

O.*2403 

0.62528 

2.28 

0. 18674E*04 

0.17005E-C5 

ICS. 

0.028? 

0 . C3®0 

0.58438 

0 . 535*2 

2.20 

0. 1  ?83?E*04 

0.  151 285*05 

150. 

0.0C80 

0  028C 

0.54*02 

0.546*5 

2.72 

0.16*385*04 

0.1315AE*0S 

ITS. 

0.0372 

0.027C 

0.5030* 

0.5C965 

1.65 

0.12794E*04 

0 . 1 16*15*05 

COO. 

0.0267 

0.0263 

0.47076 

0.67124 

2.17 

0.20578E*04 

0. 1041 1E*05 

CCS. 

0.0258 

0.0256 

0 .42402 

0.434*5 

3. *7 

0 . 28?6 IE  »06 

0.85536E*04 

cso. 

0.024? 

0.0250 

0.2*825 

0. 3**24 

6  67 

0.20612E*04 

0.545935*04 

C75 . 

0.024J 

0.0262 

0.26222 

0.26426 

6.47 

0.2*«5*E*04 

0.23«79E*04 

200. 

0.0227 

0.0228 

0. 228*8 

0.23002 

6.67 

0.278*25*04 

-0.597*25*03 

2C5 . 

0.02C7 

0 . 02C° 

0.2*533 

0.29*23 

2.86 

0.221 72E*04 

-0.3286*5*04 

cso. 

0.0319 

0.02C1 

0.C6261 

0.26223 

3 .  OC 

0.172985*04 

-0.57061E*04 

275. 

0.0211 

0.0214 

0.C3075 

0.22125 

2.13 

0. 1 106?E*06 

-0.74439E*04 

4C0. 

0.0206 

0.0207 

0. 1**60 

0. 1**88 

1.17 

0.22530E*C2 

-0.8S508E*04 

ACS  . 

O.OC97 

0.020C 

0.1b*C4 

0.16*1* 

-0  20 

-0 . 72630E*03 

-0.88761E*04 

450. 

0.02*1 

0.02*4 

0.12*4* 

0.12*02 

-2.00 

-0.17234E*04 

-0.813«8E*04 

6  75. 

O.OC85 

0.028* 

0.11037 

0.10*64 

-2.15 

-0.21821E‘04 

-0.641645*06 

500. 

O.OC7? 

0.0C78 

0.08187 

0.08108 

-2.37 

-0 . 2 1526E*04 

-0.42334E*04 

525 . 

0.0C74 

0.0C70 

0.052*7 

0.05326 

-3.06 

-0.1 552*E*06 

-0.20798E*04 

sso. 

0.0C66 

O.OCoC 

O.OCooC 

0.02625 

-1.57 

-0 . 52  724E  *03 

-0.527245*03 

575. 

O.0C60 

0 .0254 

c.ccoco 

0  00000 

0.00 

0.O00O0E*00 

0 . OOOOOE*QO 

TOTAL 

TRANSPORT  s  0. 

933665*04 

0  velocity  CALCULATED  assuming  a  level  cf  no  motion  of  STS .  DB 
DISTANCE  BETWEEN  STATIONS  •  £.678  KM 


o  MMi  SECSTRrMIC  VELOCITY  BETWEEN  STATIONS  *  AND  10 


0  2(C?) 

AV  DVN{  *) 

AV  CVN<  10) 

SUM  DVN(  9) 

SUM  DYMC  10) 

VELOCITY (CM/S) 

TRANSPORT (M3/S) 

SUM  OF  TRANS 

0. 

0.C55® 

0.0525 

0.16819 

0.12**® 

-18.57 

-0  1°6®*E*05 

-0.213®8E*05 

25- 

0.0482' 

0 . 064  7 

0.0*123 

0.08*1® 

-11.11 

-0.95985E*04 

-0. 1 1900E*05 

50. 

0.0421 

0 . 04  1 5 

0.04307 

C. 04147 

-3.50 

-0.2201 1E*04 

-0.2201 1E*04 

75. 

0.06C8 

0.014] 

0  00000 

0.00000 

0.00 

O.OOCOOE-OO 

0.000005*00 

0 

TOTAL 

TRANSPORT  *  -0.2] 

1298E*CS 

0  VELOCITY  CALCULATED  ASSUMING  A  LEVEL  0^  NO  MOTION  OF  7S.  DB 
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The  following  pages  contain  the  mean  Gcostrophic  velocities  for  the  inland  section, 
towards  the  head  of  the  Monterey  Submarine  Canyon  (stations  10  through  17).  Veloci¬ 
ties  were  derived  by  calculating  the  mean  mass  field  (mean  temperatures  and  salinities) 
from  data  collected  during  the  OS-11  May  1988  NPS  Student  Cruise,  and  then  using  the 
thermal  wind  relation  to  calculate  the  velocities. 


DISTANCE  BETWEEN  STATIONS  :  5.256  KH 


l 


.  CEOSTROPMJC  VELOCITY  AND  TRANSPORT  FOR  DATA  OF  THE  CRUISE  MAY  1989 

(UNITS  IN  DYNAMIC  METERS) 


0 

mmm 

«•  C-EOSTRPHIC 

VELOCITY 

BETWEEN 

STATIONS  17  AND  i6 

0  Z(DP) 

AV 

DVNt 

1/)  AV  DNN( 

16) 

SUM  DYN( 

17)  SUM  DVNC  16) 

VELOCITVCCM/S)  TRANSPORT (MI/S) 

SUM  OF  tran: 

0. 

0 

05*1 

0.0554 

0.05611 

0.05525 

-2  04  -0.1 085 1E*04 

-0.1 085 16*04 

25 . 

C 

0 

0C58 

0.045? 

0.00000 

0.00000 

TOTAL 

0.00  0 . 0C000E  *00 

TRANSPORT  *  -0.108516*04 

0 . QOOOOE*QO 

0  VELOCITY  CALCULATED  ASSUMING  A  LEVEL  OF  NO  MOTION  OF  25.  DB 
DISTANCE  BETWEEN  STATIONS  :  4.246  KM 


0 

MU! 

ceosTRpHic 

VELOCITY 

BETWEEN  STATIONS  16  AND  15 

o  r(tp) 

AV  PYN( 

16)  AV  DYN( 

15) 

SUM  DYN(  16) 

SUM  DYM(  15) 

VELOCITY (CM/S) 

TRANSPORT (MJ/S  ) 

SUM  OF  TRAN! 

0. 

C .0554 

0.0544 

0.10101 

0.10024 

-1.55 

-0.4°84J£«05 

-0. 332°8E*02 

25. 

0.0457 

0.0460 

0.04566 

0 . 045*3 

0.75 

0.4t5 1 36*03 

0.465136*03 

5C . 

0.01*6 

0.0415 

0.00000 

0.00000 

0  00 

0.00000E*C0 

0 .00000E*C0 

0 

TOTAL 

TRANSPORT  t  -0.332 

«8E*02 

0  VELOCITY  CALCULATED  ASSUMING  A  LEVEL  OF  NO  MOTION  OF  50.  DB 


D1 

STANCE  BE 

TWEEN  STATIONS  :  4.9*3 

KM 

c 

■  ■  * 

GEOSTPPMIC  VELOCITY 

BETWEEN  STATIONS  15  AND  14 

0  2 ( DB ) 

AV  D>N( 

15)  AV  DYN<  14) 

SuM  DVN(  15) 

SUM  DYN(  14) 

VEL0CITY(CM/3) 

TRANSPORT (M2/S) 

SUM  OF  tran: 

0 

0  0544 

0.0521 

0.43*00 

0 

0.77 

0 . 24220E  *C4 

0.236°5E*05 

25. 

0.0460 

0.0457 

0.28465 

0.28623 

5.22' 

0 . 4605*E  *04 

0.212726*05 

50. 

0  0415 

0.0418 

0  228*7 

0.34041 

6.16 

0  4562*E *04 

0.166666*05 

75. 

0. C2°l 

0.02«7 

0.29727 

0 .2°862 

5.12 

0 . 323656*04 

0. 121026*05 

too. 

0  0284 

0.038* 

0.25800 

0.258°l 

2  89 

0.20456E*C4 

0.88659E*04 

125. 

0  C379 

0.0378 

0.21071 

0.22033 

2.17 

0. 1°4<>4E*04 

0.682026*04 

150. 

0  0272 

0.0273 

0.18183 

0.18258 

2.64 

0 .200476*04 

0.487096*04 

175. 

0.02*7 

0.0270 

0.14457 

0.14523 

2.31 

0. 15496E*04 

0.286626*04 

200. 

0.0264 

0.0266 

0.10784 

0.10827 

1.52 

0. °26 18E*02 

0.151676*04 

225. 

0.0259 

0.02*1 

0.07147 

0.07160 

0.77 

0.251 026*03 

0.390496*03 

250. 

0.0255 

0.0256 

0.03554 

0.02556 

0.10 

0.3°465E*02 

0.394656*02 

275. 

0.025C 

0.0349 

0.00000 

0.00000 

0.00 

0.000006*00 

0.000006*00 

0 

TOTAL 

TRANSPORT  c  0.2 

36«SE*0S 

0  VELOCITY  CALCULATED  ASSUMING  A  LEVEL  OF  NO  MOTION  OF  27$.  DP 
DISTANCE  BETWEEN  STATIONS  :  3.237  KM 
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0 

•  — 

■■  GEOSTRPH1C  VELOCITY 

BETWEEN  STATIONS  14  AND  13 

0  2(DB ) 

AV  DYM( 

14)  AV  DYN(  13) 

SUM  DYM(  14) 

SUM  DYNC  13) 

VELOCI Ts (CM/3) 

TRANSPORT (MJ/S ) 

SUM  OF  TRAM! 

0. 

0.0531 

0.0538 

0.420^2 

0-44082 

6.17 

0.35307E*04 

0 .43S5°E*  03 

25. 

O.OnS? 

0.0450 

0.38*13 

o.se*°? 

3.35 

0.237 15E*04 

-0 . 30°52E*04 

50. 

0.0418 

0.0423 

0.34041 

0.34120 

3.04 

0. 158°0E*Q4 

-0.546676*04 

75. 

O.O307 

0.0403 

0.2*862 

0 . 2°8°4 

1 .24 

0. 1J02°E*02 

-0.705566*04 

100. 

0.0186 

0  0387 

0 .258°! 

0.2585° 

-1  .20 

-0.113156*04 

-0.707076*04 

125 . 

0 . 0378 

0.0378 

0 . 22nS3 

0.21  °8S 

-1.85 

-0.1467*6*04 

-0.593?2E*04 

ISO. 

0.0373 

0.0373 

0.18258 

0.18203 

-2.11 

-0.15:836*04 

-0.447136*04 

175. 

0.0370 

0.0368 

0.14523 

0.14470 

-2.04 

-0 . 1 J43°E*  04 

-0.2?33QE*04 

200. 

0.03*6 

0.03*4 

0.10827 

0.10786 

-1.58 

-0 . 95443E*  03 

-0 . 158? 1E*04 

r:s. 

0 . 03b  1 

0.03*0 

0.0716° 

0.07144 

-0  °° 

-0.500  78E*03 

-0.63466E*03 

250. 

0.0356 

0.0355 

0.03556 

0.03547 

-0.36 

-0. 13388E*03 

-0.133886*03 

275. 

0.034* 

0.0348 

0.00000 

0.00000 

0  00 

0 . 00000E*  00 

0 . OOOOOE*OQ 

0 

TOTAL 

TRANSPORT  «  0 .  4355°E*03 

0  VELOCITY  CALCULATED  ASSUMING  a  LEVEL  OF  no  MOTION  of  275.  D8 
DISTANCE  BETWEEN  STATIONS  :  2.°67  KM 


0  CEOSTRPHIC  VELOCITY  BETWEEN  STATIONS  13  AND  12  m.m 

0  r(DP)  AV  DVM(  13)  AV  DVN(  12)  SUM  D/N(  IS)  SUM  DYNC  12)  VELOCITY(CH/S)  TRANSPORT (M3/3 >  SUM  OF  TRANS. 


0. 

C .  05  38 

0 ,  C55  * 

0.U18? 

0. ]448Q 

7.36 

0 . 60  *8  IE  *04 

0 . °0320E*  04 

25 . 

0  0459 

0.C466 

0.08805 

0 . 08°28 

3.00 

0.235s7E*04 

0 , 2°5  3°E*  04 

55. 

C.C423 

0.0*2? 

0.04225 

0  04267 

1.02 

0  5°72C6*03 

0 . 5°720E*03 

75. 

0.0403 

0.0301 

0.00000 

0.00000 

0  00 

0 . CCOOOE  *  00  ' 

0 .  OOOOCE*  00 

TOTAL 

TRANSPORT  s 

0.?0320E*04 

0  VELOCITY  CALCULATED  ASSUMING  A  LEVEL  OF  NO  NOTION  OF  75.  DB 
DISTANCE  BETWEEN  STATION’S  *  4 , 6  °3  KM 


0  GECSTRPMIC  VELOCITY  BETWEEN  STATIONS  12  AND  11 


0  2JD3) 

AV  O'MC  12) 

AV  DYNC  11 )  SUM  OvM<  12) 

S’JM  DYNC  ID 

VELOCI TY(CM/S 

)  TRANSPORT (M3/S) 

SUM  OF  TRANS 

0. 

0 .055* 

O.CS*2  0.1448° 

0-14612 

3.07 

C.26°02£*04 

0.36582E*04 

25 

0.04** 

c.o*.;:'  0 . 08°28 

0.08*°3 

1.62 

0.°4S22E*03 

0 . ?5°°5E*03 

50. 

0.0427 

0.C427  0.042*7 

0.04268 

0.03 

0. 147386*02 

0.147386*02 

75. 

0.C3C1 

0 .0303  0.00000 

0.00000 

0.00 

O.OOOOOE*OD 

0.000006*00 

C 

total 

TRANSPORT  *  0. 

36582E*04 

0  VSL 

OCiTY  CALCULA 

TED  ASSUMING  A  LEVEL  OF  NO  MOTION  OF  75.  DB 

c;: 

;TANCE  BETWEEN 

STATIONS  :  4.603  KM 

C 

GEC 

iSTcpmic  velc 

JCITY  BETWEEN  S 

TAT 

IONS  1 1  AND 

10 

0  21 

!  DP ) 

AV  DYNC  II) 

AV  DVN(  10) 

SUM  DYNC 

U) 

SUM  DYN( 

1C) 

VELOCI TV (CM/S) 

TRANSPORT (K3/S) 

SUM  OF  TRANS 

0. 

0.05*2 

0.C535 

C.  14612' 

0.13°*° 

-15.8° 

-O.I45°OE*CS 

-0.2342?E*05 

25. 

0  04  72 

0.0447 

0 .  C8°®2 

0.0861° 

-°  .23 

-0. 70°°7E*C4 

-0.883866*04 

eC . 

0.0427 

0.0415 

0.04268 

0.04147 

-2.00 

-0. 1 738°E*C4 

-0. 1 7!8°E*04 

75. 

C .0303 

0.0141 

O.OOOUO 

0.00000 

0  00 

0 . 00000E  *00 

0 . 00000E*00 

Q 

TOTAL 

TRANSPORT  *  -0.2 

342*6*05 

c 

velocity  CALCULAT 

ED  ASSUMING 

A  LEVEL  tr  NO 

MOT 

ION  OF  75. 

DB 

DI 

STANCE  BE  TWEE’) 

STATIONS  • 

4.646  KM 

APPENDIX  B.  MEAN  ACOUSTIC  DOPPLER  CURRENT  PROFILER 

(ADCP)  VELOCITIES 


ADCP  STATISTICS 


(Note:  Negative  velocities  are  into  the  Bay, 

positive  velocities  are  out  of  the  Bay.  ) 

GROUP:  1 

MEAN  DISTANCE  FROM  NORTH  (KM):  A.  313 


MEAN 

STANDARD 

STANDARD 

DF,PTH(M) 

VELOCITY(CM/S) 

DEVI ATI0N( CM/S) 

ERROR(CM/S) 

22 

-2.26 

10.  54 

4.  71 

38 

-0.95  7.94 

GROUP:  2 

MEAN  DISTANCE  FROM  NORTH  (KM):  8.620 

3.55 

MEAN 

STANDARD 

STANDARD 

DEPTH(M) 

VELOCITY(CM/S) 

DEVIATION  CM/S) 

ERROR ( CM/ S) 

22 

-4.48 

10.  99 

4.92 

38 

-3.02 

10.  84 

4.  85 

54 

-1.29  7.84 

GROUP:  3 

MEAN  DISTANCE  FROM  NORTH  (KM):  12.950 

3.50 

MEAN 

STANDARD 

STANDARD 

DEPTH(M) 

VELOCITY(CM/S) 

DEVIATION  CM/S) 

ERROR(CM/S) 

22 

-9.46 

10.  67 

4.  77 

38 

-6.40 

9.  40 

4.20 

54 

-5.94 

9.  57 

4.  28 

70 

-4.71  5.60 

GROUP:  4 

MEAN  DISTANCE  FROM  NORTH  (KM):  17. 324 

2.50 

MEAN 

STANDARD 

STANDARD 

DEPTH(M) 

VELOCITY(CM/S) 

DEVIATION  CM/S) 

ERROR(CM/S) 

22 

-14.97 

3.  86 

1.  72 

38 

-13. 22 

3.  22 

1.44 

54 

-10. 58 

5.  02 

2.  24 

90 


70 

-8.57 

5.41 

2.42 

86 

-4.  28 

4.  97 

2.  22 

102 

-2.  44 

3.  94 

1.  76 

118 

-0.  86 

5.45 

2.44 

134 

1.  03 

6.  66 

2.:  98 

150 

1.  64 

7.  57 

3.  39 

166 

1.  94 

6.  89 

3.  08 

182 

2.  10 

7.  11 

3.  18 

198 

1.  78 

6.  70 

3.  00 

214 

1.31 

6.  07 

2.  71 

230 

0.  63 

5.  23 

2.  34 

246 

-0.  81 

4.  70 

2.  10 

262 

0.07 

4.  78 

2.  14 

278 

0.  70 

5.  29 

2.  36 

294 

-0.  10 

4.  67 

2.  09 

310 

-1.  92 

3.  50 

1.  56 

326 

-2.  97 

2.  82 

1,26 

342 

-2.96 

2.  01 

0.  90 

358 

-3.  24 

1.  45 

0.  65 

374 

-4.  27 

1.  40 

0.  63 

390 

-4.52 

1.  99 

0.  89 

406 

-4.  62 

2.  36 

1.  06 

422 

-4.  64 

3.42 

1.53 

GROUP:  5 

MEAN  DISTANCE  FROM 

NORTH  (KM):  21.809 

MEAN 

STANDARD 

STANDARD 

DEPTH(M) 

VELOCITYCCM/S) 

DEVIATION  CM/S) 

ERROR(CM/S) 

22 

-11.47 

3.81 

1.  70 

38 

-9.  81 

3.09 

1.  38 

54 

-8.  36 

4.  50 

2.  01 

70 

-6.42 

2.  36 

1.  05 

86 

-4.  99 

2,87 

1.  28 

102 

-8.  37 

3.  62 

1.  62 

118 

-8.  86 

5.  29 

2.  36 

134 

-8.  47 

6.  76 

3.  02 

150 

-6.  66 

5.  64 

2.52 

166 

-4,24 

3.  23 

1.45 

182 

-3.  19 

3.59 

1.  60 

198 

-1.  99 

3.  66 

1.64 

214 

-1.  17 

2.  85 

1.  27 

230 

-1.48 

2.  54 

1.  14 

246 

-1.  71 

3,15 

1.41 

262 

-2,06 

3.  79 

1.  69 

278 

-0.  63 

3.  45 

1.54 

294 

-0,59 

3.  69 

1.  65 

310 

-1.59 

3.  19 

1.42 

326 

-2.03 

2.42 

1.08 

342 

-2.  11 

2.  16 

0.  97 

358 

-2.  60 

2.  44 

1.  09 

374 

-2.  30 

1.  93 

0.  86 

390 

-2.  03 

2.  56 

1.  14 

406 

-2.  24 

4.  21 

1.  88 

91 


422 


-3.65 


3.  89 


1.  74 


GROUP:  6 

MEAN  DISTANCE  FROM  NORTH  (KM):  25.809 


MEAN 

STANDARD 

STANDARD 

DEPTH(M) 

VELOCITY(CM/S) 

DEVIATION  CM/S) 

ERROR(CM/S) 

22 

-8.  78 

5.  77 

2.58 

38 

-8.  19 

4.  36 

1.  95 

54 

-6.91 

4.  76 

2.  13 

70 

-6.35 

4.  97 

2.  22 

86 

-6.  62 

5.45 

2.44 

102 

-6.49 

4.  20 

1.88 

118 

-3.62 

2.  67 

1.  19 

134 

-2.  34 

3.54 

1.59 

150 

-1.25 

4.87 

2.  18 

166 

-1.31 

6.53 

2.  92 

182 

-0.84 

7.  10 

3.  17 

198 

-0.48 

6.  26 

2.  80 

214 

-0.  74 

6.  84 

3.  06 

230 

0.  10 

7.  37 

3.  30 

246 

-0.  70 

6.  68 

2  99 

262 

-3.  02 

7.26 

3.  25 

278 

-3.  34 

6.  10 

2.  73 

294 

-2.  71 

4.  39 

1.96 

310 

-1.91 

3.  84 

1.  72 

326 

-2.56 

3.31 

1.48 

342 

-2.  69 

2.47 

1.  10 

358 

-2.  76 

3.  62 

1.  62 

374 

-2.56 

4.  07 

1.  82 

390 

-2.  61 

4.  28 

1.91 

406 

-2.  10 

4.  06 

1.  82 

422 

-1.  14 

3.  19 

1.42 

GROUP:  7 

MEAN  DISTANCE  FROM 

NORTH  (KM):  29.986 

MEAN 

STANDARD 

STANDARD 

DEPTH(M) 

VELOCITY(CM/S) 

DEVIATION(CM/S) 

ERROR(CM/S) 

22 

5.46 

3.46 

1.55 

38 

5.  07 

6.  89 

3.08 

54 

1.  79 

5.97 

2.  67 

70 

-1.94 

4.  89 

2.  19 

86 

-5.  37 

5.  89 

2.  63 

102 

-10. 12 

7.  91 

3.54 

118 

-13. 71 

9.  17 

4.  10 

134 

-12. 88 

6.  73 

3.  01 

150 

-11.  07 

6.  02 

2.69 

166 

-10. 56 

5.  84 

2.  61 

182 

-10. 25 

4.53 

2.  03 

198 

-8.  94 

4.  05 

1.  81 

214 

-8.69 

3.50 

1.  57 
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230 

-8.  96 

4.  15 

1.  86 

246 

-8.  77 

4.43 

1.  98 

262 

-8.  17 

3.  29 

1.47 

278 

-8.  14 

2.  88 

1.  29 

294 

-8.  59 

2.  34 

1.  05 

310 

-7.  85 

1.  62 

0.  72 

326 

-6.  83 

2.  34 

1.04 

342 

-6.  00 

3.  68 

1.64 

358 

-5.  20 

4.  11 

1.  84 

374 

-5.51 

3.  39 

1.52 

390 

-6.  21 

2.  70 

1.  21 

406 

-5.  83 

2.  51 

1.  12 

422 

-5.21  2.87 

GROUP:  8 

MEAN  DISTANCE  FROM  NORTH  (KM):  34.476 

1.28 

MEAN 

STANDARD 

STANDARD 

DEPTH(M) 

VELOCITY(CM/S) 

DEVIATION CM/S) 

ERROR(CM/S) 

22 

10.48 

7.  81 

3.49 

38 

6.  17 

6.  14 

2.  75 

54 

3.  32 

6.  80 

3.04 

70 

3.  92 

7.  81 

3.49 

86 

-1.56 

8.  79 

3.93 

102 

-4.  17 

10.50 

4.  70 

118 

-4.  34 

10.  49 

4.69 

134 

-2.65 

9.28 

4.  15 

150 

-1.89 

9.  24 

4.  13 

166 

-2.09 

9.52 

4.26 

182 

-2.  35 

8.  21 

3.67 

198 

-2.  36 

6.45 

2.89 

214 

-1.  00 

4.  76 

2.  13 

230 

0.  03 

3.53 

1.58 

246 

-0.  98 

3.  66 

1.  64 

262 

-1.  08 

3.  77 

1.  69 

278 

-1.41 

3.  41 

1.53 

294 

-1.  14 

2.  61 

1.  17 

310 

-1.  28 

1.  97 

0.  88 

326 

-1.  10 

1.  37 

0.  61 

342 

-1.44 

2.  25 

1.  00 

358 

-1.  78 

2.  84 

1.  27 

374 

-2.  08 

2.  83 

1.  27 

390 

-1.  89 

3.  13 

1.40 

406 

-1.  04 

3.  49 

1.  56 

422 

-0.  56 

GROUP:  9 

MEAN  DISTANCE  FROM 

3.  36 

NORTH  (KM):  38.597 

1.  50 

MEAN 

STANDARD 

STANDARD 

DEPTH(M) 

VELOCITY(CM/S) 

DEVIATION  CM/S) 

ERROR(CM/S) 

22 

7.  02 

7.  88 

3.52 
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38 

-0.  32 

6.  39 

2.  86 

54 

-0.  98 

5.  19 

2.  32 

70 

-3.  28 

7.51 

3.  36 
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APPENDIX  C.  MONTEREY  BAY  MODEL  BOUNDARY  CONDITIONS 


MEAN  (UNCORRECTED)  VELOCITIES  ALONG  MOUTH  SECTION 


POINT 

DIST(KM)  FROM  SOUTH 

VI 

V2 

DEPTH(M) 

0. 33573E+02 

0. 16357E+01 

0. OOOOOE+OO 

72 

2 

0.  29217E+02 

0. 14175E+01 

0. 00000E+00 

91 

3 

0. 24861E+02 

0. 12107E+02 

0. 00000E+00 

165 

4 

0. 21431E+02 

0. 12209E+02 

0, 39263E+01 

483 

5 

0. 19022E+02 

0. 12831E+01 

-0. 46778E+01 

971 

6 

0. 14852E+02 

0. 84049E+01 

0.  78001E+00 

1528 

7 

0.  99097E+01 

-0. 54945E+01 

0.  34386E+01 

1501 

8 

0. 65941E+01 

-0. 64503E+01 

0.  12990E+01 

940 

9 

0. 26277E+01 

-0. 11061E+02 

0.  OOOOOE+OO 
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MEAN  V  VELOCITIES  PROJECTED  (INTERPOLATED)  ON  MODEL  GRID 

GRID 

POINT 

DIST(KM)  FROM  SOUTH 

VI 

V2 

MODEL  DEPTH(M) 

2 

0. 28500E+02 

0. 14175E+01 

0.  OOOOOE+OO 

90 

3 

0.  21000E+02 

0. 12107E+02 

0.  OOOOOE+OO 

162 

4 

0.  19000E+02 

0. 12209E+02 

0.  39263E+01 

450 

5 

0. 12500E+02 

0. 12831E+01 

-0.  46778E+01 

900 

6 

0. 11500E+02 

0. 84049E+01 

0.,  78001E+00 

1080 

7 

0.  75000E+01 

-0. 54945E+01 

0.  34386E+01 

900 

8 

0. 65000E+01 

-0. 64503E+01 

0.  12990E+01 

540 

9 

0. 20000E+01 

-0. 11061E+02 

0.  OOOOOE+OO 

157 

POINT 

MODEL  BOUNDARY  U  (INTERPOLATED) 
DIST(KM)  FROM  SOUTH  VI 

VELOCITIES 

V2 

65 

0. 32000E+02 

0.59060E+00 

0. OOOOOE+OO 

64 

0. 31500E+02 

0. 70873E+00 

0. OOOOOE+OO 

63 

0. 31000E+02 

0. 82685E+00 

0. OOOOOE+OO 

62 

0. 30500E+02 

0. 94497E+00 

0. OOOOOE+OO 

61 

0. 30000E+02 

0. 10631E+01 

0. OOOOOE+OO 

60 

0. 29500E+02 

0. 11812E+01 

0. OOOOOE+OO 

59 

0. 29000E+02 

0.  12993E+01 

0. 00000E400 

58 

0. 28500E+02 

0. 14175E+01 

0. OOOOOE+OO 

57 

0, 28000E+02 

0. 21301E+01 

0. OOOOOE+OO 

56 

0. 27500E+02 

0. 28427E+01 

0. OOOOOE+OO 

55 

0. 27000E+02 

0..  35554E+01 

0. OOOOOE+OO 

54 

0. 26500E+02 

0. 42680E+01 

0. OOOOOE+OO 

53 

0.,  26000E+02 

0. 49806E+01 

0. OOOOOE+OO 

52 

0. 25500E+02 

0. 56933E+01 

0. OOOOOE+OO 

51 

0. 25000E+02 

0. 64059E+01 

0. OOOOOE+OO 

50 

0..  24500E+02 

0. 71185E+01 

0. OOOOOE+OO 

49 

0.;  24000E+02 

0.  78312E+01 

0. OOOOOE+OO 

48 

0.,  23500E+02 

0. 85438E+01 

0. OOOOOE+OO 

47 

0. 23000E+02 

0. 92565E+01 

0. OOOOOE+OO 

46 

0..  22500E+02 

0. 99691E+01 

0. OOOOOE+OO 
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45 

0. 22000E+02 

0.  10682E+02 

44 

0. 21500E+02 

0.  11394E+02 

43 

0. 21000E+02 

0.  12107E+02 

42 

0. 20500E+02 

0.  12132E+02 

41 

0. 20000E+02 

0.  12158E+02 

40 

0. 19500E+02 

0. 12183E+02 

39 

0. 19000E+02 

0. 12209E+02 

38 

0. 18500E+02 

0. 11368E+02 

37 

0. 18000E+02 

0. 10528E+02 

36 

0. 17500E+02 

0. 96874E+01 

35 

0. 17000E+02 

0. 88470E+01 

34 

0. 16500E+02 

0. 80065E+01 

33 

0. 16000E+02 

0. 71661E+01 

32 

0. 15500E+02 

0. 63257E+01 

31 

0. 15000E+02 

0. 54853E+01 

30 

0. 14500E+02 

0. 46448E+01 

29 

0. 14000E+02 

0. 38044E+01 

28 

0. 13500E+02 

0. 29640E+01 

27 

0. 13000E+02 

0. 21236E+01 

26 

0. 12500E+02 

0. 12831E+01 

25 

0. 12000E+02 

0. 48440E+01 

24 

0. 11500E+02 

0. 84049E+01 

23 

0. 11000E+02 

0. 66675E+01 

22 

0. 10500E+02 

0. 49301E+01 

21 

0. 10000E+02 

0. 31926E+01 

20 

0. 95000E+01 

0. 14552E+01 

19 

0. 90000E+01 

-0. 28223E+00 

18 

0. 85000E+01 

-0. 20197E+01 

17 

0. 80000E+01 

-0. 37571E+01 

16 

0. 75000E+01 

-0. 54945E+01 

15 

0. 70000E+01 

-0. 59724E+01 

14 

0. 65000E+01 

-0. 64503E+01 

13 

0. 60000E+01 

-0. 69626E+01 

12 

0.  55000E+01 

-0. 74750E+01 

11 

0. 50000E+01 

-0. 79873E+01 

10 

0. 45000E+01 

-0. 84996E+01 

9 

0. 40000E+01 

-0. 90119E+01 

8 

0.  35000E+01 

-0. 95243E+01 

7 

0.  30000E+01 

-0. 10037E+02 

6 

0. 25000E+01 

-0. 10549E+02 

5 

0. 20000E+01 

-0. 11061E+Q2 

4 

0. 15000E+01 

-0. 82959E+01 

3 

0. 10000E+01 

-0. 55306E+01 

2 

0.  50000E+00 

-0. 27653E+01 

1 

0. OOOOOE+OO 

0. OOOOOE+OO 

0.  00000E+00 
0.  00000E+00 
0.  OOOOOE+OO 
0.  98157E+00 
0.  19631E+01 
0.  29447E+01 
0.  39263E+01 
0.  32644E+01 
0.  26026E+01 
0.  19407E+01 
0.  12789E+01 
0.  61703E+00 
-0.  44824E-01 
-0.  70668E+00 
-0.  13685E+01 
-0. 20304E+01 
-0. 26922E+01 
-0. 33541E+01 
-0. 40159E+01 
-0. 46778E+01 
-0. 19489E+01 
0. 78001E+00 
0. 11123E+01 
0. 14446E+01 
0. 17770E+01 
0. 21093E+01 
0. 24416E+01 
0. 27739E+01 
0. 31063E+01 
0. 34386E+01 
0. 23688E+01 
0. 12990E+01 
0. 11547E+01 
0. 10103E+01 
0. 86600E+00 
0. 72167E+00 
0. 57733E+00 
0. 43300E+00 
0. 28867E+00 
0. 14433E+00 
0. 00000E+00 
0. OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
0. OOOOOE+OO 
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DERIVATION  OF  MASS  BALANCED  FLOW  FIELD 


RAW  VALUES  MASS  BALANCED  VALUES 

(M**2/S)  (M**2/S) 


GRID  POINT 

VI 

V2 

VI 

V2 

1 

0.0 

0.0 

0.  0 

0.0 

3 

0.5 

0.  0 

0.  3 

0.0 

5 

0.  6 

0.  0 

0.4 

0.  0 

7 

0.7 

0.0 

0.  5 

0.  0 

9 

1.4 

0.0 

0.  9 

0.  0 

11 

2.  1 

0.  0 

1.4 

0.  0 

13 

2.  8 

0.  0 

1.  8 

0.  0 

15 

3.6 

0.  0 

2.4 

0.  0 

17 

4.3 

0.  0 

2.  8 

0.  0 

19 

5.0 

0.  0 

3.  3 

0.  0 

21 

5.  7 

0.  c 

3.  7 

0.  0 

23 

6.  1 

0.  0 

4.  0 

0.  0 

25 

6.  1 

4.  9 

4.  0 

3.  6 

27 

5.  7 

11.  0 

3.  7 

8.  1 

29 

4.  8 

9.6 

3.  1 

7.  1 

31 

4.  0 

3.  9 

2.  6 

2.9 

33 

3.  2 

-5.4 

2.  1 

-7.9 

35 

2.  3 

-16.  0 

1.5 

-23.4 

37 

1.5 

-26.  0 

1.  0 

-38.0 

39 

0.6 

-40.  0 

0.4 

-58.4 

41 

4.2 

7.4 

2.8 

5.5 

43 

2.5 

16.0 

1.6 

11.  8 

45 

0.  7 

26.  0 

0.  5 

19.  2 

47 

-1.0 

34.  0 

-1.5 

25.  2 

49 

-2.  7 

40.  0 

-4.  2 

29.  6 

51 

-3.  2 

11.  0 

-4.  9 

8.  1 

53 

-3.  7 

5.  1 

-5.  7 

3.  8 

55 

-4.  2 

2.  5 

-6.5 

1.  8 

57 

-4.  8 

0.  9 

-7.4 

0.  6 

59 

-5.  3 

0.  3 

-8.  2 

0.  2 

61 

-4.  1 

0.  0 

-6.  3 

0.  0 

63 

0.  0 

0.  0 

0.  0 

0.  0 

NET 

TRANSPORT 

39.4 

85.  1 

0.  0 

0.  0 

LINEAR  CORRECTION  COEFFICIENTS  APPLIED  TO: 

NEGATIVE  VI  VALUES  =  1.5450 
POSITIVE  VI  VALUES  =  0. 6550 
NEGATIVE  V2  VALUES  =  1.4605 
POSITIVE  V2  VALUES  =  0. 7399 
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MODEL  BATHYMETRY  ALONG  OPEN  BOUNDARY 


DISTANCE(KM)  FROM  SOUTH  MODEL  DEPTH  (M) 


31.5 

90.  0 

31.0 

90.  0 

30.5 

90.  0 

30.  0 

90.  0 

29.5 

90.  0 

29.  0 

90.  0 

28.5 

90.  0 

28.  0 

90.  0 

27.5 

90.  0 

27.  0 

90.  0 

26.5 

90.  0 

26.  0 

90.0 

25.5 

90.0 

25.0 

90.  0 

24.5 

90.0 

24.0 

94.5 

23.5 

99.  0 

23.  0 

103.5 

22.5 

108.  0 

22.  0 

126.0 

21.  5 

144.  0 

21.  0 

162.  0 

20.5 

180.  0 

20.  0 

270.  0 

19.  5 

360.  0 

19.  0 

450.  0 

18.5 

540.  0 

18.  0 

585.0 

17.5 

630.  0 

17.  0 

675.0 

16.5 

720.  0 

16.  0 

720.0 

15.5 

720.0 

15.0 

720.  0 

14.5 

720.0 

14.  0 

765.0 

13.5 

810.  0 

13.  0 

855.  0 

12.5 

900.  0 

12.  0 

990.  0 

11.  5 

1080. 0 

11.  0 

1170. 0 

10.5 

1260. 0 

10.  0 

1260. 0 

9.5 

1260. 0 

9.  0 

1260. 0 

8.5 

1260.0 

8.  0 

1080.0 

7.5 

900.  0 

7.0 

720.0 

6.5 

540.  0 

6.0 

461.2 
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5.5 

382.  5 

5.0 

303.  7 

4.5 

225.  0 

4.  0 

213.  7 

3.  5 

202.5 

3.  0 

191.  2 

2.5 

180.  0 

2.  0 

157.5 

1.5 

135.  0 

1.  0 

112.5 

0.5 

90.  0 
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APPENDIX  D.  MEAN  BRUNT- VAISALA  FREQUENCIES  AND  ROSSBY 

RADII 


STATION:  1  LAT:  36  55.2  LON:  122  5.8 

BOTTOM  DEPTH  =  45.5 

MEAN  BRUNT-VAISALA  FREQ  (CPS)  =  0. 1174E-01 
EXTERNAL  ROSSBY  RADIUS  (KM)  =  241.85 

FIRST  INTERNAL  ROSSBY  RADIUS  (KM)  =  1.  95 


STATION:  2  LAT:  36  52.8  LON:  122  5.9 

BOTTOM  DEPTH  =  75.  2 

MEAN  BRUNT-VAISALA  FREQ  (CPS)  =  0. 1063E-01 
EXTERNAL  ROSSBY  RADIUS  (KM)  =  311. 12 

FIRST  INTERNAL  ROSSBY  RADIUS  (KM)  =  2.92 


STATION:  3  LAT:  36  50.5  LON:  122  5.9 

BOTTOM  DEPTH  =  91.  1 

MEAN  BRUNT-VAISALA  FREQ  (CPS)  =  0. 9942E-02 
EXTERNAL  ROSSBY  RADIUS  (KM)  =  342. 38 

FIRST  INTERNAL  ROSSBY  RADIUS  (KM)  =  3.  30 


STATION:  4  LAT:  36  48.  1  LON:  122  5. 8 

BOTTOM  DEPTH  =  186.  3 

MEAN  BRUNT-VAISALA  FREQ  (CPS)  =  0. 7315E-02 
EXTERNAL  ROSSBY  RADIUS  (KM)  =  489. 65 
FIRST  INTERNAL  ROSSBY  RADIUS  (KM)  =  4.97 


STATION:  5  LAT:  36  46.8  LON:  122  5.8 

BOTTOM  DEPTH  =  669.  8 

MEAN  BRUNT-VAISALA  FREQ  (CPS)  =  0. 4895E-02 
EXTERNAL  ROSSBY  RADIUS  (KM)  =  928. 38 

FIRST  INTERNAL  ROSSBY  RADIUS  (KM)  =  11.96 


STATION:  6  LAT:  36  45.5  LON:  122  5.8 

BOTTOM  DEPTH  =  1023. 6 

MEAN  BRUNT-VAISALA  FREQ  (CPS)  =  0. 4220E-02 
EXTERNAL  ROSSBY  RADIUS  (KM)  =  1147.67 
FIRST  INTERNAL  ROSSBY  RADIUS  (KM)  =  15. 76 
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STATION:  7  LAT:  36  42.3  LON:  122  5.8 

BOTTOM  DEPTH  =  1582. 0 

MEAN  BRUNT-VAISALA  FREQ  (CPS)  =  0. 3597E-02 
EXTERNAL  ROSSBY  RADIUS  (KM)  =  1426.77 
FIRST  INTERNAL  ROSSBY  RADIUS  (KM)  =  20. 76 


STATION:  8  LAT:  36  41.1  LON:  122  3.6 

BOTTOM  DEPTH  =  891.3 

MEAN  BRUNT-VAISALA  FREQ  (CPS)  =  0.4457E-02 
EXTERNAL  ROSSBY  RADIUS  (KM)  =  1070. 94 
FIRST  INTERNAL  ROSSBY  RADIUS  (KM)  =  14.49 


STATION:  9  LAT:  36  40. 3  LON:  122  2.  1 

BOTTOM  DEPTH  =  594.  6 

MEAN  BRUNT-VAISALA  FREQ  (CPS)  =  0. 5098E-02 
EXTERNAL  ROSSBY  RADIUS  (KM)  =  874.71 

FIRST  INTERNAL  ROSSBY  RADIUS  (KM)  =  11. 06 


STATION:  10  LAT:  36  38. 8  LON:  121  59. 1 
BOTTOM  DEPTH  =  83.  2 

MEAN  BRUNT-VAISALA  FREQ  (CPS)  =  0. 1041E-01 
EXTERNAL  ROSSBY  RADIUS  (KM)  =  327. 13 

FIRST  INTERNAL  ROSSBY  RADIUS  (KM)  =  3. 16 


STATION:  11  LAT:  36  41.2  LON:  121  58.2 
BOTTOM  DEPTH  =  93.  1 

MEAN  BRUNT-VAISALA  FREQ  (CPS)  =  0. 1005E-01 
EXTERNAL  ROSSBY  RADIUS  (KM)  =  346. 09 
FIRST  INTERNAL  ROSSBY  RADIUS  (KM)  =  3.41 


STATION:  12  LAT:  36  43.6  LON:  121  57.4 
B OTTO 'l  DEPTH  =  93.  1 

MEAN  BRUNT-VAISALA  FREQ  (CPS)  =  0. 9903E-02 
EXTERNAL  ROSSBY  RADIUS  (KM)  =  346. 09 

FIRST  INTERNAL  ROSSBY  RADIUS  (KM)  =  3.36 


STATION:  13  LAT:  36  46.  1  LON:  121  57.9 
BOTTOM  DEPTH  =  648.  2 

MEAN  BRUNT-VAISALA  FREQ  (CPS)  =  0. 4821E-02 
EXTERNAL  ROSSBY  RADIUS  (KM)  =  913.24 

FIRST  INTERNAL  ROSSBY  RADIUS  (KM)  =  11.40 
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STATION:  14  LAT:  36  47.6  LON:  121  58.6 
BOTTOM  DEPTH  =  297.  5 

MEAN  BRUNT-VAISALA  FREQ  (CPS)  =  0.5949E-02 
EXTERNAL  ROSSBY  RADIUS  (KM)  =  618. 68 

FIRST  INTERNAL  ROSSBY  RADIUS  (KM)  =  6.46 


STATION:  15  LAT:  36  49. 3  LON:  121  59.  1 
BOTTOM  DEPTH  =  428.  2 

MEAN  BRUNT-VAISALA  FREQ  (CPS'*  =  0.  5447E-02 
EXTERNAL  ROSSBY  RADIUS  (KM)  =  742. 25 

FIRST  INTERNAL  ROSSBY  RADIUS  (KM)  =  8.51 


STATION:  16  LAT:  36  51.9  LON:  121  59.9 
BOTTOM  DEPTH  =  59.4 

MEAN  BRUNT-VAISALA  FREQ  (CPS)  =  0. 1177E-01 
EXTERNAL  ROSSBY  RADIUS  (KM)  =  276.58 

FIRST  INTERNAL  ROSSBY  RADIUS  (KM)  =  2.  55 


STATION:  17  LAT:  36  54.1  LON:  122  0.7 

BOTTOM  DEPTH  =  27.  8 

MEAN  BRUNT-VAISALA  FREQ  (CPS)  =  0. 1490E-01 
EXTERNAL  ROSSBY  RADIUS  (KM)  =  189. 10 

FIRST  INTERNAL  ROSSBY  RADIUS  (KM)  =  1.51 
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